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ABSTRACT 


^  , 

Previous  investigations  of  longshore  currents  have  included  various'  simplifying  as- 

I  o<~ 

sumptions  and  restricttons^tneluding*  a  planar  beach,  a  steady  and  depth  uniform  flow, 
spatiallv-invariant  bed  shear  stress  and  turbulent  momentum  exchange,  and  the  exclu¬ 
sion  of  surface  wind  stress.  These  assumptions  are  quantitatively  investigated  by  calcu¬ 
lating  the  relative  importance  of  each  term  in  the  longshore  momentum  balance  with  an 
emphasis  on  the  relative  importance  of  wind  forcing  across  the  barred  nearshore  Wind 
and  wave  forcing  of  longshore  currents  across  a  barred  beach  are  examined  using  both 
a  numerical  model  and  field  measurements.  A  local  momentum  balance  was  measured 
at  various  locations  across  the  surf  zone  during  the  SUPERDUCK  experiment  held  at  rt->  US 
the  U.S.  Army  Corps  of  Engineers,' ^CERC  Field  Research  Facility,  Duck,  X'.C.riU.S.A-.*-' 
in  October  19S6.  A  moveable  sled  was  instrumented  with  pressure,  current,  and  wind 
sensors  to  measure  the  various  terms  in  the  longshore  momentum  equation.  Stability- 
dependent  atmospheric  drag  coefficients  for  the  surf  zone  are  determined  from  wind 
stress  measurements  acquired  just  beyond  the  surf  zone  and  wind  speed  measurements 
acquired  from  an  anemometer  atop  the  9  m  sled  mast.  Breaking  waves  were  visually 
identified  and  electronically  marked  on  the  data  tapes.  Mean  current  velocities  at  three; 


effects  of  barred  topography.  The  predicted  longshore  current  profiles  are  compared  with 
field  data  and  bed  shear  stress  coefficients  are  determined  from  model-fitting.  Significant 
spatial  variability  of  longshore  current  [<7(1.0  ms-1)!]  is  observed  across  the  bar  with 
depth  uniform  longshore  flow.  Wind  force  is  found  to  be  a  first-order  term  along  with 
wave  force  under  certain  wind  and  wave  conditions.  Spatially-dependcnt  bed  shear 
stress  coefficients  relative  to  the  bar  are  determined  from  a  local  momentum  balance. 
Mean  bed  shear  stress  coefficients  were  0.0t.»4  ±  0.0010,  0.002  ±  0.0006,  and  0.001  ± 
0.0003  for  offshore  the  bar.  on  top  and  immediately  before  the  bar.  and  in  the  trough. 
Mean  bed  shear  stress  coefficients  determined  from  model-fitting  and  from  the  momen¬ 
tum  balance  are  in  excellent  agreement.  Longshore  current  velocities  are  modeled  within 
±  20%  of  observed  velocities  across  the  barred  beach  at  SUPERDUCK  with  a  mean 
least  square  error  of  7 - 1 5 °  o . 
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I.  INTRODUCTION 


Previous  investigations  of  longshore  currents  have  utilized  many  simplifying  re¬ 
strictions  and  assumptions.  Major  assumptions  and  restrictions  include  a  planar  beach, 
a  steady  and  depth  uniform  longshore  current,  spatially-invariant  bed  shear  stress  and 
turbulent  mixing,  and  the  exclusion  of  surface  wind  stress.  These  assumptions  are  often 
too  restrictive  to  permit  realistic  results.  This  research  quantitatively  investigates  these 
assumptions  by  measuring  the  relative  importance  of  each  term  in  the  longshore  mo¬ 
mentum  balance  equation  (see  equation  (1)). 

Modern  theories  for  mean  longshore  currents  are  based  on  a  longshore  momentum 
balance  (Bowen,  1969;  Thornton,  1970;  and  Longuet-Higgins,  1970a,  1970b)  with  a  ra¬ 
diation  stress  approach  (Longuet-Higgins  and  Stewart,  1964).  For  a  plane-parallel 
beach,  the  complete  longshore  momentum  equation  is: 


My  |  SSVX 
ct  cx 


b  .  y) 

t,  +  t; 
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term:  1  2  3  4  5 


where  a  right-handed  Cartesian  system  is  adopted  with  x  increasing  offshore  (see 
Figure  1),  and 

term  1  =  temporal  change  of  the  time-averaged  depth-integrated  mean  mo¬ 
mentum  per  unit  area  in  the  longshore  direction  due  to  both  steady  (My) 
and  unsteady  (low  (My), 

term  2  =  cross-shore  gradient  of  the  wave-induced  longshore  momentum  flux 
(also  called  the  wave-induced  radiation  stress  gradient), 
term  3  =  bottom  shear  stress  in  the  longshore  direction, 
term  4  =  surface  wind  stress  in  the  longshore  direction,  and 
term  5  =  cross-shore  gradient  of  the  longshore  momentum  flux  due  to  turbu¬ 
lence. 

Most  investigators  (Longuet-Higgins,  1970a;  Jonsson  et  a].,  1974;  Kraus  and  Sasaki, 
1979;  Thornton  and  Guza,  19S6;  Losada  et  al.,  1987;  Baum  and  Basco,  19S7;  among 
others)  apply  the  previously  mentioned  assumptions  and  restrictions  and  then  reduce  the 
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Figure  1.  Section  and  plan  definition  for  the  SUPERDUCK  experiment  (modi¬ 
fied  from  Birkemeicr  and  Dalrymplc,  1975). 


equation  for  longshore  momentum  to  a  simplified  balance  between  the  wave-induced 
radiation  stress  gradient,  the  bottom  shear  stress,  and  the  turbulence-induced  radiation 
stress  gradient.  The  most  restrictive  of  these  assumptions  is  that  of  a  planar  beach,  since 
barred  beaches  are  more  common  than  planar  beaches  (Wright  et  al.,  1986).  Temporal 
variability  of  longshore  currents,  often  exceeding  ±  100%  of  the  mean  flow,  have  been 
measured  in  the  field  (Dette  and  Fuhrboter,  1974;  Wood  and  Meadows,  1975;  Meadows, 
1976).  The  depiction  of  spatially-in variant  bed  shear  stress  and  turbulent  mixing  across 
the  dynamic  surf  zone,  especially  over  a  barred  beach,  may  be  unrealistic.  Additionally, 
recent  investigations  addressing  combined  currents  and  waves  indicate  that  the  bed  shear 
stress  coefficient,  which  is  directly  proportional  to  bed  shear  stress  itself,  may  be  1 .5  to 
3.5  times  larger  in  magnitude  (Grant  et  al.,  1984)  than  the  invariant  values  universally 
used  by  longshore  current  investigators. 

Longshore  current  forcing  is  studied  utilizing  the  data  sets  acquired  during  the 
SUPERDUCK  experiment.  This  experiment  was  held  at  the  U.S.  Army  Corps  of  Engi¬ 
neers,  Waterways  Experiment  Station's  Coastal  Engineering  Research  Center  (CERC). 
Field  Research  Facility  (FRF)  in  Duck,  XC  in  September  and  October  19S6.  A  unique 
instrumentation  system  designed  specifically  for  this  experiment  was  utilized  to  measure 
a  local  momentum  balance  at  various  locations  across  the  surf  zone  of  this  single-barred 
beach  under  various  wave,  current,  and  wind  conditions.  Instruments  were  mounted  on 
a  mobile  sled  to  measure  the  various  terms  in  the  momentum  balance  of  equation  (1). 
Two  triangular,  differential  pressure  sensor  slope  arrays  were  spaced  approximately  2  m 
apart.  Coupled  with  an  absolute  pressure  reference,  the  slope  arrays  can  be  used  to 
measure  the  gradient  of  the  wave-induced  momentum  flux.  Marsh-McBirney 
electromagnetic  (EM)  current  meters  located  at  the  centroid  of  the  slope  arrays  can  be 
used  to  determine  both  mean  momentum  as  well  as  the  total,  or  combined,  momentum 
flux  due  to  the  wave-induced  and  turbulent  components.  The  difference  between  the 
current  meter  determination  of  the  total  momentum  flux  (SVJt)  and  the  slope  array  de¬ 
termination  of  the  wave-induced  momentum  flux  (S>x)  yields  a  measure  of  the  depth- 
integrated  Reynolds  stress,  or  turbulent  component  of  the  momentum  flux.  The  surface 
stress  can  be  inferred  from  the  wind  velocity  measured  using  a  wind  anemometer 
mounted  atop  the  9  m  sled  mast. 

Three  methods  for  calculating  the  bed  shear  stress  coefficient  are  attempted  - 
measuring  the  momentum  balance,  model-fitting  to  field  data,  and  the  mean  profile 
method.  By  operating  in  areas  where  a  plane-parallel  beach  assumption  is  valid  and  by 
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precisely  orienting  dynamic  measurements  perpendicular  to  the  local  bottom  contours, 
all  the  terms  contributing  to  the  longshore  momentum  can  be  measured.  Having  meas¬ 
ured  all  the  terms  in  the  balance,  the  unspecified  bed  shear  stress  coefficient  is  then  de¬ 
termined.  Theoretical  mean  longshore  current  distributions  generated  by  a  numerical 
model  are  compared  with  field-measured  longshore  current  speeds.  The  bed  shear  stress 
coefficients  in  the  numerical  model  are  then  adjusted  until  the  modeled  and  observed 
current  distributions  coincide  in  a  least  square  sense.  Preliminary  indications  from  an 
on-going  surf  zone  investigation  suggest  that  the  mean  longshore  velocity  profile  in  the 
vertical  appears  to  be  logarithmic  (personal  communication,  R.  T.  Guza,  1987).  Hence, 
the  bed  shear  stress  coefficient  might  also  be  determined  by  the  profile  method  of 
boundary'  layer  physics.  However,  logarithmic  profiles  were  not  found  in  the  SUPER- 
DUCK  data. 

The  sled  was  designed  to  make  transects  anywhere  along  the  shoreline  and  as  far 
offshore  as  the  mast  height  would  allow.  Sled  mobility  was  provided  by  U.S.  Army 
Corps  of  Engineers  equipment.  The  sled  was  tethered  by  a  300  m  chain  to  a  heavy-duty 
forklift  on  the  beach.  The  sled  was  towed  offshore  by  the  Coastal  Research  Amphibious 
Buggy  (CRAB)  to  the  desired  location  beyond  the  surf  zone  and  then  detached  from  the 
sled.  The  forklift  dragged  the  sled  onshore  along  its  pre-planned  transect  of  data  col¬ 
lection  locations. 

A  previous  difficulty  with  directly  measuring  momentum  flux  in  the  field  was  the 
precise  orientation  required  of  the  instrumentation  (Thornton  and  Guza,  19S1).  In  the 
SUPERDUCK  experiment,  precise  orientation  of  the  sled  (and  attached  slope  arrays) 
was  achieved  to  within  0.5°  using  laser  ranging  to  triangulate  on  two  prisms  mounted 
2.4  m  apart  on  a  mast  spreader  above  the  sled.  A  daily  calibration  of  current  meter 
orientation  was  accomplished  by  initially  positioning  the  sled  beyond  the  surf  zone  for 
the  first  data  run  of  a  surf  zone  transect.  Based  on  conservation  of  radiation  stress 
outside  the  surf  zone,  the  current  meters  were  numerically  oriented  so  that  their  respec¬ 
tive  radiation  stress  values  were  identical.  This  orientation  was  then  held  constant  for 
the  remaining  sled  locations  along  that  transect.  This  calibration  procedure  is  essential 
to  obtaining  realistic  physical  results  from  the  data  and  is  a  distinguishing  feature  of  the 
research. 

A  primary  improvement  between  the  SUPERDUCK  experiment  and  previous 
large-scale  nearshore  experiments  was  the  extensive  bathymetric  data.  Dai.y 
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bathymetric  profiles  of  the  entire  experiment  area  were  accurately  measured  with  the 
CRAB  and  the  laser  ranging  system. 

Pertinent  literature  is  reviewed  in  the  following  chapter.  A  numerical  model  is  de¬ 
veloped  in  Chapter  III  to  compare  wind  and  wave  forces  for  various  wind  speeds  and 
directions.  The  SUPERDL'CK  field  experiment  is  described  in  Chapter  IV.  Data  anal¬ 
ysis  of  the  field  experiment  is  extensive  and  therefore  is  addressed  in  separate  chapters. 
Preliminary  data  analysis  is  described  in  Chapter  V  and  includes  an  overview  of  the  ob¬ 
served  mean  current  structure  over  the  bar.  Since  surface  wind  stess  is  of  major  impor¬ 
tance  to  this  research  and  few  wind  stress  studies  have  been  conducted  in  the  surf  zone, 
Chapter  VI  is  entirely  devoted  to  the  parameterization  of  the  surf  zone  surface  wind 
stress  for  this  field  experiment.  Calculation  of  the  longshore  momentum  balance,  de¬ 
termination  of  the  relative  contributions  of  each  term  in  the  balance,  and  calculation  of 
spatially-variable  bed  shear  stress  coefficients  are  addressed  in  Chapter  VII.  A  longshore 
current  model  for  planar  beaches  is  modified  to  account  for  barred  topography  and  is 
applied  to  the  data.  The  model's  longshore  current  results  are  compared  to  observations 
and  bed  shear  stress  coefficients  are  determined  by  model-fitting  in  Chapter  VII 1. 
Summary  and  conclusions  are  presented  in  Chapter  IX.  Because  of  the  large  quantity 
of  data  presented  graphically  for  ea  h  day  of  data  analysis,  only  a  single  typical  day 
(October  16)  is  illustrated  within  the  text.  Data  analysis  figures  for  other  days  are  ref¬ 
erenced  in  the  text  by  figure  number  but  are  collectively  grouped  in  Appendix  A. 

The  data  obtained  from  this  field  experiment  represents  the  first  quantitative  data 
set,  including  directional  wave,  surface  wind  stress,  bathymetry,  and  breaking  wave 
measurements,  for  the  study  of  longshore  current  over  barred  beaches.  In  most  previous 
studies,  wind  forcing  of  longshore  currents  has  not  been  considered  due  to  the  absence 
of  quantitative  data.  This  research  quantitatively  determines  the  role  of  wind  forcing  in 
drising  longshore  currents  from  both  field  measurements  and  a  simple  numerical  model. 
Wind  forcing,  as  well  as  a  variable  bed  shear  stress  coefficient,  can  represent  measurable 
and  reasonably  significant  terms  in  the  momentum  balance  equation,  and  therefore  must 
be  taken  into  account  in  any  quantification  of  longshore  currents. 


II.  LITERATURE  REVIEW 

A.  INTRODUCTION 

One  of  the  earliest  works  on  coastal  hydrodynamics  and  processes  was  that  of 
Johnson  (1919).  He  qualitatively  discussed  wave  set-up,  longshore  currents,  and  rip 
currents  (undertow)  as  hydraulic  currents  due  to  waves.  He  utilized  incident  wave  angle 
as  his  criterion  for  rip  currents  (waves  approaching  perpendicular  to  shore)  and  for 
longshore  currents  (waves  approaching  oblique  to  shore). 

In  World  War  II,  amphibious  landing  craft  had  difficulty  landing  their  troops  in  the 
presence  of  a  longshore  current.  The  current  often  represented  a  greater  hazard  than  the 
breakers  themselves  by  causing  the  craft  to  broach  inside  the  breaker  zone.  The  L’.S. 
Navy's  concern  for  this  effect  prompted  Putnam  et  al.  (1949)  to  conduct  one  of  the  first 
applied  research  studies  on  longshore  currents. 

Since  that  time,  various  theories  have  been  postulated  to  account  for  the  basic 
physics  of  longshore  currents  and  their  predicted  velocity  distribution  across  the  surf 
zone.  These  theories  include  conservation  of  energy,  conservation  of  mass,  empirical 
formulation,  and  conservation  of  momentum  . 

B.  RADIATION  STRESS 

Modern  longshore  current  theory'  is  based  on  conservation  of  momentum.  One  of 
Putnam  et  al.'s  (1949)  approaches  was  the  application  of  conservation  of  momentum. 
The  momentum  flux  of  an  obliquely  incident  wave  was  averaged  over  a  wave  period  to 
give  a  mean  momentum  flux  entering  the  surf  zone.  The  longshore  component  of  the 
change  (due  to  wave  breaking)  in  this  momentum  flux  provided  a  driving  force  for 
longshore  currents.  They  assumed  a  steady-state  longshore  current,  and  therefore  the 
driving  force  was  balanced  by  a  bottom  frictional  force. 

The  momentum  flux  approach  was  formalized  and  the  physics  clarified  when  the 
concept  of  radiation  stress  was  introduced.  Longuet-Higgins  and  Stewart  (1964)  defined 
radiation  stress  as  the  "excess  flow  of  momentum  due  to  the  presence  of  the  waves." 
Elaborations  on  the  concept  of  radiation  stress  are  contained  in  Longuet-Higgins  and 
Stewart  (1960,  1961,  1962,  1963,  1964),  Bowen  (1969),  Longuet-Higgins  (1970a,  1970b, 
1972a,  1972b),  and  Thornton  (1970). 

Theory  was  further  formalized  with  the  derivation  of  time-averaged  and  depth- 
integrated  conservation  of  mass  and  conservation  of  momentum  which  separates  the 
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mean  and  fluctuating  momentum  contributions  (Phillips,  1977;  Nlei,  1983).  This  ap¬ 
proach  has  been  used  extensively  in  the  study  of  wave  set-up  and  set-down,  horizontal 
longshore  current  profiles,  rip  currents,  and  surf  beat  generation.  This  approach  forms 
the  foundation  of  modern  longshore  current  theory,  as  well  as  the  foundation  for  the 
momentum  balance  of  equation  (1)  on  which  this  dissertation  is  based. 

When  there  is  an  oblique  wave  approach  to  the  shoreline,  the  longshore  component 
of  radiation  stress  (or  longshore  momentum  flux  due  to  unsteady  motion)  is  given  by 
(Longuet-Higgins,  1970a): 


J.v* 


pu'v'dz 


where 


)/  =  the  sea  surface. 

h  =  water  depth  as  measured  from  the  still  water  line  (SWL), 
p  =  water  density, 

u',v'  =  fluctuating  (x.  y)  components  of  the  flow  velocity,  and 

the  overbar  denotes  a  time-averaging  operator. 


(2) 


The  basic  concept  is  that  the  longshore  current  velocity  is  related  to  the  change  in  in¬ 
coming  wave  momentum,  which  in  turn  relates  to  the  radiation  stress.  This  change  of 
momentum  flux  in  the  cross-shore  direction  is  proportional  to  the  rate  of  wave  energy 
dissipation,  where  dissipation  is  due  primarily  to  surface  breaking  but  can  include  bot¬ 
tom,  and  internal  (turbulent)  dissipation.  By  using  linear  wave  theory,  radiation  stress 
for  monochromatic  waves  can  be  related  to  the  incident  wave  angle  (a)  at  the  measure¬ 
ment  location  by: 

Syx  —  En  sin  a  cos  a  (3) 


where  E  is  total  wave  energy  density,  a  is  the  wave  incident  angle,  or  angle  between  the 
wave  crest  and  a  line  parallel  to  the  shoreline,  and  n  is  the  ratio  of  wave  group  velocity 
to  wave  phase  velocity  expressed  in  terms  of  wave  number  (k)  as: 


n 


2k  h  \ 
sinh(2/c/z)  J 


(4) 
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Radiation  stress  has  also  been  applied  to  random  waves  (Collins,  1972;  Battjes,  1972a; 
Thornton  and  Guza,  19S6)  and  non-linear  wave  theory-  (James,  1974a,  1974b)  to  solve 
for  the  longshore  current. 

The  conservation  laws  of  mass,  horizontal  momentum  llux,  and  radiation  stress  by 
Phillips  (1977)  and  Mei  (1983)  are  both  widely  applied  and  similar.  Both  develop  the 
equations  for  mean  water  motion  by  integrating  over  depth  and  averaging  over  a  wave 
period.  They  differ  in  their  approach  to  the  mean  flow  due  to  unsteady  motion  (i.e. 
waves).  Separating  the  mean  and  fluctuating  horizontal  velocity  components  («,•)  as: 

Ui  =  Ut  +  Uj  i=  1,2,  (5) 

Mei(19S3)  defines  the  mean  flow  as: 

1  f n 

U/ixj-'j)  =  — — -  u,dz  i  —  1,2  (6) 

>]  +  h  J  ~h 


and 


*  = 0  i-l,2  (7) 

Thus  he  incorporates  the  mean  flow  due  to  unsteady  motion  in  the  mean  velocity  term, 
l),  resulting  in  a  mean  momentum  per  unit  area  defined  as: 

M,  =  Mt  =  pljij  +  h )  (S) 

Phillips  (1977)  utilizes  equation  (5)  and  (6),  however  lie  does  not  include  the  mean 
flow  due  to  unsteady  motion  in  the  term  Ut.  He  evaluates  the  time-averaged  depth- 
integrated  current  flow  by  partitioning  the  integration  of  u  (equation  7)  into  two  sepa¬ 
rate  regions: 


(9) 


where  the  first  term  on  the  right-hand  side  is  zero  by  definition.  This  results  in  a  mean 
momentum  per  unit  area  of 
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(10) 


Mt  =  p  ut(n  +  h)  +  p 


utdz 


i  =  1.2 


Phillips  approach  is  adopted  for  this  research. 

Since  the  momentum  balance  of  equation  ( 1 )  represents  the  basis  for  this  research, 
each  remaining  term  --  surface  wind  stress,  bottom  shear  stress,  turbulent  momentum 
flux,  and  temporal  variability  —  will  now  be  examined. 


C.  SURFACE  WIND  STRESS 

Surface  wind  stress  is  a  shear  stress  imparted  to  the  water  column  when  the  wind 
flows  across  the  w'ater  surface,  For  simplicity,  surface  wind  stress  has  been  historically 
neglected  as  a  forcing  function  for  longshore  currents.  Although  no  physically-based 
model  has  specifically  evaluated  the  wind  stress  contribution  to  longshore  currents,  there 
is  empirical  evidence  suggesting  wind  stress  can  be  important,  and  therefore  should  be 
included  as  a  driving  force  for  longshore  currents.  Nummedal  and  Finley  (1978)  used  a 
stepwise  multiple  regression  procedure  utilizing  surf  zone  environmental  parameters  to 
explain  the  variability  of  longshore  current  velocities.  Although  this  procedure  was 
empirical  and  not  theoretical,  the  longshore  current  component  of  the  wind  velocity  was 
found  to  account  for  more  of  the  observed  current  variance  than  any  single  parameter, 
or  combination  of  parameters,  descriptive  of  the  breakers. 

Fox  and  Davis  (1973)  were  able  to  relate  longshore  currents  with  barometric  pres¬ 
sure  at  Lake  Michigan.  Although  wind  is  certainly  related  to  barometric  pressure,  they 
did  not  expand  their  work  to  directly  deal  with  wind.  Komar  (1976)  mentioned  that 
winds  blowing  in  the  longshore  direction  may  contribute  to  the  observed  flow. 

Birkemeier  and  Dalrymple  (1975)  are  one  of  the  few’  nearshore  modelers  who  did 
incorporate  a  wind  forcing  term  (albeit,  an  underweighted  term)  in  their  longshore  cur¬ 
rent  model.  Notwithstanding  the  different  spatial  scales  involved,  it  is  noted  that  lake 
(e.g.,  Sheng  et  al.,  1978)  and  offshore  modelers  do  consider  wind  stress  whereas 
nearshore  modelers  do  not.  Nearshore  modelers  either  assume  that  wave  stress  domi¬ 
nates  over  surface  wind  stress  or  that  no  wind  is  present. 

Hubertz  (1986)  calculated  the  spectra  of  the  mean  longshore  wind  component  and 
mean  longshore  current  for  a  four-month  period  at  the  FRF  in  Duck,  N  O.  Longshore 
current  spectral  peaks  were  present  at  periods  of  4.6  days,  23.8  hours,  and  12.5  hours. 
Longshore  wind  component  spectral  peaks  were  present  at  5.9  days  and  24.4  hours.  He 
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related  the  two  longest  period  current  peaks  to  wind  forcing  by  synoptic  weather  systems 
and  by  sea  breeze  conditions,  respectively.  The  12.5  hours  peak  was  attributed  to  tides. 

There  are  several  methods  for  calculating  surface  wind  stress  on  the  water  surface. 
These  methods  are: 

•  the  direct  method  (also  called  the  covariance  or  eddy  correlation  method), 

•  the  profile  method, 

•  the  dissipation  method,  and 

•  the  drag  coefficient  method  (also  called  the  bulk  aerodynamic  method). 

Each  method  has  its  own  inherent  advantages,  disadvantages,  and  accuracy.  The  drag 
coefficient  method  was  selected  for  this  research  based  on  equipment  availability,  the 
rough  environment  and  platform  involved  (surf  zone  measurements  by  a  mobile  sled), 
and  the  availability  of  wind  stress  measurements  acquired  by  other  field  investigators  at 
the  SUPERDUCK  experiment. 

The  drag  coefficient  method  formulation  for  surface  wind  stress  (Taylor,  1916)  in  the 
longshore  direction  is: 

r]-PaCd\lV\lVy  (11) 

where  =  time-averaged  surface  wind  stress  vector  in  the  longshore  (y)  direction 
(hereafter  referred  to  as  surface  wind  stress), 
pa  -  atmospheric  density, 

Cd  =  an  atmospheric  drag  coefficient  (generally  of  O(10~3)  for  a  height  of  10 
m  above  the  water  surface)  for  the  specific  height  above  the  water  surface 
at  which  W  is  measured, 

W  =  the  total  wind  speed,  and 

lVy  =  the  longshore  component  of  the  wind  velocity. 

The  absolute  value  is  used  in  equation  (11)  to  retain  directionality. 

A  large  quantity  of  research  has  been  directed  at  determining  values  for  Q  for  land 
and  the  open  ocean.  Garratt  (1977)  provides  an  excellent  review  for  investigations  of 
Q  for  land  and  the  open  ocean  prior  to  1977.  More  recent  investigations  include  Wu 
(1980,  19S2),  Safaie  (1984),  Geernaert  and  Katsaros  (1986),  Geernaert  et  al.  (1986, 
1987a,  1987b).  and  Hsu  (1986a,  1986b,  1987).  Unfortunately,  few  investigations  have 
dealt  with  determining  Q  near  the  land-sea  interface  and  in  the  surf  zone.  Q  magnitudes 
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of  0(1.7  x  10-3  to  2.4  x  10~3)  were  determined  for  nearly  calm  conditions  outside  the 
surf  zone  of  a  Chukchi  Sea  barrier  island  (Walters,  1975)  and  Cd  magnitudes  of  0( 
2.0  x  ID-3  to  7.0  x  10-3)  for  the  upper  swash  zone  of  a  Florida  beach  (Hsu.  197U). 
Walter's  (1975)  investigation  was  conducted  under  adiabatic  and  onshore  wind  condi¬ 
tions,  whereas  Hsu's  (1970)  investigation  included  adiabatic  as  well  as  non-adiabatic 
conditions  with  onshore  wind  conditions. 

Hsu  (1973)  determined  surf  zone  wind  stress  values  for  onshore  winds  of  0.5  to  S.5 
ms-1  under  adiabatic  conditions  for  a  two-week  period  at  Santa  Rosa  Island,  FL.  Q 
magnitudes  of  0(1.0  x  10-2  to  1.5  x  10-3)  are  obtained  upon  conversion  of  his  surface 
wind  stress  values.  No  breaker  wave  height  values  were  indicated;  however,  it  is  assumed 
that  they  were  small  due  to  the  experiment's  location  (northwest  Florida  coast),  month 
(June),  wind  speeds,  and  a  photograph  accompanying  the  published  results. 

Yugts  and  Cannemeijer  (19S1)  calculated  Cd  for  non-adiabatic  and  onshore  wind 
conditions  from  field  measurements  at  a  Danish  beach.  They  instrumented  two  20  m 
masts  with  multi-level  cup  anemometers  and  temperature  sensors  and  placed  one  at  the 
high-tide  line  and  one  directly  inland  at  a  distance  of  127  m.  They  concluded  from  ap¬ 
proximately  100  half-hour  profiles  that  the  drag  coefficient  during  unstable  conditions 
over  the  surf  zone  is  the  same  as  over  the  ocean.  They  calculated  Q  magnitudes  of  0( 
1.3  x  10~3  to  1.5  x  10-3)  for  wind  speeds  at  10  m  height  from  3-10  ms-1. 

Geernaert  et  al.  (I9S7)  measured  wind  stress  using  a  sonic  anemometer  from  a  re¬ 
search  platform  in  the  North  Sea.  They  found  an  empirical  relation  between  the  neutral 
drag  coefficient  ( Cdn )  and  the  ratio  cju.,  where  c0  is  the  dominant  long  wave  phase  speed 
and  ».  is  the  surface  wind  friction  velocity.  They  applied  this  relationship  to  other  data 
sets  taken  in  varying  water  depths.  They  found  increased  Cdn  values  for  shallower  depths 
(Figure  2).  As  wind  speed  increased  to  25  ms-’,  the  Cdn  value  for  10  m  depth  was  in¬ 
creased  by  a  factor  of  2  (to  0(3.0  x  10~3))  over  that  for  deep  water. 

D.  BOTTOM  SHEAR  STRESS 

Putnam  and  Johson  (1949  )  conducted  one  of  the  earliest  investigations  of  nearshore 
bottom  stress  or  bottom  friction.  They  assumed  sinusoidal  oscillatory  bottom  motion, 
an  invariant  friction  coefficient,  negligible  perpendicular  flow  at  the  bottom  due  to 
percolation,  a  planar  bottom  slope,  and  a  constant  water  density.  Based  on  their  inves¬ 
tigations,  bottom  shear  stress  (r*)  may  be  represented  by: 

/  =  /uy|  (12) 
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Figure  2.  Calculated  dependence  of  C ^  on  the  10  m  height  Mind  speed  for  five 
depths  (from  Geerneart  et  al,  19S7b). 

where  Cy  =  the  dimensionless  bed  shear  stress  coefficient  due  to  waves  and  currents 
and  generally  of  O(10-2),  and 

ut  =  resultant  vector  sum  of  the  mean  longshore  current  velocity  (V)  and  the 
instantaneous  wave  orbital  speed  near  the  bottom  (see  Figure  3). 

As  with  wind  stress,  the  absolute  value  is  used  in  equation  (12)  to  retain  directionality. 
This  equation  has  essentially  formed  the  basis  for  modern  bottom  friction  studies.  In 
general,  later  investigators  have  worked  toward  its  modification  as  well  as  solving  for 
empirical  values  for  cy. 

One  can  decompose  u  and  v  into  a  steady  or  mean  velocity  (U,  V)  and  an  unsteady 
or  fluctuating  velocity  (u\  v').  The  unsteady  velocity  may  be  further  separated  into 
wave-induced  ( u ,  v)  and  turbulence-induced  (u",  v")  components  (Mei,  19S3): 

u  =  U  +  u'  —  U  +  u  +  u"  (13) 


v  =  V  +  v'  —  V  +  v  +  v" 


(14) 


Incorporating  the  above  decomposition,  the  bottom  shear  stress  is  represented  as  a 
vector  oscillating  in  both  direction  and  magnitude  due  to  the  interaction  of  a  quasi- 
steadv  current  with  superimposed  wave  motion.  From  conservation  of  mass  flux  and  the 
assumption  of  straight-and-parallel  bathymetric  contours,  the  mean  cross-shore  current 
velocity  (L')  is  assumed  zero.  ut  is  then  expressed  as: 

A  AAA 

ut  =  uxi  +  (V  +  uy]j  =  ui  +  yj  —  u  +  v  (15) 


where  ux  =  cos  a, 

Uy  =  uw  sin  a, 

uw  =  wave  orbital  speed  just  above  the  bottom  and  immediately  outside  the 
bottom  boundary  layer,  and 

A 

/ ,  &if  unit  vectors  in  the  x  and  y  directions,  respectively. 

Bottom  shear  stress  (t*)  can  be  written  as: 

"6  b  7  .  b  *  .... 

Z=Txl+TyJ  (16) 

and  from  equations  (12)  and  (15): 

rb  -  pCf\u,\u  +  pcf\ut\v  (17) 

The  longshore  component  of  time-averaged  bottom  shear  stress  can  then  be  written  as: 

t y  =  pcf  I  u,  |  v  =  pcjiii1  +  v2)T v  (18) 

The  time-averaging  procedure  must  be  performed  after  evaluating  the  velocity  compo¬ 
nents. 

The  bed  shear  stress  coefficient  is  a  function  of  bottom  roughness  and  flow  intensity 
near  the  bed.  Bottom  roughness  is  a  function  of  both  sand  characteristics  and  bed  form, 
where  density  and  grain  size  define  sand  characteristics  and  bed  form  is  described  by  the 
existence  of  ripples,  which  are  in  turn,  a  function  of  the  flow  field.  Hence,  the  bed  shear 
stress  coefficient,  bottom  roughness,  and  flow  intensity  are  all  interrelated. 

Often  in  geophysical  research,  laboratory  results  complement  field  measurements 
and  assist  in  building  a  particular  data  base  as  well  as  increasing  our  knowledge  of  the 
phenomena  being  studied.  Unfortunately,  very  few  of  the  large  number  of  laboratory 
studies  of  cf  are  appropriate  for  extrapolation  to  the  field.  The  flow  field  in  nature  will 
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Figure  3.  Vector  Relationship  of  /7,  due  to  presence  of  staves  and  currents. 
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almost  always  be  in  the  rough  turbulent  flow  regime  (Jonsson,  1966).  Since  the 
Reynolds  number  is  the  appropriate  dynamical  similitude  scale  parameter,  large 
Reynolds  numbers  are  required  to  be  in  this  flow  regime.  Hence,  only  laboratory  studies 
performed  in  large  wave  flumes  or  oscillating  water  tunnels,  where  longer  wave  periods 
can  be  obtained,  appear  appropriate. 

Numerous  investigators  have  solved  for  cy  under  waves  (Table  1)  with  results  rang¬ 
ing  from  O(10_1  to  10-3)  for  different  beaches.  Most  have  sought  an  invariant  value  for 
cy  based  on  a  given  set  of  conditions  (e.g.,  bottom  roughness,  beach  slope,  or  for  only 
one  particular  beach),  without  regard  for  any  spatial  or  temporal  variability.  Vitale 
(1979)  states  that  a  common  engineering  practice  is  to  assume  Cy  =  0.01.  Johns  and 
Jefferson  (1980)  showed  the  variation  of  bottom  stress  within  a  wave  cycle  in  a  numer¬ 
ical  analysis.  They  found  increasing  Cy  with  increasing  bottom  roughness  and  decreasing 
wave  period. 

There  are  four  methods  of  bottom  shear  stress,  and  subsequent  Cy ,  determination: 
—  momentum  balance,  model-fitting,  the  profile  method,  and  direct  shear  measurements. 
The  first  three  are  indirect  methods  and  are  addressed  below.  The  direct  shear  meas¬ 
urement  method  involves  directly  measuring  the  shear  force,  such  as  with  pressure  plates 
fixed  on  the  bottom  or  measuring  the  Reynolds  stress  at  the  bed.  This  has  not  been 
done  in  a  surf  zone  and  is  therefore  omitted  from  further  discussion. 

1.  Momentum  Balance 

Forcing  terms  and  assumptions  are  specified  (often  based  on  data  availability) 
to  form  an  appropriate  longshore  momentum  balance  equation.  The  bottom  stress  term 
and  its  accompanying  bed  shear  stress  coefficient  are  then  determined  as  residual  values. 
The  strength  of  this  method  depends  on  how  complete  the  momentum  balance  is  and 
how  realistic  are  the  restrictions  and  assumptions.  Thornton  and  Guza  (19SI)  utilized 
this  method  at  Torrev  Pines  Beach,  CA  (a  near-planar,  unbarred  beach)  in  connection 
with  the  Nearshore  Sediment  Transport  Study  (N’STS)  experiments.  Assuming  steady- 
state  conditions,  neglecting  wind  forcing,  and  operating  in  an  area  of  straight-and- 
parallel  contours,  they  assumed  a  simple  balance  where  the  cross-shore  gradient  of  the 
total  longshore  momentum  flux  (both  wave-induced  and  turbulent  components)  was 
balanced  by  the  mean  longshore  bottom  shear  stress.  For  simplicity,  they  neglected  the 
temporal  variation  term.  They  found  considerable  temporal  and  spatial  variation  of  cf 
with  a  mean  value  of  0.01  +  0.01  .  However  .hey  stated  that  the  results  were  to  be  used 
with  caution  due  to  the  difficulties  in  directly  measuring  the  radiation  stress  ( S)X )  and  its 
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Table  1.  BED  SHEAR  STRESS  COEFFICIENTS  PROPOSED  AND  UTILIZED 
BY  VARIOUS  INVESTIGATORS 


Investigator 

cf 

Comments 

Bretschneider  and  Reid  (1954) 

0.01 

Kishi'  (1954) 

0.03  -  0.04 

Iwagaki  and  Tsuchiya1  (1966) 

0.01  -  0.04 

Iwagaki  and  Kakinuma  (1967) 

0.01  -  1.16 

various  data  sets 

Jonsson  (1966) 

0.01 

Hasselman  and  Collins  (1968) 

0.015 

Longuet-Higgins  (1970a) 

0.01 

James2  (1972) 

0.0025  -  0.001 

Jonsson  et  al.  (1974) 

0.05  -  0.10 

Kamphius  ( 1975) 

0.004  -  0.250 

^bottom  roughness) 

Komar2  (1975,  1976) 

0.017  ±  0.0225 

Huntley  (1976) 

0.0026  +  0.0006 

non-linear  eqn. 

Liu  and  Dalrymple  (1978) 

0.04 

Hsiao  and  Shemdin  (1978) 

0.002  -0.100 

various  data  sets 

Madsen  et  al.2  ( 197S) 

0.01 

Kraus  and  Sasaki  (1979) 

0.011  -0.024 

Vitale  (1979) 

0.003  -  0.250 

Johns  and  Jefferson  (1980) 

0.0096  -  0.00319 

Grant  et  al.  (1984) 

0.005  -  0.01 1 

Wind  and  Vreugdenhil  (1986) 

0.0015 

Guza  et  al.  (198~6b) 

o.oos 

Thornton  and  Guza  (1986) 

0.006  +  0.0007 

•  if  required,  the  investigator's  friction  factor  was  converted  to  a  bed  shear  stress 

coefficient  from  c.-=  (friction 

1  from  Sonu  (1975) 

2  from  Basco  (1982) 

factor  2). 

gradients.  These  difficulties  were  due  to  angularity  errors  caused  by  lack  of  instrument 
resolution,  installation  misalignment,  and  errors  associated  with  choosing  the  correct 
alongshore  direction.  All  of  these  errors  prevented  accurate  directional  resolution  of 
waves  with  small  angles  of  incidence  (<  3°).  Resolution  of  these  difficulties  was  ad¬ 
dressed  in  the  SUPERDUCK  experiment  by  utilizing  a  precise  surveying  system  and 
selecting  days  for  analysis  when  the  mean  wave  incident  angle  was  large  (>  8°).  In  the 
following  determination  of  cj  by  the  momentum  balance  method,  both  the  temporal 
variability  and  wind  stress  terms  in  equation  (1)  are  included. 
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2.  Model-fitting 

Theoretical  horizontal  mean  longshore  current  distributions  generated  by  a  nu¬ 
merical  model  are  compared  with  field-measured  longshore  current  speeds.  The  Cyin  the 
numerical  model  is  then  adjusted  so  that  the  two  profiles  agree.  The  resulting  numerical 
model  values  for  Cy  can  then  be  assumed  to  apply  for  the  field  measurements.  Selecting 
the  appropriate  physics  for  the  numerical  model  is  obviously  important  for  this  method. 
Thornton  and  Guza  (1986)  used  this  method  with  a  narrow-band,  random-wave  trans¬ 
formation  model  for  data  collected  in  another  NSTS  experiment  at  Leadbetter  Beach, 
CA.  A  hierarchy  of  cy  solutions  were  derived  for  linear  and  non-linear  formulations  for 
bottom  shear  stress,  and  with  and  w’ithout  turbulent  momentum  exchange.  Thornton 
and  Guza  (1986)  found  a  least-squares  "best  fit"  value  of  cy  =  0.006  ±  0.0007  using  the 
full  non-linear,  and  presumably  more  correct,  formulation  for  the  bed  shear  stress. 

3.  Profile  Method 

The  profile  method  (e.g.,  see  Roll,  1965;  Schlichting,  1968)  has  been  utilized  in 
boundary  layer  physics  for  years  to  determine  friction  or  shear  velocities  and  roughness 
heights.  It  assumes  a  constant  stress  layer  with  a  mean  logarithmic  vertical  flow  struc¬ 
ture  and  requires  at  least  two  simultaneous  speed  values  at  different  heights  above  the 
same  location. 

Grant  et  al.  (1984)  conducted  high-quality  near-bottom  boundary  layer  meas¬ 
urements  at  a  continental  shelf  location  in  90  m  water  depth.  Even  at  this  substantial 
depth,  they  found  the  near-bottom  flow  field  to  be  composed  of  oscillatory  currents  due 
to  swell  superimposed  on  mean  currents,  a  situation  similar  to  Figure  3.  The  mean, 
near-bottom  ( <2  m)  velocity  profiles  were  found  to  be  logarithmic.  Using  the  profile 
method,  mean  stress  values  were  3  to  7  times  larger  than  expected  from  predictions  of 
mean  current  alone  using  a  typical  smooth-bottom  drag  coefficient,  and  1.5  to  3.5  times 
larger  than  expected  for  predictions  using  a  drag  coefficient  based  on  the  observed  rough 
bottom.  This  significant  increase  in  mean  stress  magnitudes  was  attributed  to  the  pres¬ 
ence  of  the  swell  wave  oscillatory  currents.  These  results  may  apply  to  the  surf  zone,  if 
the  surf  zone's  entire  vertical  extent  of  chaotic  flow  may  be  likened  to  a  fully  developed, 
neutral,  turbulent  boundary  layer  in  the  longshore  direction. 

Analogous  to  surface  wind  stress  (e.g.,  see  Roll,  1965;  Schlichting,  1968),  the 
longshore  bottom  shear  stress  may  be  written  as: 

T_v  =  pv.2  (19) 
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where  v.  is  the  longshore  shear  or  friction  velocity  associated  with  the  mean  flow. 
Combining  equations  (18)  and  (19): 

2  2 

v»  v. 

cf  =  =====  =  ====•  (20) 

(l«;U')  (u2  +  V2)iV 

The  time-averaged  velocity  profile  in  a  neutral  turbulent  boundary'  layer  over  a 
topographically  simple  bottom  is  logarithmic  and  may  be  expressed  as: 

[^]  <J1> 

where  k  is  von  Karman's  constant  and  and  za  is  the  roughness  height.  By  simultaneously 
measuring  current  speed  at  two  different  heights  (V(z,)  and  V(z2)),  two  equations  in  the 
form  of  (21)  may  be  written.  The  two  unknowns,  v.  and  z0,  are  first  solved,  and  then  cf 
is  obtained  from  equation  (20).  Grant  et  al.  (1984)  employed  a  similar  profile  method 
in  their  deep  water  (90  m  depth)  investigation  where  friction  velocity  was  measured  in  a 
downstream  direction  and  u  and  v  were  measured  in  the  streamwise  and  cross-stream 
directions,  respectively. 

E.  TURBULENT  MOMENTUM  FLUX 

The  horizontal  transfer  of  momentum  due  to  turbulent  mixing  processes  is  not  well 
understood.  When  the  classical  analytical  solutions  for  the  horizontal  profile  of  the 
longshore  current  were  initially  solved  using  monochromatic  breaking  waves  (Bowen, 
1969;  Longuet-Higgins,  1970a,  1970b;  Thornton,  1970),  a  velocity  maximum  and  dis¬ 
continuity  was  observed  at  the  breaker  line,  with  zero  longshore  velocity  at  the  shoreline 
as  well  as  beyond  the  breaker  line.  This  occurred  since  there  were  no  driving  forces 
outside  the  sharply  defined  breaker  line  to  generate  a  longshore  current  and  there  was 
maximum  forcing  occurring  at  the  breaker  line.  However,  field  and  laboratory  data 
showed  a  smooth  transition  of  the  longshore  current  from  outside  the  surf  zone  to  a 
maximum  inside  the  breaker  line.  This  observation  was  attributed  to  turbulent  mixing 
(lateral  shear  stress).  Analytical  solutions  for  longshore  currents  on  a  planar  beach  were 
proposed  with  the  assumptions  of  isotropic  turbulence,  a  planar  beach,  monochromatic 
waves,  and  the  parameterization  of  a  lateral  shear  stress  in  terms  of  an  eddy  viscosity 
coefficient  multiplied  by  a  horizontal  velocity  gradient.  Values  for  the  eddy  viscosity 
coefficient  were  obtained  by  model-fitting. 


Whereas  the  bed  shear  stress  coefficient  is  a  physical  quantity  related  to  bottom 
roughness  and  flow  intensity,  the  eddy  viscosity  coefficient  is  a  function  of  the  flow  field 
and  is  a  parameterization  used  to  remove  the  longshore  current  discontinuity  at  the 
breaker  line.  An  alternative  approach  to  a  monochromatic  wave  model  is  to  describe 
breaking  wave  heights  in  terms  of  a  probability  distribution.  This  results  in  a  more  re¬ 
alistic  region  of  wave  breaking  rather  than  a  distinct  breaker  line.  This  spreading  out 
of  the  breaker  line  results  in  a  smooth  horizontal  profile  of  the  longshore  current,  thus 
eliminating  the  need  to  include  an  eddy  viscosity  (Collins,  1970;  Battjes,  1972a). 

Thornton  and  Guza  (1986)  utilized  a  random  wave  model  and  model-fitting  to  de¬ 
termine  values  of  the  bed  shear  stress  coefficient  and  eddy  viscosity  coefficient  for  both 
a  linear  and  non-linear  expression  for  the  bed  shear  stress.  Eddy  viscosity  coefficients  ( 
<0  were  represented  by  an  adjustable  coefficient,  N,  in  the  range  of  0  <  N  <  0.016  , 
expressed  by  (Longuet-Higgins,  1970a,  1970b): 

Z  =  K\x\(gh)2  x  >  0,  (22) 

where  g  is  the  acceleration  due  to  gravity.  Optimal  cf  and  N  coefficients  were  determined 
by  least-squares  fitting  between  calculated  and  measured  longshore  current  speeds 
shoreward  of  the  mean  breaker  line.  Using  a  non-linear  bed  shear  stress  formulation,  a 
Cj  value  of  0.006  was  determined.  Eddy  viscosity  was  found  not  to  be  important  for  the 
near  planar  topography  of  their  field  experiment. 

The  last  term  in  equation  (1)  represents  an  expression  for  horizontal  diffusion  or 
lateral  mixing  of  turbulent  momentum  in  terms  of  a  depth-integrated  Reynolds  stress. 
The  longshore  component  of  momentum  flux  due  to  unsteady  motion  (equation  1)  is  the 
depth-integrated  covariance  between  the  fluctuating  velocity  components  of  u  and  v. 
measured  by  the  current  meters.  By  partitioning  the  unsteady  motion  as  that  due  to 
waves  and  turbulence  (equations  13  and  14),  the  excess  momentum  flux  (equation  2)  can 
be  separated  into  a  wave-induced  ( Syx )  and  turbulent  component  (S"rJC): 

~  p»7 _  P'7 

Syx  *  Syx  +  S"yx  =  L  puvdz  +  I  pu'v’dz  (23) 

It  has  been  assumed  for  this  separation  that  the  characteristic  time  scales  of  the  two 
fluctuations  arc  vastly  different  (i.e.,  ut  and  u",  are  uncorrelated  and  n,-  u"j  =  0).  inte¬ 
grated  viscous  stresses  are  negligible,  and  that  the  horizontal  scales  of  depth  and  mean 
velocity  are  comparable.  The  turbulent  component  of  the  excess  momentum  flux  due 


to  unsteady  motion  may  be  obtained  by  comparing  the  total  radiation  stress  (Syx) 
measured  by  a  current  meter  to  the  wave-induced  radiation  stress  (Svx)  measured  by  a 
slope  array. 

F.  TEMPORAL  VARIABILITY 

Most  previous  investigators  have  assumed  steady-state  conditions  for  longshore 
momentum.  However,  significant  temporal  variations  in  longshore  current  velocity 
measurements  have  been  observed  (Putnam  et  al.,  1949;  Inman  and  Quinn,  1951;  Inman 
and  Bagnold,  1963;  Dette  and  Fuhrboter,  1974;  Dette,  1974;  Huntley  and  Bowen,  1975; 
Wood  and  Meadows,  1975;  Meadows,  1976,  1978;  Guza  and  Thornton,  1978;  Guza  et 
al.,  1986b;  and  others).  The  temporal  variability  of  mean  currents  (dMyISt)  will  be  in¬ 
cluded  in  the  analysis  where  My  is  the  time-averaged  and  depth-integrated  mass  flux  per 
unit  area,  or  momentum  per  unit  area,  in  the  longshore  direction  due  to  steady  {XL)  and 
unsteady  (My)  flow: 

My  =  My  +  My'  =  pV(7i  +  h)+  r  pv'dz  (24) 

Inman  and  Bagnold  (1963)  ascribed  temporal  variability  oflongshore  currents  to  the 
wave  groupness.  Meadows  (1976)  concluded  that  the  total  observed  longshore  current 
velocity  was  composed  of  three  components: 

^  ~  ^  steady  "F  ^  wave  "F  ^  long 

where  Fsteady  is  the  steady  longshore  component,  Fwave  is  the  coincident  wave  period 
fluctuating  component,  and  Llong  is  the  long  period  (as  related  to  surf  beat  or  low  mode 
edge  waves)  fluctuating  component.  He  defined  the  steady  state  component  as  the  ve¬ 
locity  component  whose  period  of  variation  was  greater  than  200  seconds.  The  selection 
of  200  seconds  was  not  explained.  Guza  and  Thornton  (1978)  concluded  that  a  temporal 
averaging  period  for  mean  longshore  currents  was  not  well-understood.  Various  re¬ 
searchers  utilize  different  averaging  times  to  determine  the  steady  state  component, 
based  more  on  record  length  or  data  analysis,  than  on  rationale. 

G.  VERTICAL  DISTRIBUTION  OF  MEAN  LONGSHORE  FLOW 

The  equations  by  Phillips  (1977)  and  Mei  (19S3)  integrate  the  water  motion  over 
depth.  For  simplicity,  the  longshore  current  is  assumed  uniform  from  the  surface  to  the 
bottom.  No  longshore  current  theories  account  for  a  depth-dependent  longshore  cur- 
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rent.  Theoretical  work  on  the  depth-dependency  of  longshore  currents  (Stive  and  de 
Vriend,  1987)  and  cross-shore  currents  (Svendsen  et  al.,  19S7)  has  only  recently  been 
undertaken.  Most  laboratory  and  field  investigations  have  not  measured  the  depth- 
dependence  of  longshore  currents,  and  as  a  result,  there  is  little  field  data  on  vertical 
current  structure. 

Greenwood  and  Sherman  (1983)  measured  vertical  structure  of  the  longshore  cur¬ 
rent  at  the  Wendake  Beach,  Ontario  field  experiment.  They  calculated  mean  flows  for 
three  elevations  as  well  as  a  depth-averaged  mean  longshore  current  velocity.  Their 
maximum  mean  variability  was  only  10%  of  the  depth-averaged  flow.  Thus  their  data, 
as  well  as  Meadows'  (1978)  and  Hubertz's  (1986)  data,  suggests  that  modeling  longshore 
currents  as  uniform  over  depth  is  appropriate.  Although  they  assumed  a  logarithmic 
boundary  layer  velocity  profile  in  order  to  calculate  a  boundary  roughness  length  {z0). 
their  logarithmic  plots  of  only  two  points  arc  inconclusive  regarding  the  existence  of  a 
logarithmic  vertical  velocity  profile  in  the  surf  zone.  However,  their  linear  plots  of  ver¬ 
tical  profiles  show  a  "consistent  deformation  of  the  velocity  field  strongly  suggestive  of 
the  influence  of  bed  friction,"  and  hence  a  logarithmic  profile  within  at  least  0.5  -  1.0  m 
of  the  bed  remains  a  possibility  for  their  data. 

Wright  et  al.  (1986)  obtained  field  data  of  vertical  profiles  from  0.1  to  1.0  m  above 
the  bottom  of  longshore  currents  in  a  bar  trough.  Although  the  longshore  velocities 
were  on  the  same  order  of  magnitude  as  the  wave  orbital  velocities,  the  bottom  bound¬ 
ary  layer  was  dominated  by  the  waves  acting  in  the  shore-normal  direction,  orthogonal 
to  the  longshore  currents.  The  velocity  profiles  were  not  logarithmic  and  it  appeared 
that  they  were  strongly  influenced  by  wave-generated  vortices  from  the  bed  ripples 
present. 

Theoretical  investigations  by  Lundgren  (1972),  Smith  (1977),  and  Grant  and 
Madsen  (1979),  laboratory  studies  by  Bakker  and  Van  Dorn  (197S)  and  Grant  and 
Madsen  (1982),  and  field  measurements  by  Grant  et  al.  (1984)  assume  that  the  mean 
velocity  profile  is  comprised  of  an  inner  region  (close  to  the  bed  which  is  directly  affected 
by  wave-induced  turbulence)  and  an  outer  region  above  it.  and  that  the  outer  region  has 
a  logarithmic  profile  which  is  also  affected  by  wave-induced  turbulence.  Preliminary 
indications  from  an  on-going  surf  zone  field  investigation  suggest  that  the  mean 
longshore  velocity  profile  in  the  vertical  might  be  logarithmic  (personal  communication, 
R.T.  Guza.  1987).  The  existence  of  a  logarithmic  profile  would  allow  bed  shear  stress 
coefficients  to  be  determined  by  the  profile  method. 
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H.  BARRED  BEACHES 


Although  barred  beaches  are  more  common  (Wright  et  al.,  19S6),  far  less  is  known 
about  surf  zone  currents  on  barred  beaches  than  on  planar  beaches.  Allender  and 
Ditmars  (1981)  found  little  spatial  variability  in  the  horizontal  longshore  current  profile 
across  the  entire  surf  zone  of  a  steep,  barred  beach  during  storm  conditions  on  Lake 
Michigan.  Xo  significant  cross-shore  flow  or  vertical  variability  of  current  flow  was 
observed.  Greenwood  and  Sherman  (1983,  1986)  studied  surf  zone  currents  in  Lake 
Huron  off  a  steep,  multi-barred  beach  at  Georgian  Bay.  They  found  significant  cross¬ 
shore  variability  in  the  longshore  current,  and  up  to  a  four-fold  variation  during  storm 
conditions.  These  two  contrasting  conclusions  are  indicative  of  the  embryonic  state  of 
knowledge  concerning  the  surf  zone  dynamics  of  barred  beaches.  The  SL'PERDUCK 
experiment  was  conducted  on  a  barred  beach.  Data  analysis  of  the  surf  zone  currents 
will  be  used  to  improve  the  understanding  of  surf  zone  dynamics  at  a  barred  beach. 
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III.  WIND  FORCING  OF  LONGSHORE  CURRENTS  FROM  A 
FINITE-DEPTH  SPECTRAL  WAVE  MODEL  APPROACH 

To  assess  the  importance  of  the  wind  contribution  to  longshore  currents,  wave  and 
wind  forcing  terms  in  the  longshore  momentum  equation  are  integrated  across  the  surf 
zone  in  a  simple  model  and  compared.  Wave  and  wind  forcing  are  coupled  by  assuming 
the  waves  driving  the  longshore  currents  are  generated  by  the  local  wind.  The  TMA 
wind-wave  shallow  water  wave  spectrum  is  used  as  input  to  the  radiation  stress. 

A.  THE  TMA  SPECTRUM 

The  TMA  spectrum  (Bouws  et  al.,  1985)  was  developed  from  data  obtained  from 
two  extensive  shallow-water  wave  measurement  programs.  MARSEN  (Marine  Remote 
Sensing  Experiment  at  the  North  Sea)  and  ARSLOE  (Atlantic  Ocean  Remote  Sensing 
Land-Ocean  Experiment),  and  data  from  a  prolonged  storm  near  lightship  TEXEL  in  the 
southern  North  Sea.  The  MARSEN  and  ARSLOE  sites  were  both  on  the  continental 
shelf  with  depths  less  than  40  m.  The  combined  data  set  was  named  TMA,  after  the  3 
sub-data  sets,  and  includes  water  depths  from  6  -  42  m,  wind  speeds  up  to  30  ms-1, 
bottom  composition  from  fine  to  coarse  sands,  and  beach  slopes  ranging  from  1:150  to 
nearly  flat.  The  uniqueness  of  the  TMA  data  set  to  this  research  is  that  the  ARSLOE 
data,  which  comprises  two-thirds  of  the  nearly  3000  wind-wave  spectra,  were  measured 
at  the  CERC  Field  Research  Facility  in  Duck,  N.C. 

Bouws  et  al.  (19S5)  proposed,  and  then  verified,  that  a  finite-depth  wind-wave 
spectrum  could  be  approached  through  application  of  the  similarity  principles  by 
Kitaigorodskii  et  al.  (1975).  The  TMA  spectrum  is  an  extension  of  the  self-similar 
JONSWAP  (Joint  North  Sea  Wave  Project)  spectral  shape  (Hasselmann  et  al..  1976)  to 
include  the  effects  of  finite  water  depth  embodied  in  the  Kitaigorodskii  et  al.  (1975) 
finite-depth  transformation  factor,  <PA.(a>,  h)  : 

BtmaW  =  ZpxiifJm*)  'bjifJm’  Vj,  <V  °b)  h)  (26) 


where 

Epm  if)  =  ag2/‘5(2-)-4expC-,-25(f )"]  (27) 
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yj>  G aPb ) 


(28) 


<t>K{o),  h) 


iK<o,h))~3 

{k(a),  oo))-3 


dk(co,h) 


cA(a>,  oo) 


and 


a 

a) 

f 

fm 

Vj 


=  Phillips  constant, 

=  2rS, 

=  frequency, 

=  peak  spectral  frequency, 

=  peak  enhancement  factor,  and 

=  aa  =  peak  spectral  width  on  low  frequency  side  of fm, 
=  ab  —  peak  spectral  width  on  high  frequency  side  of fm. 


(29) 


£pm (ffma)  is  the  Pierson-Moskowitz  (1964)  spectral  shape  for  fully-developed  deepwater 
waves,  and  oaob)  is  the  JOXSWAP  spectral  shape  function. 

Thompson  and  Vincent  (1983)  developed  an  approximation  for  Oa(<u,/i)  using  linear 
theory  and  redefining  the  two  independent  variables  (a)  ,  h)  as  a  single  variable  (o oh): 

1  2  _ , 

—  i  -j 

<DA(tu, /i)  =  0(ft)A)  =  J  l  (30) 

U  _ -1- (2  -  coft)2  «a>i3 


where 


=  2nJ[j-)T  (31) 

However,  the  approximation  (equation  30)  is  incorrect  for  wh  >  2  and  therefore  is  mod¬ 
ified  here  as: 
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<l>A(to,  h) 


1  2 

->  &  h 

“h  ^  1 

1  <  uh  <, 

~ 

1 

(oh>2 

0>(wA)  approaches  a  value  of  one  in  deep  water  and  a  value  of  zero  as  the  depth  de¬ 
creases  (Figure  4). 


Figure  4.  Kitaigorodskii  et  al.  (1975)  finite-depth  spectral  transformation  function, 
4>x(c/,A)  (from  Bouws  et  al.,  1985), 


Bomvs  et  al.  (1985)  developed  the  following  empirical  parameter  expressions  for  in¬ 
put  into  the  TN1A  data  set: 


oi  =  0.007S  K04’ 

(33) 

yj  =  2.47  K0-39 

(34) 

ir?n 

K  = 

(35) 
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a 


(36) 


0-07  /S/J 
0-09  f>fm\ 

where  IV]0  =  wind  speed  at  z  =  10  m  and 

km  =  wave  number  for  waves  at  peak  frequency. 

These  expressions  are  valid  for  all  water  depths  in  the  TMA  data  set.  Thus,  the  TMA 
spectrum  requires  specification  of  Wl0  ,  h,  and  fm. 

B.  NUMERICAL  APPROACH 

To  compare  wind  and  wave  forcing  for  various  conditions  using  finite-depth  spectra, 
it  is  assumed  that; 

1)  Wave  conditions  are  solely  the  result  of  a  local  wind  (i.e.  there  is  no  swell 
present),  and  therefore  the  wind  direction  (6)  and  mean  incident  wave  direction 
are  identical,  and 

2)  Sea  conditions  are  fully  arisen  (i.e.,  the  sea  is  not  duration-limited). 
Comparisons  are  made  for  inputs  of  gentle  and  steep  planar  beach  slopes  (tan  /?  =  .01 
and  .10),  wind  speeds  (10-30  ms-!),  all  wind  directions,  and  peak  spectral  frequencies 
commensurate  with  wind  speed. 

Surface  wind  stress  is  determined  from  equation  (11).  Q  is  determined  from  meas¬ 
urements  taken  by  Geernacrt  et  al.  (19S7bl  in  open  waters  of  the  North  Sea  for  a  depth 
of  10  m  (Figure  2).  These  Q  values  are  increased  by  33%  to  account  for  surf  zone 
roughness.  The  33%  increase  will  be  addressed  later  in  the  Chapter  VI.  Peak  frequency 
is  determined  for  the  Pierson- Moskowitz  spectrum  (Carter,  1982): 

„  1  27 

fm  -  -777-  (37) 

' '  !0 

Peak  frequency  is  assumed  not  to  shift  during  shoaling  (Hughes,  1984). 

After  input  of  tan  /?  and  IF10,  peak  wavenumber  (km)  is  calculated  using  the  linear 
theory7  relationship  starting  at  an  initial  water  depth  of  20  m.  The  TMA  spectrum  is 
determined  from  equation  (26).  The  RMS  wave  height  is  obtained  by  integrating  the 
TMA  spectrum: 


(^rmsVviA  ~  ^  ?,W)  2 

0.0  J 


(38) 
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where  only  sea  waves  (maximium  frequency  of  0.5  Hz)  are  considered. 

Thornton  and  Guza  (1982,  1983)  found  an  envelope  curve  relating  //rms  to  depth  in 
the  inner  surf  zone: 

HTms a  0.42  h  (39) 


The  Hrms  determined  by  equation  (38),  which  includes  the  effects  of  variable  depth,  is 
checked  against  equation  (39)  to  ascertain  if  wave  breaking  due  to  shoaling  had  com¬ 
menced.  If  the  ratio  Hrmslh  is  less  than  required  for  surf  zone  wave  breaking,  water 
depth  is  reduced  and  a  new'  Htms  determined.  This  process  for  calculating  breaker  depth 
iterates  until  equation  (39)  is  satisfied.  Surf  zone  width  (xb)  is  then  determined  from: 


h 

tan  ft 


(40) 


The  radiation  stress  gradient  (8Syx  /  dx)  and  surface  wind  stress  (rff)  are  both  inte¬ 
grated  with  respect  to  xb  to  allow  comparison  of  a  wave  force  (S^)  and  a  wind  force 
(t!J  xb).  with  units  of  force  per  unit  longshore  length.  The  integration  is  confined  to  xb 
since  the  radiation  stress  gradient  is  conserved  outside  the  surf  zone  where  there  is  no 
wave  breaking.  Momentum  exchange  between  the  wind-driven  longshore  current  out¬ 
side  xb  and  the  w'ind-  and  wave-driven  longshore  current  inside  jrfc  is  assumed  negligible. 
Using  equation  (2).  wave  force  is  expressed  as: 


dx  =  —  <4.MA*n(Q  PS  sin  6  cos  6  (41) 

w'here  the  energy  is  expressed  in  terms  of  variance,  which  is: 
o.s 

CTTMAxn(f)  —  X^tmaOO  n(J)Af  (42) 

0.0 

and  n(f)  is  expressed  in  a  frequency-dependent  form.  The  factor  of  1/2  is  an  adjustment 
factor  for  broad-banded  wind-driven  waves  which  incorporates  the  effect  of  short¬ 
crestedness  on  the  ratio  of  radiation  stress  to  the  average  energy  content  of  the  waves 
(Battjes,  1972b). 

Wind  force  is  constant  over  the  surf  zone  and  is  determined  from  equation  (11): 


s-  ’  j< 


xt  3S , 


y* 


dx 
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T;  x„  =  P  ty  dx  =  Pa  cj  W\  W  sin  Qxb  (43) 

C.  RESULTS  AND  DISCUSSION 

Resultant  spectra,  wind  force,  and  wave  force  for  both  a  gentle  beach  slope  and  a 
steep  beach  slope  are  illustrated  as  cases  A  through  F  in  Figure  5  and  Figure  6.  Hrms, 
h ,fm,  wind  speed,  Q,  and  xb  for  each  of  the  six  cases  are  listed  in  Table  2.  Increasing 
wind  speed  effect  for  a  gentle  beach  slope  (tan  /?  =  0.01)  is  depicted  in  cases  A,  B,  and 
C,  and  for  a  steep  beach  slope  (tan  /?  =  0.10)  in  cases  D,  E,  and  F.  Wind  and  wave 
direction  are  relative  to  normal  beach  incidence.  The  three  spectra  previously  described 
(Pierson-Moskowitz,  JON’S  WAP,  and  TMA)  are  included  to  illustrate  the  effect  of 
shoaling  on  the  wind-wave  spectra,  as  well  as  the  expected  shift  to  lower  peak  spectral 
frequency  as  wind  speed  increases.  The  Hm%  and  depth  (h)  at  commencement  of 
breaking  also  increase  with  increasing  wind  speed. 

Table  2.  CASE  PARAMETERS  FOR  WIND  VERSUS  WAVE  FORCE  COM¬ 


PARISONS  (see  Figure  5  and  Figure  6). 


Case 

tan  /? 

wind 

speed 

(ms->) 

Hrms 

(m) 

h(m) 

xb(m) 

fm(  Hz) 

Q 

A 

0.01 

10 

0.7 

1.6 

160 

0.13 

0.0021 

B 

0.01 

20 

3.0 

7.0 

700 

0.06 

0.0033 

C 

0.01 

30 

7.6 

18.0 

1800 

0.04 

0.0053 

D 

0.10 

10 

0.7 

1.6 

16 

0.13 

0.0021 

E 

0.10 

20 

3.0 

7.0 

70 

0.06 

0.0033 

F 

0.10 

30 

7.6 

18.0 

180 

0.04 

0.0053 

The  wind  force  is  directly  proportional  to  the  width  of  the  surf  zone  (xb),  which  is  a 
function  of  tan  /?.  As  a  beach  steepens,  surf  zone  width  narrows  and  thus  there  is  a  re¬ 
duced  area  for  surface  wind  stress  to  act.  Wave  force  is  independent  of  tan  /?. 

Wind  direction  dependency  is  inherent  in  equations  (41)  and  (43)  where: 

Syx  oc  sin  8  cos  8  =  sin(20)  (44) 
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Figure  5.  Comparison  of  Mind  and  wave  force  for  a  gentle  beach  slope  (tan  /J  = 
0.01)  with  associated  wave  spectra. 
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Figure  6.  Comparison  of  wind  and  nave  force  for  a  steep  beach  slope  (tan  fi  =  0.10) 
with  associated  wave  spectra. 
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(45) 


Ty  oc  sin  6 


Thus  the  maximum  force  occurs  at  ±  90°  (alongshore)  for  wind  and  ±  45°  for  waves. 
Wind  force  increases  more  rapidly  than  wave  force  for  increasing  6.  When  winds  are 
blowing  offshore  (i.e.  absolute  wind  directions  >  90°  ),  waves  do  not  form  because  of 
fetch  limitations  and  thus  only  wind  force  exists,  creating  a  weak,  but  totally  wind-driven 
longshore  current. 

The  dependency  of  wind  and  wave  force  on  wind  speed  is  not  as  readily  apparent 
because  of  the  interrelationship  of  wind  speed  with  the  various  TMA  parameterizations 
(equations  33  to  37).  The  absolute  force  of  both  wind  and  waves  increases  as  wind  speed 
increases.  The  increase  in  wave  force  is  due  to  the  increase  in  spectral  energy.  The  in¬ 
crease  in  wind  force  is  through  the  wind  speed  and  the  increase  in  surf  zone  width  caused 
by  larger  waves  (generated  from  the  TMA  spectrum)  breaking  farther  offshore.  Thus 
the  wind  force's  dependency  on  wind  speed  is  indirectly  dependent  on  the  spectral  wave 
energy. 

Wave  and  wind  force  are  compared  for  increasing  wind  speed  and  constant  wave 
incident  angle  (Figure  7).  Since  wave  force  is  independent  of  ,  only  wind  force  is  il¬ 
lustrated  for  both  gentle  (solid  line)  and  steep  (dashed)  beach  slopes.  The  amount  of 
increase  for  both  wind  and  wave  force  over  the  wind  speed  range  of  5  ms-1  to  30  ms-1 
is  nearly  four  orders  of  magnitude.  Wind  force  is  seen  to  increase  more  rapidly  than 
wave  force  for  increasing  wind  speed  and  constant  incident  wave  angle.  Varying  the  in¬ 
cident  wave  angle  results  in  the  same  conclusion;  therefore  only  6  equal  to  20°  is  illus¬ 
trated. 

The  more  rapid  rate  of  increase  for  wind  force  is  quantified  in  terms  of  a  wind  force 
to  wave  force  ratio: 


Syx 


(46) 


The  ratio  of  forces  is  depicted  in  Table  3  (tan  /?  =  0.01)  and  Table  4  (tan  /?  =  0.10)  for 
the  quadrant  of  wind  direction  from  0  -  90°.  The  ratio  is  symmetric  for  wind  directions 
of  0  to  —90°  and  is  a  function  of  slope,  wind  speed,  and  wind  direction.  For  steep  beach 
slopes,  wind  force  is  much  less  than  wave  force  and  thus  may  be  neglected.  However, 
for  gentle  slopes,  wind  force  increases  as  the  wind  becomes  more  alongshore  in  direction. 
For  a  given  wind  direction,  an  increase  of  10  ms-1  in  wind  speed  results  in  a  50%  in- 
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Figure  7.  Comparison  of  wind  and  wave  force  for  increasing  wind  speed  at  constant 
wind  direction  ( 0  =  20°). 
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crease  in  the  ratio.  Although  this  ratio  for  gentle  beach  slope  is  generally  less  than  unity 
for  all  wind  speed  and  direction  combinations,  wind  force  is  never  less  than  an  order  of 
magnitude  compared  to  wave  force.  Wind  to  wave  force  ratio  versus  beach  slope  for  a 
typical  mid-latitude  cold  frontal  passage  is  depicted  in  Figure  8.  Wind  and  wave  frontal 
conditions  are  defined  as  a  wind  speed  of  15  ms-1  and  6  =  20°.  The  wind  to  wave  force 
ratio  increases  with  decreasing  beach  slope.  The  ratio  is  approximately  10%  for  a  beach 
slope  of  0.04  and  approaches  40-50%  for  beach  slopes  less  than  0.01.  Therefore  wind 
force  should  be  considered  in  studying  longshore  currents  during  strong  winds  on  gentle 
beaches  and  during  offshore  winds.  In  Chapter  VII,  wind  to  wave  force  ratios  from  the 
field  experiment  will  be  compared  to  the  results  in  Table  3  and  Table  4. 


Table  3.  RATIO  OF  WIND  FORCE  TO  WAVE  FORCE  FOR  A  GENTLE 
BEACH  SLOPE  (tan  /?  =  0.01) 


Wind 

Speed 

(ms-1) 

Wind  Direction  (relative  to  normal  incidence) 

0° 

10° 

O 

O 

ri 

30° 

40° 

50° 

60° 

70° 

80° 

O 

O 

Os 

10 

0.00 

0.15 

0.15 

0.17 

0.19 

0.22 

0.29 

0.83 

INF 

20 

0.00 

0.23 

0.24 

0.26 

0.29 

0.35 

0.44 

0.65 

1.28 

INF 

30 

0.00 

0.34 

0.35 

0.38 

0.43 

0.52 

0.66 

0.97 

1.91 

INF 

•  INF  -  ratio  is  infinite  since  there  is  no  wave  force 

•  waves  are  assumed  to  be  locally  generated  and  to  approach  the  shoreline  from 
the  wind  direction 
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Table  4.  RATIO  OF  WIND  FORCE  TO  WAVE  FORCE  FOR  A  STEEP  BEACH 
SLOPE  (tan  j?  =  0.10) 


Wind 

Speed 

(ms-1) 

Wind  Direction  (relative  to  normal  incidence) 

0° 

10° 

20° 

o 

O 

40° 

50° 

60° 

o 

o 

00 

o 

o 

90° 

10 

0.00 

0.01 

0.02 

0.02 

0.02 

0.02 

0.03 

0.04 

0.08 

INF 

20 

0.00 

0.02 

0.02 

0.03 

0.03 

0.03 

0.04 

0.07 

0.13 

INF 

30 

0.00 

0.03 

0.04 

0.04 

0.04 

0.05 

0.07 

0.10 

0.19 

INF 

•  INF  -  ratio  is  infinite  since  there  is  no  wave  force 

•  waves  are  assumed  to  be  locally  generated  and  to  approach  the  shoreline  from 
the  wind  direction 
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Wind  to  Wave  Force  Ratio 


Figure 


Beach  Slope 


Wind  to  wave  force  ratio  versus  beach  slope  for  a  tvpical  mid-latitude  cold 
frontal  passage.  The  wind  and  wave  incident  angle  (6)  is  relative  to 
beach  normal. 


I 
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IV.  FIELD  EXPERIMENT 


A.  INTRODUCTION 

To  develop  an  improved  understanding  of  coastal  processes  (currents,  waves, 
sediment  transport,  and  nearshore  geomorphologv)  under  a  wide  variety  of  conditions, 
a  one  month-long  field  experiment  was  held  in  the  autumn  of  1986  at  the  L'.S.  Army 
Corps  of  Engineers,  Waterways  Experiment  Station's  Coastal  Engineering  Research 
Center  (CERC),  Field  Research  Facility  (FRF)  at  Duck,  N.C.  Duck  is  on  the  Atlantic 
Ocean  side  of  the  North  Carolina  Outer  Banks  and  is  located  approximately  10  miles 
north  of  Kitty  Hawk,  N.C.  This  major  field  experiment  was  entitled  SLPERDUCK  and 
included  nearly  20  scientific  investigations. 

B.  SITE  DESCRIPTION 

The  FRF  site  (Figure  9)  is  located  along  a  100-km  unbroken  stretch  of  shoreline 
of  an  extensive  barrier  island  formation  known  as  North  Carolina's  ''Outer  Banks.''  The 
barrier  island  upon  which  the  FRF  site  is  located  is  bordered  by  the  Atlantic  Ocean  to 
the  east  and  Currituck  Sound  to  the  west.  Except  for  a  few  fishing  piers,  there  are  no 
major  coastal  structures  or  littoral  barriers  along  the  entire  reach  of  shoreline.  The  site 
has  a  tidal  range  of  0.5  to  2.0  m  and  regular  offshore  bathymetry  free  of  features  which 
may  affect  wave  climate.  The  beach  is  characterized  by  a  mean  foreshore  slope  of  1:10, 
a  single  longshore  linear  bar  system  which  can  become  three-dimensional,  and  a  mean 
slope  of  1:100  offshore  of  the  bar.  Shoreline  excursion  due  to  tidal  fluctuations  is  typi¬ 
cally  10  m.  Mean  sediment  grain  size  is  approximately  0.5  mm. 

FRF  personnel  support  includes  a  staff  of  nine  (scientists,  technicians,  divers,  and 
heavy  equipment  and  computer  operators).  Extensive  equipment  support  includes  a 
LARC-V  amphibious  vehicle,  a  Coastal  Research  Amphibious  Buggy  (CRAB),  auto¬ 
mated  data  processing  (ADP)  equipment,  a  heavy-duty  forklift,  four-wheel  drive  and 
all-terrain  vehicles  (ATV's),  laboratories,  electronic  workshops,  a  maintenance  garage, 
and  high-resolution  Zeiss  Elta-2  electronic  surveying  systems. 

Several  planning  workshops  were  held  months  prior  to  the  experiment  to  synergize 
the  efforts  of  all  the  principal  investigators.  The  FRF  staff  and  Sl’PERDLCK  partic¬ 
ipants  jointly  installed  more  than  70  electronic  sensors  in  the  nearshore  environment  for 
the  experiment.  The  overall  set-up  of  the  various  instruments  at  the  FRF  site  is  illus¬ 
trated  in  Figure  10.  The  longshore  current  meter  array  was  approximately  500  m  in 
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length  and  25  m  offshore.  Its  primary  purpose  was  for  the  investigation  of  infragravity 
waves,  however  its  data  will  also  be  used  in  this  research  to  investigate  temporal  and 
spatial  variability  of  longshore  momentum  flux  and  to  provide  an  extra  data  station  for 
sled  transects  of  the  surf  zone.  Unfortunately,  the  cross-shore  array  of  current  meters 
and  pressure  sensors  was  not  functional  for  the  experiment.  The  third  cross-shore  sen¬ 
sors  were  sonar  bottom  profilers.  A  high-resolution  linear  pressure  sensor  array  in  8  m 
depth  provided  offshore  directional  wave  data.  High  resolution  bathymetry  data  were 
collected  on  a  nearly  daily  basis  between  profile  lines  165  and  275  by  laser  positioning 
the  CRAB  with  the  Zeiss  surveying  system  as  it  methodically  traversed  the  area  along 
the  profile  lines.  Example  products  of  the  surveying  effort  are  included  in  Figure  1 1  and 
Figure  12.  A  major  difference  between  this  and  previous  large-scale  nearshore  exper¬ 
iments  was  the  extensive  bathymetric  data  acquired  using  the  CRAB. 

C.  EQUIPMENT  DESIGN  AND  OPERATION 

Only  a  brief  description  of  the  sled  and  its  instrumentation  is  provided  here. 
Martens  and  Thornton  (1987)  describe  in  detail  the  design,  instrument  calibration,  and 
specification  of  the  sled  system.  The  all-aluminum  sled  (Figure  13)  was  constructed  by 
the  CERC  Field  Research  Facility  staff  in  Duck,  N.C.,  based  on  a  Naval  Postgraduate 
School  design.  The  sled  foundation  is  two  runners  3  m  in  length  on  which  was  placed 
a  3  m  x  3  m  aluminum  sled  frame  with  a  9  m  aluminum  sailboat  mast.  The  frame  was 
approximately  one  quarter  of  a  meter  above  the  sled  runners.  All  the  instrumentation 
was  then  attached  to  this  frame.  Significant  additional  weight  was  placed  on  the  runners 
to  provide  stability  and  prevent  dragging  in  the  robust  surf  zone.  The  sled  was  designed 
with  a  low  profile  to  minimize  flow  disturbance. 

The  differential  pressure  slope  array  concept  was  patterned  after  Bodge  and  Dean 
(1984)  with  some  major  modifications.  Two-inch  O.D.  stainless  steel  threaded  pipe  and 
pipe  fittings  were  assembled  in  the  shape  of  a  capital  '  H"  to  form  two  adjacent  triangular 
slope  arrays  with  a  common  vertex  at  the  center  (see  Figure  13).  The  five  differential 
pressure  ports  were  located  at  the  four  ends  of  the  "I"  and  one  at  the  center.  The  pipe 
cavity  was  filled  with  a  premium  grade  vacuum  pump  oil,  similar  to  mineral  oil,  which 
had  a  low  viscosity  coefficient  and  was  assumed  incompressible.  Four  Scnsym  model 
LX06002D  monolithic  differential  pressure  transducers  were  installed  at  the  ends  of  the 
"I",  just  inboard  of  the  diaphragms  which  were  exposed  to  the  sea.  These  differential 
pressure  transducers  provide  an  output  voltage  proportional  to  the  applied  differential 
pressure  with  a  response  time  of  0.1  seconds.  The  pressure  at  the  center  diaphragm  was 
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SUPERDUCK  BATHYMETRY  16  OCT  86 


Figure  11.  Three-dimensional  bathymetry  example,  16  October  1986.  Short  ver¬ 
tical  lines  represent  locations  of  cross-shore  and  longshore  instruments. 
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hydraulically  transmitted  to  the  positive  port  of  each  differential  pressure  sensor,  while 
the  sea  pressure  at  that  location  was  sensed  by  the  negative  port  of  the  sensor.  A 
Par o scientific  model  8020  digiquartz  pressure  sensor  was  installed  adjacent  to  the  center 
diaphragm  to  measure  the  absolute  pressure. 

Current  velocities  were  measured  using  three  Marsh-McBirney  bi-directional 
electromagnetic  current  meters  with  spherical  (4-cm  diameter)  probes  mounted  on  the 
sled.  These  current  meters  are  ideally  suited  for  surf  zone  studies  by  providing  accurate 
measurements  in  free  stream  turbulence  as  well  as  combined  steady;  oscillatory  currents 
(Guza  et  al.,  1986a).  Both  variable  vertical  stacking  and  slope  array  comparisons  were 
considered  in  positioning  the  current  meters  on  the  sled  frame.  Two  current  meters  were 
mounted  at  the  centroids  of  the  two  triangular-shaped  slope  arrays  at  elevations  0.7  and 

l. 0  m  above  the  bed,  with  the  third  current  meter  midway  between  at  an  elevation  of  1.5 

m.  During  sled  deployment,  the  current  meters  were  always  positioned  on  the  "up- 
current”  side  of  the  sled  so  as  to  avoid  flow  contamination  by  the  sled  mast. 

The  current  meters'  directional  electrodes  were  aligned  with  the  longitudinal  axis  of 
the  sled  frame  by  eye  and  are  considered  to  be  within  1  —  2°  of  the  true  axis.  During 
data  analysis,  a  computational  rotation  is  done  for  true  alignment  based  on  conservation 
of  momentum  flux  outside  the  surf  zone. 

A  single-point  R.M.  Young  model  05103  wind  monitor  was  mounted  atop  the  mast 
to  provide  "in-situ"  measurement  of  vector  wind  in  the  surf  zone  with  simultaneous  wave 
momentum  measurements. 

All  equipment  were  battery-powered  by  12  sealed  lead-acid  batteries  housed  in  a 
watertight  canister  mounted  on  the  sled.  The  analog  data  were  PCM  (pulse  code  mod¬ 
ulated)  encoded  and  then  telemetered  ashore  by  a  system  described  by  Lowe  et  al. 
(1972).  The  data  were  then  decoded  and  recorded  on  9-track  magnetic  tape.  Data  ac¬ 
quisition  times  for  all  the  various  SUPERDUCK  experiments  were  referenced  to  a  pre¬ 
cise  Eastern  Standard  Time  (EST). 

Visual  observations  of  waves  breaking  on  the  sled  mast  were  marked  by  an  electrical 
trigger  which  registered  a  pulse  on  the  data  tapes.  Although  extremely  tedious,  this  ef¬ 
fort  allows  determination  of  the  occurrence  of  wrave  dissipation  events. 

Precise  orientation  of  the  sled,  with  respect  to  true  north  and  the  SUPERDUCK 
experiment  area,  was  obtained  with  the  use  of  the  FRF's  Zeiss  Elta-2  electronic  survey¬ 
ing  system.  Two  reflective  prisms  (5-cm  diameter)  were  mounted  at  each  end  of  a  2.4 
m  longitudinal  mast  spreader  located  approximately  7  m  above  the  sled.  From  triangu- 


lation,  the  relative  positioning  of  these  two  prisms  and  their  absolute  positioning  within 
the  FRF  coordinate  system  provided  the  required  precise  orientation  (within  0.5°)  of  the 
current  meters  and  slope  array.  Relative  three-dimensional  positions  between  various 
points  on  the  sled  and  the  prisms  were  obtained  with  separate  Zeiss  shootings  prior  to 
and  after  the  sled's  deployment  in  the  SL'PERDLCK  experiment. 

Prior  to  selecting  an  area  of  the  surf  zone  to  transect,  the  previous  day's  bathymetric 
profiles  and  contours  were  studied  and  changes  in  morphology  noted.  The  geographic 
area  with  the  highest  degree  of  straight-and-parallel  isobaths  was  selected  for  each  day's 
operations.  The  runs  were  made  during  daylight  and  scheduled  for  1-2  hours  before  and 
after  a  tidal  change,  so  as  to  minimize  the  tidal  effect  and  variation  in  mean  shoreline 
location.  This  timing  also  coincided  with  the  recording  of  all  other  SL'PERDLCK  data, 
which  was  accomplished  within  ±  2  hours  of  a  tidal  change  on  a  24-hour  basis.  The  sled 
was  initially  towed  offshore  by  the  FRF's  CRAB  to  a  position  just  beyond  the  breaker 
zone  for  the  first  run  of  a  surf  zone  transect.  Subsequent  runs  were  made  at  the  point 
of  maximum  breaking,  on  top  the  nearshore  bar,  and  in  the  nearshore  trough.  Each  run 
was  at  least  35-minute  in  duration.  The  sled  was  brought  inshore  for  each  subsequent 
run  of  a  transect  by  towing  the  sled  with  a  chain  tether  attached  to  a  four-wheel  drive 
forklift  located  on  the  beach.  Zeiss  Elta-2  laser  surveying  shots  of  the  sled's  prisms  were 
taken  at  the  beginning  and  end  of  each  run,  with  many  runs  having  additional  surveying 
shots  during  the  run  itself. 

In  addition  to  the  sled's  anemometer,  simultaneous  meteorological  data  of  air  tem¬ 
perature,  sea  surface  temperature,  atmospheric  pressure,  and  wind  velocity  were  re¬ 
corded  by  other  SL’PERDLCK  investigators  at  the  end  of  the  pier  and  near  the  FRF 
laboratory. 
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V.  PRELIMINARY  DATA  ANALYSIS 

Data  analysis  of  the  SL'PERDUCK  field  experiment  is  addressed  in  the  next  four 
chapters.  Preliminary  data  analysis  is  addressed  in  this  chapter  and  includes  data 
screening  for  data  gaps  and  signal  interference,  instrument  orientation,  three- 
dimensional  depiction  of  the  mean  current  structure  over  the  bar,  and  a  discussion  of 
homogeneity  and  stationaritv  for  the  experiment.  Since  surface  wind  stress  is  a  major 
concern  and  few  wind  stress  studies  have  been  conducted  in  the  surf  zone,  Chapter  VI 
is  entirely  devoted  to  the  parameterization  of  surf  zone  surface  wind  stress.  The 
longshore  momentum  balance  is  addressed  in  Chapter  VIII.  The  relative  importance 
of  each  momentum  term  is  calculated.  Spatially-variable  c/s  are  determined  as  residuals 
from  the  momentum  balance.  A  longshore  current  model  for  planar  beaches  is  modified 
to  account  for  barred  topography  and  its  results  are  compared  with  observations  in 
Chapter  VIII. 

A.  INITIAL  DATA  SCREENING 

The  analog  data  were  PCM  decoded  and  then  digitized  at  8  Hz.  Time  series  plots 
for  all  16  channels  of  information  were  produced  for  the  entire  data  collection  period. 
These  records  were  visually  perused  for  data  gaps  and  signal  interference.  Basic  statistics 
(mean,  variance,  skewness  and  kurtosis)  were  calculated  and  formated  in  tabular  form. 
Temporal  plots  of  two-minute  means  and  variances  were  produced.  This  information 
was  used  to  eliminate  bad  data  collection  periods  and  data  channels.  Data  collection 
achieved  a  9S.5%  successful  collection  rate  with  over  40  hours  of  data  collected  over  a 
10-day  period.  Nearly  all  of  the  eliminated  data  were  identified  at  the  beginning  or  end 
of  a  data  run  and  therefore  temporal  continuity  was  maintained. 

The  wind  direction  portion  of  the  R.M.  Young  wind  monitor  was  not  calibrated 
before  its  unexpected  destruction  by  wave  forcing  on  18  October.  Therefore  the  mean 
wind  direction  measured  at  the  sled  is  unavailable.  The  mean  wind  direction  measured 
at  the  end  of  the  FRF  pier  (SethuRaman  et  al.,  1987)  was  used  for  this  research.  Mean 
wind  speed  after  18  October  was  also  taken  from  the  end-of-pier  data. 

Occasionally,  one  or  two  current  meters  w-ould  be  partially  or  totally  above  water 
during  a  sled  data  collection  period,  such  as  when  the  sled  was  atop  the  bar  at  low  tide. 
These  periods  were  identified  from  visual  inspection  of  the  time  series,  supplemented  by 
field  logbook  entries.  These  data  were  not  utilized  in  the  radiation  stress  calculations. 
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B.  SLED  INTER-ORIENTATION. 

A  portable  Zeiss  laser  ranging  system  was  used  to  precisely  measure  the  sled 
instrumentation  positions  relative  to  the  sled.  A  hand-held  reflective  prism  was  placed 
at  a  designated  measurement  point  on  the  sled  (e.g.,  the  right  rear  corner  of  the  sled 
frame  top)  and  that  position  designated  as  a  three-dimensional  reference  point  from 
which  all  other  sled  points  of  interest  (prisms,  current  meter  probes,  slope  array 
diaphragms,  etc.)  were  measured.  A  complete  laser  measurement  of  the  sled  was  con¬ 
ducted  on  26  September  1986  and  20  October  1986.  Sled  modifications,  such  as  instru¬ 
ment  adjustment,  mast  replacement,  etc.,  were  noted  in  the  field  logbook  so  as  to  enable 
selection  of  the  proper  survey  day  for  subsequent  calculations.  The  Zeiss  measurements 
were  supplemented  with  engineer  tape  measurements. 

The  sled  instruments  were  aligned  to  a  three-dimensional  local  cartesian  coordinate 
system  referenced  to  the  top  of  the  sled  frame  by  a  rectangular  coordinate  rotation: 

*iww-  *oid  cos  0  +  Told  sin  0 
Tnew  =  -''old  sin  0  +  Jold  COS  0 

where  0  is  the  angle  of  rotation.  A  table  of  absolute  three-dimensional  distances  be¬ 
tween  all  points  of  interest  was  then  generated.  These  distances  were  repeatedly  used 
throughout  the  analysis. 

During  field  operations,  three-dimensional  instrument  coordinates,  with  respect  to 
the  FRF  coordinate  system  were  obtained  from  the  sled  prism  three-dimensional  coor¬ 
dinates  by: 

•  relating  the  mast  spreader,  slope  array,  and  current  meter  system  longitudinal  and 
vertical  axes  to  the  sled  frame  longitudinal  and  vertical  axes, 

•  determining  instrument  pitch  and  yaw  angles  every  64  seconds  from  linear  in¬ 
terpolation  of  the  two-to-five  survey  fixes  for  a  sled  run  collection  period.  All  data 
were  then  adjusted  every  64  seconds  for  the  current  sled  pitch  and  yaw  angles. 

•  adjusting  current  meter  and  anemometer  heights  above  the  bottom  even  64  sec¬ 
onds  for  sled  pitch  and  yaw  angles,  and 

•  correcting  the  distance  between  current  meters  for  the  sled's  pitch  and  yaw  angles. 

Since  all  data  were  taken  from  a  mobile  sled  rather  than  from  fixed-bed  instrumen¬ 
tation,  it  was  imperative  that  the  data  were  accurately  translated  into  a  coordinate  sys¬ 
tem  relative  to  the  local  bathymetric  contours.  Use  of  the  Zeiss  laser  ranging  system  to 
accurately  measure  the  sled  points  of  interest,  to  fix  the  sled  in  the  experiment  area,  and 
to  provide  comprehensive  bathymetric  data  coverage,  enabled  this  translation. 
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C.  BEACH  ANGLE  ORIENTATION. 

Location  coordinates  obtained  from  Zeiss-sighting  of  the  sled’s  prisms  were  con¬ 
verted  into  sled  pitch  and  yaw  angles  so  as  to  correct  the  sensor  data  with  respect  to  the 
horizontal  and  vertical  and  the  FRF  coordinate  system.  Pitch  and  yaw  angles  were 
calculated  even'  64  seconds  from  linear  interpolation  between  the  two  most  current  Zeiss 
sightings  to  account  for  any  sled  movement  during  the  data  run.  The  instrument  array 
orientation  was  then  numerically  rotated  relative  to  the  measured  bathymethy.  The  se¬ 
lection  of  64  seconds  coincided  with  the  ensemble  averaging  time  for  the  spectral  calcu¬ 
lations. 

Sled  data  were  acquired  from  11-18  October  and  20-21  October  1986.  Locations  of 
sled  data  runs  within  the  SUPERDL'CK  experiment  site  are  illustrated  in  Figure  14. 
Incomplete  survey  sightings  of  the  sled  prisms  due  to  salt  incrustation  on  the  prism  lens 
and  Zeiss  operator  error  prevented  determination  of  precise  sled  location  on  13.  14  and 
20  October.  Therefore  these  three  days  of  data  were  not  used. 

A  complete  bathymetric  survey  was  not  done  on  11,  17  and  21  October  due  to  high 
surf  conditions  or  equipment  availability.  For  sled  collection  days  when  no  coincident 
bathymetry  was  surveyed,  the  previous  or  next  day’s  bathymetry  was  appropriately  se¬ 
lected  based  on  minimizing  any  storm-induced  morphological  change  (see  Table  5). 

A  mean  bathmetric  contour  direction  was  determined  by  applying  a  least-squares 
regression  fit  to  the  survey  data  offshore  and  located  within  approximately  I0U  -  150  m 
of  the  nearshore  bar.  This  is  the  area  of  maximum  bathymetric  gradient  with  minimal 
slope  farther  offshore.  The  bathymetry  inside  the  bar  was  more  complex  than  offshore 
the  bar.  however  the  shoreward  distance  over  which  it  could  affect  wave  action  was  small 
0(30-50  m).  The  most  influential  bathymetry  was  assumed  to  be  that  located  on  and 
offshore  the  bar. 

The  longshore  length  scale  over  which  this  regression  fit  was  applied  was  weighted 
(2.3  -  1/3),  with  emphasis  on  the  portion  of  the  beach  which  favored  the  approaching 
waves  (e.g.,  for  waves  approaching  from  the  northeast,  the  northern  offshore  contours 
were  favored  by  a  factor  of  2).  The  total  distance  assessed  was  varied  from  100  to  250 
m  (or  to  maximum  extent  of  precision  bathymetry)  with  resultant  mean  orientation 
varying  less  than  0.5°  between  various  distances  considered.  The  mean  orientation  of 
the  four  "best-fit’’  straight  lines  (corresponding  to  0.5  m  contour  intervals)  was  chosen 
as  the  mean  bathymetric  contour.  The  beach  orientation  angle  was  the  angle  between 
this  mean  contour  and  the  coordinate  system’s  y-axis.  The  beach  orientation  angle  was 
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Figure  14.  Sled  data  run  locations.  Open  circles  are  sled  run  locations  and  clc 
circles  are  fixed-bed  current  meter  locations. 
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then  assigned  as  a  constant  for  that  transect  and  used  in  conjunction  with  instrument 
orientation  which  is  described  later. 
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Table  5.  BATHYMETRY  SELECTION  AND  BEACH  ORIENTATION  AN¬ 
GLES 


Sled  data  dav 

Bathvmetrv  dav  Profile  1 

Beach  orientation  anele  (°)2 

11 

12 

250 

0.5 

12 

12 

250 

0.5 

15 

15 

255 

4.0 

16 

16 

197 

-0.5 

17 

16 

195 

-0.5 

18 

18 

195 

1.5 

21 

22 

170 

0.8 

1  Profile  locations  are  shown  in  Figure  10. 

2  Beach  orientation  angle  is  measured  between  the  FRF  coordinate  system's 

y-axis  and  the  mean  bathymetric  contour  from  on  top  and 
offshore  of  the  bar.  The  angle  was  measured  over  a  longshore  length 
scale  vary  ing  between  100-250  m. 


A  combined  angle  (8C)  was  used  to  convert  from  the  yaw  angle  measured  by  the 
|  Zeiss  to  the  correct  rotation  angle  upon  which  to  coordinate  transform  the  u,v  velocities. 

This  combined  angle  was  the  sum  of  the  sled  yaw  angle,  the  beach  orientation  angle,  the 
angle  between  the  current  meter  longitudinal  axis  and  the  sled  centerline,  and  the  angle 
between  the  mast  spreader  longitudinal  axis  and  the  sled  centerline.  This  combined  an- 
|  gle  (6.)  was  used  in  equation  (47),  restated  here  with  velocities  substituted  for  coordinate 

locations,  to  perform  the  coordinate  transformation: 

Wnew  —  woid  8C  +  sin  dc  (4S) 

lnew  ^old  ^ht  4"  lold  ^0$  ® c 

I 

Fixed-bed  current  meter  velocities  were  also  converted  to  the  FRF  coordinate  sys¬ 
tem  using  equation  (48)  where  8C  is  field-measured  orientation  angles  with  respect  to  the 
FRF  coordinate  system. 

»  D.  CURRENT  METER  ORIENTATION  PRIOR  TO  A  TRANSECT. 

Since  the  radiation  stress  gradient  is  conserved  outside  the  surf  zone,  any  difference 
in  measured  S^X(J]  between  the  current  meters  on  the  sled  is  considered  to  be  due  to  er¬ 
rors  in  current  meter  alignment  with  respect  to  each  other  and  to  the  sled,  current  meter 

I 


49 


I 


calibration,  or  to  imprecise  orthogonal  orientation  of  the  dynamic  measurements  with 
respect  to  the  assumed  straight-and-parallel  bathymetric  contours.  The  contours  outside 
the  bar  were  nearly  straight-and-parallel  (see  Figure  12  and  Appendix  B).  Therefore, 
current  meter  alignment  or  gain  was  considered  in  error.  Orientation  of  the  current 
meters'  axes  with  the  sled's  axes  was  considered  to  be  within  1  —  2°. 

To  correct  for  the  radiation  stress  difference  between  two  current  meters  outside  the 
surf  zone,  one  current  meter  was  assumed  to  be  true  in  both  gain  and  alignment.  The 
other  current  meter  was  then  rotated  numerically  using  equation  (48),  or  adjusted  in 
gain,  until  the  radiation  stress  gradient  was  minimized. 

Two  current  meters  and  two  adjustment  criteria  (gain  and  rotation)  represent  four 
initial  parameters  for  adjustment: 

•  assume  the  front  current  meter  is  correct  and  adjust  the  rear  current  meter  in  gain 
(RGAIN)  or  rotation  (RTVV1ST),  or 

•  assume  the  rear  current  meter  is  correct  and  adjust  the  front  current  meter  in  gain 
(FGA1N)  or  rotation  (FTWIST). 

Because  the  sled  instruments  undergo  significant  stress  in  the  surf  zone,  all  sled  transects 
commenced  with  an  orientation  run  outside  the  surf  zone  to  allow  for  daily  reorien¬ 
tation.  A  numerical  orientation  prior  to  each  transect's  calculations  was  employed  dur¬ 
ing  data  analysis.  The  resulting  possible  orientation  adjustments  for  each  transect  and 
parameter  (RGAIN.RTWIST,  etc.)  are  shown  in  Table  6. 

The  adjustments  and  resulting  total  radiation  stress  (Sy[)  were  found  to  be  linearly 
related,  therefore  an  iterative  method  of  linear  interpolation,  or  "regula  falsi"  method 
(Gerald,  1978)  was  utilized  to  determine  the  orientation  required  to  minimize  Sv£  be¬ 
tween  the  current  meters.  Convergence  to  within  a  radiation  stress  gradient  tolerance 
of  10~3  Jm~2  was  generally  obtained  within  six  iterations.  Syx(f)  spectra  with  these  ad¬ 
justments  were  perused  to  ensure  that  the  strict  numerical  definition  of  a  conserved  ra¬ 
diation  stress  gradient  achieved  a  realistic  and  comparative  spectral  shape  for  the  current 
meters. 

The  selection  of  orientation  by  gain  or  rotation  does  have  constraints  when  dealing 
with  bi-directional  radiation  stress  spectra,  which  commonly  occurred  during  SUPER- 
DUCK.  Examples  of  typical  bi-directional  spectra  are  schematically  illustrated  in 
Figure  15.  The  spectra  indicate  swell  with  a  10-second  mean  period  from  the  southeast, 
while  wind  waves  with  a  5-second  period  and  significantly  more  energy  approached  from 
the  northeast.  Spectra  from  the  front  current  meter  are  assumed  correct  and  arc  de¬ 
picted  by  a  solid  line,  whereas  the  spectra  from  an  incorrectly  oriented  or  calibrated  rear 


current  meter  (and  requiring  gain  or  rotation  to  make  its  spectra  coincide  with  the  front 
current  meter  spectra)  are  depicted  by  a  dashed  line.  In  the  first  example,  gain  is  the 
obvious  correction  whereas  in  the  second  example,  rotation  is  the  obvious  correction. 

All  three  current  meters  were  calibrated  in  a  flow  tank  before  and  after  the  exper¬ 
iment  by  different  laboratories.  The  agreement  in  pre-  and  post-calibration  was  ±  3% 
in  gain  and  0.01  —  0.04  ms~*  in  offset.  This  agreement  is  considered  excellent  for  meters 
exposed  to  the  robust  surf  zone.  With  this  agreement,  coupled  with  the  subjective  ana¬ 
lyses  of  the  different  spectral  shapes  corresponding  to  different  orientation  combina¬ 
tions,  and  the  field  accuracy  of  current  meter  alignment  (±  2°),  it  was  decided  to  utilize 
current  meter  rotation  as  the  orientation  criteria.  Rear  current  meter  rotation 
(RTWIST)  was  arbitrarily  selected  over  front  current  meter  rotation  for  the  entire  data 
analysis,  except  for  Run  18-1  for  which  front  gain  was  used.  Resultant  RTWIST  angles 
are  within  the  values  of  field  accuracy  alignment.  Sensitivity  to  this  rotation  and  the 
effect  of  other  rotations  (FTWIST  only  or  both  FTWIST  and  RTWIST)  is  addressed  in 
the  error  analysis  discussion  in  Chapter  VII. 

After  completing  this  rotation  for  the  rear  current  meter,  the  side  current  meter  was 
rotated  (STWIST)  in  a  similar  manner  and  its  Sy[  values  served  as  a  first-order  com¬ 
parison  to  the  other  two  current  meters  values  for  Sy[.  These  orientation  rotations  were 
performed  after  the  6C  coordinate  transformation  were  fixed  for  the  calculations  inside 
the  surf  zone  for  that  particular  transect.  This  post-calibration  procedure  provided  a 
means  of  improving  the  accuracy  of  the  current  meter  calibration  and  alignment.  The 
procedure  is  essential  to  obtaining  realistic  physical  results  from  the  data  and  is  a  dis¬ 
tinguishing  feature  of  this  research. 

E.  MEAN  CURRENT  OBSERVATIONS 
1.  Overview 

Three-dimensional  portrayal  of  the  mean  current  structure  is  illustrated  in  Fig¬ 
ure  16  and  Appendix  A  (Figure  35  through  Figure  40).  All  current  vectors  represent 
34.1 -minute  means.  The  numbers  indicate  the  data  run  location.  The  short  dotted  lines 
above  each  run  location  indicate  the  mean  water  level  (MWL)  measured  by  the  sled's 
pressure  sensor  for  that  location  and  time.  The  different  arrowheads  delineate  the  height 
of  measurement  as  illustrated  in  the  three-dimensional  box.  MSL  is  indicated  by  the 
long  dashed  line.  Only  current  meters  which  were  completely  submerged  for  the  entire 
run  duration  arc  depicted.  Thus  some  run  locations  have  only  one  or  two  current  vec¬ 
tors. 
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BEFORE 


AFTER 


Figure  15.  Constraints  on  current  meter  rotation  or  gain  correction  for  bi¬ 
directional  spectra  outside  the  surf  zone. 
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Table  6.  INSTRUMENT  ORIENTATION  PARAMETERS.  Rear  and  front  cur¬ 
rent  meter  gains  and  rotations  ("TWIST")  required  to  conserve  the  radi¬ 
ation  stress  gradient  outside  the  surf  zone. 


Run 

Rear  Gain(%) 

Front  Gain  (%) 

Rear  Twist  (dee) 

Front  Twist  (dee) 

12-1 

-5.25 

5.54 

-0.81 

0.80 

15-1 

1.04 

-1.03 

-0.29 

0.27 

16-2 

4.57 

4.38 

-2.86 

2.34 

17-1 

-0.58 

0.58 

0.12 

-0.11 

18-1 

5.20 

-4.94 

-2.14 

1.89 

21-1 

1.40 

-1.44 

0.02 

-0.02 

21-6 

- 

- 

-0.41 

0.40 

•  indicates  no  iterative  convergence  or  acceptable  spectrum  shape. 

•  'Twist's',  or  current  meter  rotations,  are  in  units  of  degrees. 


The  sled  was  first  deployed  on  1 1  October  during  the  first  of  two  SUPERDUCK 
storms.  High  surf  (//mo  greater  than  3  m)  prevented  sled  deployment  beyond  the  surf 
zone,  therefore  there  was  no  instrument  orientation  run  for  that  day.  The  high  surf  also 
prevented  completion  of  a  minigrid  survey  by  the  CRAB.  Therefore  the  bathymetric 
profile,  as  well  as  the  instrument  orientation  parameter  was  obtained  from  the  12  Octo¬ 
ber  data.  The  unrealistic  Cj  values  later  calculated  for  1 1  October  highlight  the  necessity 
for  an  instrument  orientation  prior  to  a  transect.  The  1 1  October  data  set  was  therefore 
deleted  from  further  consideration. 

The  12  October  data  revealed  a  meandering  rip  current.  Two  transects  were 
completed  with  Runs  12-1  and  12-4  being  the  orientation  runs  outside  the  surf  zone. 
Run  12-4's  data  were  not  used  because  of  a  Zeiss  system  malfunction.  Therefore  Run 
12-1  is  the  only  data  run  outside  the  surf  zone  in  Figure  36.  The  strong  offshore  flow, 
0(0.5  ms-1),  at  all  three  levels  above  the  bed  was  substantiated  by  comparison  with 
current  vectors  from  the  longshore  array.  Although  this  rip  current  eliminates  homo¬ 
geneity  for  SyJ  calculations,  it  provides  a  unique  data  set  for  a  rip  current  study  separate 
from  this  research. 

October  15  is  the  first  day  with  an  established  bar-trough  system.  Run  15-1  is 
the  orientation  run  outside  the  surf  zone.  Strong  flow,  0(1.0  ms-1),  is  apparent  inside 
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Figure  16.  Mean  current  structure  across  the  bar  on  16  October  1986.  Numbers 
represent  data  runs,  short  dotted  lines  are  MWL's,  and  distance 
alongshore  is  arbitrary. 
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the  surf  zone  with  a  50%  reduction  in  speed  outside  the  surf  zone.  All  speeds  are  nearly 
uniform  in  depth  with  slight  speed  reductions  for  sensor  elevations  closer  to  the  bed. 
Multi-level  current  directions  in  the  trough  and  outside  the  surf  zone  are  co-linear. 
Angular  deviation  of  currents  with  height  (maximum  20°)  are  apparent  at  the  location 
of  maximum  breaking,  with  the  lower  sensor  elevations  exhibiting  a  bias  toward  offshore 
flow  near  the  bed. 

October  16's  data  have  the  strongest  velocities  of  all  the  sled  data  with  speeds 
exceeding  1.0  ms-1.  Run  16-2  was  the  orientation  run  outside  the  surf  zone.  Current 
directions  and  speeds  were  nearly  uniform  with  depth,  with  current  speeds  exhibiting  a 
slight  decrease  closer  to  the  bed.  Strongest  flows  were  on  the  offshore  side  of  the  bar. 

Significantly  weaker  currents  were  measured  on  17  October.  Run  17-1  was  the 
orientation  run  outside  the  surf  zone.  The  three  trough  runs  all  produced  negligible 
current.  A  combination  of  low  tide  and  low  wave  heights  [//rm3  0(O.2ms-1)]  on  this  day 
may  have  caused  most  waves  to  break  on  the  bar  and  thus  cause  the  negligible  current 
in  the  trough.  The  onshore  flows  at  Run  17-4  and  17-5  were  due  to  the  current  meters' 
location  in  the  upper  water  column  of  the  shallow  water. 

Strong  flows,  0(1.0  ms-1),  were  again  measured  on  18  October.  Runs  18-1  and 
18-6  were  orientation  runs  outside  the  surf  zone  for  two  separate  transects.  Run  18-6 
was  the  start  of  a  second  transect  which  was  abruptly  terminated  when  the  sled  mast 
snapped  in  the  high  surf  during  Run  18-7.  A  nearly  uniform  horizontal  velocity  gradient 
(dV'Idx)  is  observed  with  a  maximum  current  located  offshore  and  just  before  the  bar 
crest.  Vector  uniformity  is  observed, 

Two  transects  were  conducted  on  21  October.  Runs  21-1  and  21-6  were  orien¬ 
tation  runs  outside  the  surf  zone.  Unfortunately  current  velocities  during  both  transects 
were  too  weak  and  variable  in  direction  for  meaningful  analysis.  Hence  this  data  set  was 
eliminated  from  further  consideration,  leaving  only  four  (15-18  October)  of  the  original 
nine  transects  for  further  analysis. 

2.  Discussion 

Profiles  of  longshore  currents  across  the  barred  beach  generally  showed  signif¬ 
icant  spatial  variability.  Maximum  flow  is  observed  at  the  steepest  offshore  slope  of  the 
bar,  generally  just  below  the  bar  crest,  while  weaker  currents  are  found  outside  the 
breaker  zone  and  in  the  trough.  A  typical  distribution  of  longshore  currents  across  a 
barred  beach  is  illustrated  in  Figure  16.  Longshore  currents  are  generated  after  wave 
breaking  (Run  3  and  4)  and  are  extended  outside  the  surf  zone,  possibly  by  turbulent 
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momentum  exchange.  The  vertical  distribution  of  longshore  and  cross-shore  velocities 
are  compared  to  normalized  water  depth  (Figure  17).  Current  measurements  for  one 
run  are  denoted  by  a  continuous  line  connecting  only  those  current  meters  which  were 
entirely  submerged  for  the  entire  run.  Thus  some  shallow  locations  have  only  one  or  two 
current  means  illustrated.  The  plot  of  October  16  longshore  velocities  serves  as  a  legend 
and  has  numbers  beside  the  mean  current  symbols.  The  numbers  represent  the  chrono¬ 
logical  order  of  runs  in  a  transect,  with  the  most  offshore  run  being  labeled  as  number 
1. 

The  magnitude  of  the  cross-shore  flow  is  generally  much  less  than  the  longshore 
flow  and  is  of  0(0.2  ms-1).  The  longshore  flow  velocities  on  three  of  the  four  data  days 
are  of  0(1.0  ms-1).  The  cross-shore  velocities  exhibit  a  relative  onshore  flow  in  the  upper 
third  of  the  water  column  and  relative  offshore  flow  in  the  lower  third  of  the  water  col¬ 
umn.  The  lower  mean  currents  after  wave  breaking  tend  to  flow  more  offshore,  thus 
exhibiting  a  limited  "spiral"  distribution  with  height  and  indicating  the  presence  of  an 
undertow  flow  (an  offshore  flow  near  the  bed).  The  observations  are  in  agreement  with 
cross-shore  theory  (e.g.,  Svendsen  ct  al.,  1987).  The  longshore  velocities  indicate  nearly 
depth  uniform  flow  with  a  slight  velocity  increase  with  elevation.  This  depth  uniform 
flow  confirms  the  appropriateness  of  the  depth-integration  assumption  for  longshore 
flow.  Longshore  vertical  current  structure  did  not  exhibit  a  logarithmic  profile,  therefore 
the  profile  method  for  the  determination  of  bed  shear  stress  coefficients  is  not  applicable. 

3.  Stationarity  and  Homogeneity 

Stationarity  over  the  transect  was  investigated  by  examination  of  Sy[  calculated 
from  measurements  outside  the  surf  zone  by  a  fixed-bed  pressure  and  velocity  ('puv') 
sensor  designated  "south  tripod".  The  formulation  for  this  calculation  is  discussed  later 
in  equations  (78)  and  (79).  Sequential  34.1 -minute  spectra  and  resulting  Sy[  were  cal¬ 
culated  for  15-1S  October  (Figure  IS  and  Appendix  A  (Figure  41  through  Figure  43). 
The  short  lines  between  the  arrowheads  at  the  top  of  the  SyJ  graphs  represent  sled  run 
durations,  with  the  numbers  above  the  lines  identifying  the  sled  run  number  of  that  day. 
These  comparisons  are  considered  a  stria  indication  of  stationarity  because  of  the  sen¬ 
sitivity  of  SyT  to  changes  both  in  wave  height  and  wave  incidence  angle. 

The  'puv'  data  were  collected  by  the  FRF  for  four-hour  periods,  therefore  these 
data  were  not  completely  synoptic  with  sled  deployment  times.  The  'puv'  data  were 
digitized  at  2  Hz.  Spectral  analysis  was  accomplished  by  Bartlett  ensemble  averaging 
which  resulted  in  31  degrees  of  freedom.  Chi-squared  confidence  limits  for  95%  confi- 
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Figure  17.  Mean  cross-shore  U  and  longshore  velocity  V  profiles  versus  normalized 
water  depth  for  15-18  October.  The  16  Oct  V  profile  serves  as  a  legend 
for  the  data  symbols,  where  the  numbers  denote  chronological  order 
of  runs  in  a  transect  with  the  most  offshore  run  labeled  as  number  1. 


dence  yields  a  range  for  spectral  estimates  (Syx(J))  between  0.66  and  1.81  of  the  spectral 
estimate  calculated.  Ninety-five  percent  confidence  limits  for  the  integrated  or  total  ra¬ 
diation  stress  (SyJ,  see  equation  80)  were  based  on  the  Fisher-Z  transformation  (Miller 
and  Freund,  1985): 


where  the  correlation  coefficient  (r)  is: 


cov„ 


r  —  ■ 


(51) 


and  the  covariance  is  calculated  by  the  integrating  the  co-spectrum  of  u,v. 

The  October  15  runs  inside  the  surf  zone  (15-2  through  15-4)  appear  to  have  a 
reasonably  constant  (within  the  bounds  of  the  confidence  limits)  SyJ  value  with  no  more 
than  a  10%  variation  about  the  mean,  Both  October  16  and  17  data  sets  have  a  constant 
SyJ  with  minimal  variability.  October  18  data  also  appears  to  have  a  reasonably  con¬ 
stant  SyJ,  with  the  exception  of  Run  18-5  which  occurred  at  the  time  of  an  SyJ  increase. 
Fortunately,  this  run  is  within  the  trough,  where  flow  conditions  are  not  as  rapidly  af¬ 
fected  by  offshore  change  as  are  locations  on  the  offshore  side  of  the  bar.  Even  with  this 
SyJ  increase,  the  variation  is  no  more  than  20%  about  the  mean. 

Tidal  variation  during  an  individual  data  run  was  approximately  5%  of  the 
mean  water  depth.  Tidal  variation  over  an  entire  transect  ranged  from  14%  -  50%  of 
the  mean  transect  depth,  with  the  longest  transect  duration  having  the  largest  tidal  var¬ 
iation.  The  tidal  range  for  each  transect  is  illustrated  by  the  short  dotted  lines  in 
Figure  16  and  Appendix  A  (Figure  35  through  Figure  40). 

Another  test  for  stationarity  is  illustrated  in  Figure  19  where  Hmo  and  wave  in¬ 
cident  angle  (a)  calculated  from  the  'puv'  sensor  are  plotted  from  0500  -  2000  E.S.T.  for 
each  day  and  compared  with  sled  run  times.  Hmo  values  are  denoted  by  an  "x"  and  a 
values  are  denoted  by  an  "o".  Each  short  line  at  the  top  of  the  graphs  represents  the 
duration  of  a  sled  run,  with  the  number  above  the  line  denoting  the  sled  run  number  of 
that  day.  Comparison  of  Figure  19  with  Figure  18  and  Appendix  A  (Figure  41  through 
Figure  43)  offers  some  insight  into  the  nature  of  radiation  stress.  For  example,  the  de¬ 
crease  in  SyJ  observed  on  15  October  from  1000  -  1200  was  due  to  a  slight  reduction  in 
wave  height  and  a  significant  decrease  in  wave  incident  angle.  In  conclusion,  no  signif¬ 
icant  change  in  SyJ,  Hmo  ,  or  a  was  observed  during  the  times  of  sled  deployment. 
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Figure  18.  Offshore  S>x  spectra  and  SyJ  from  'puv'  gage  621  for  16  October 
1986.  Numbers  indicate  data  runs.  Intervals  arc  95%  confidence. 


Stationarity  requirements  are  therefore  assumed  satisfied.  It  should  be  noted  that  the 
stationarity  being  investigated  here  is  measured  over  relatively  long  temporal  scales 
(measured  in  hours)  whereas  the  temporal  term  examined  later  in  Chapter  VII  is  meas¬ 
ured  over  much  shorter  temporal  scales  (measured  in  minutes).  In  other  words,  the  as¬ 
sumption  of  stationarity  at  this  stage  of  the  data  analysis  does  not  preclude 
contributions  from  the  temporal  term  (e.g.  due  to  edge  waves)  to  the  momentum  bal¬ 
ance. 

Similarity  of  mean  current  vectors  in  the  experiment  area  was  considered  a  test 
for  homogeneity  inside  the  surf  zone.  Similarity  was  considered  satisfied  when  all  cur¬ 
rent  directions  were  in  the  same  quadrant  with  speeds  +  50%  of  an  overall  mean  for  the 
same  relative  on-offshore  position.  This  definition  was  chosen  considering  the  bounds 
of  instrument  accuracy,  current  vertical  variability,  the  number  of  instruments  in  the 
experiment  area,  and  the  variable  elevations  of  the  longshore  array  current  meters.  No 
visual  evidence  of  rip  currents  was  observed  during  15-18  October.  Two-dimensional 
plots  of  sled-determined  current  vectors  were  created  for  each  transect,  along  with  the 
fixed-bed  array  of  mean  current  vectors,  temporally  centered  for  the  sled  transect.  These 
comparisons  were  then  overlayed  on  the  coincident  bathymetry  (Appendix  B).  The 
bathymetry  allows  for  isolation  of  cases  where  bathymetric  forcing  is  dominant.  Vari¬ 
able  current  direction  is  observed  in  the  trough  on  17  October.  However  this  is  not 
considered  violation  of  homogeneity  because  of  the  extremely  weak  current  speed  in  the 
trough  on  that  day.  Perusal  of  this  data  indicates  the  assumption  of  homogeneity  for 
15-18  October  is  reasonable.  Data  sets  from  15-18  October  are  considered  to  meet  the 
dual  assumptions  of  homogeneity  and  stationarity. 
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OCT  16 


Figure  19.  Hm0  ('x')  and  wave  incidence  angle  ('o')  versus  time  for  15  -  18  October 


1986.  Numbers  indicate  data  runs. 
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VI.  SURFACE  WIND  STRESS  (zj) 

A.  INTRODUCTION 

The  objective  of  this  chapter  is  to  determine  representative  surface  wind  stress  for 
the  surf  zone  based  on  atmospheric  drag  coefficients  calculated  from  surface  wind  stress 
measured  outside  the  surf  zone  and  wind  speeds  measured  by  the  sled  anemometer  inside 
the  surf  zone.  Dr.  S.  SethuRaman  and  an  investigative  team  from  North  Carolina  State 
University  calculated  surface  wind  stress  for  the  end  of  the  FRF  pier  during  SUPER¬ 
DUCK.  Wind  speeds  were  measured  with  a  three-component  Gill  anemometer  at  a 
height  of  IS. 7  m  above  mean  sea  level  (MSL)  and  wind  stress  was  calculated  using  the 
eddy  correlation  method.  Although  this  wind  stress  measurement  was  outside  the  surf 
zone,  it  is  a  simultaneous  and  uncontaminated  measurement  over  water  and  within  500 
m  of  the  sled  transect  areas. 

A  synopsis  of  the  synoptic  scale  meteorological  conditions  during  October  1 9S6  is 
listed  in  Table  7.  Sled  deployments  weie  on  11-12,  15-1S.  and  21  October.  Predominent 
meteorological  features  were  onshore  winds  and  an  unstable  atmospheric  boundary  layer 
(as  indicated  by  the  air-sea  temperature  difference  in  Table  7). 

The  Gill  anemometer  has  three  polystyrene  propellers  which  are  rigidly  mounted 
with  their  axes  at  90°  from  each  other.  These  propellers  drive  three  d.c.  generators 
whose  output  is  directly  proportional  to  wind  speed.  The  anemometer  has  a  threshold 
of  0.1  ms-3  .  a  distance  constant  of  0.S  m.  and  can  be  sampled  at  100  Hz.  Although 
additional  wind  measurements  were  acquired  at  22  m  and  13.7  m.  the  eddy  correlation 
method  is  selected  over  the  profile  method  for  wind  stress  determination  due  to  its  in¬ 
creased  accuracy.  Air  temperature  (Ta)  and  relative  humidity  (RH)  measurements  were 
acquired  by  a  YSI  4401S  Thermilinear  thermistor  and  a  slow  response  relative  humidity 
sensor  (Model  207)  manufactured  by  Campbell  Scientific,  respectively.  Both  instruments 
were  mounted  at  IS. 7  m  above  MSU.  A  YSI  temperature  probe  measured  sea  temper¬ 
ature  ( Ts )  at  a  depth  of  0.6  m  below  MSL. 

Wind  stress  determined  by  the  eddy  correlation  method  is  computed  by: 


Table  7.  SUPERDUCK  SYNOPTIC-SCALE  METEROLOGICAL  DATA  (OC 
TOBER  1986)  (from  SethuRaman  ct  al..  19S7). 


Date  Wind  Field 

Brief  Description 

Ta  -  1\ 

Clouds  Remarks 

1 

SW  3-7.5 

High  dominates 

Clr 

2 

Yar.  3-7 

1  ligh  dominates 

5.5 

Cu.As  S-B  Hazv 

3 

SE-SW  1-6 

High  dominates 

1.7 

Ac. As  S-B 

4 

SW  7-10 

High  dominates 

7.9 

Clr 

5 

SW-NE  6-7 

Frontal  passage 

3.3 

Cu.Sc  S-B 

6 

NW-NE  eust  13 

Post-frontal  region 

-2.1 

Cu.Ac  S  Swell 

7 

NE-SE  3-5 

High  dominates 

-3.5 

Ac.Cu,  B-S 

S 

SE-SW 

High  dominates 

2.1 

Sc. As  B-0 

9 

NE-SE  1-5 

Cold  front  approach 

0.7 

Cu.Sc  S-B  Rain 

Hi 

NT  13  -  16' 

Cold  front  passage 

-2.3 

St. Sc  0  Rain 

li 

NE  11-16' 

Post-frontal  region 

-0.8 

St. Sc  O-B  Swell 

12 

NE  6-10 

High  dominates 

-0.5 

Ac. As  S 

13 

SE  3-8 

Cyclogenesis  along  shore 

1.0 

Cu.As  B  Rain 

14 

SW  6-S 

Cold  front  approach 

1.4 

TCu.St  0  Rain 

15 

NE  7-11 

Frontal  passage 

-6.8 

St.Cu  Swell 

10 

NE-NW  "  -  13' 

Post-frontal  region 

*» 

o .  ^ 

St. Sc  O-B  Swell 

I  7 

NE-NW  4-8 

High  dominates 

-2.3 

Cu.Ac  S 

IS 

NW-NE  9-14' 

High  dominates 

-a.  ' 

Cu  S-B  Swell 

19 

NE-NW  6-  13* 

High  dominates 

-2.9 

Cu  S-B 

20 

NE  5-S 

High  dominates 

-2.0 

Clr 

21 

NE-SE  1-3 

High  dominates 

-1.3 

Clr 

-i  -> 

SW  3-5 

High  dominates 

1.7 

Clr 

25 

SW  4-6 

High  dominates 

2.5 

Cs  S-B 

24 

NW-NE  3-9 

1  ligh  dominates 

O.0 

Ac.Cs  S-B 

25 

NE-E  7-10 

Warm  front  section 

0.0 

Sc. Ac  0  Ram 

26 

SE-SW  5-" 

Cold  front  approach 

4.1 

St.Cu  0 

27 

SW-NE  3-5 

(.'old  front  passage 

2.2 

Sc. As  B-0  Ram 

2S 

NE-SE  2-S 

1  ligh  dominates 

-0.6 

Cu 

29 

SW  1  -  5 

Transition  of  2  highs 

0.7 

Sc.C>;  B-O 

5 II 

SW-NE  4-7 

Front  sweeps  through 

-1.2 

Sc.Cs  B-O 

31 

NE  gust  12 

High  dominates 

-0.3 

Cu.Sc  S-B 

Comment'^: 

1. 

S:  scattered.  B:  broken.  Clr:  clear.  O:  overcast 

Cvclosenesis  occurred 

in  nearbv  recion  on  11th  and  16th 

3. 

Strong  upper-level  trough  accompanied  cold  front  during  14 

- 1 S  Oct 

4. 

*  indicates  high-energ> 

cases  of  SUPERDUCK  experiment 

5 

7'.  -  T,  was  measured  at  1400  ESI  which  was  the 

•A  J 

approximate  mid-time 

of  the  sled  transects. 

where  u.v,w  denote  wind  speed  in  the  horizontal  and  vertical  directions  and  the  primes 
(')  denote  deviations  from  the  mean.  Hourly  means  of  wind  stress  data  (SethuRaman 
et  ai.,  19S7)  allow  determination  of  an  atmospheric  shear  velocity  («.): 


15:1 


where  r  is  the  surface  wind  stress. 

Inherent  in  the  iu  determination  is  the  inclusion  of  atmospheric  stability.  The  de¬ 
scription  of  the  turbulent  regime  of  a  temperature-stratified  friction  layer  is  based  on  the 
Monin-Obukhov  similarity  theory,  which  defines  the  stability-dependent  vertical  wind 
speed  gradient  as: 


cu 


cz 


where  z 
L 

K 


=  height  above  the  surface  (m), 

=  Monin-Obukhov  length, 

=  von  Harman's  constant  (0.4),  and 
=  a  stability  function. 


(53) 


The  non-dimensional  ratio  z  L  is  a  measure  of  atmospheric  stability  with  negative 
values  denoting  unstable  conditions  and  positive  values  denoting  stable  conditions,  z  L 
is  defined  as; 


where  T,.  =  virtual  temperature  (degrees  Kelvin),  and 

T"v  =  virtual  temperature  scaling  parameter. 


(54; 


Tv  is  often  considered  to  be  the  mean  virtual  potential  temperature  within  the  surface 
layer  (Geernaert  et  al.,  1987a)  rather  than  its  more  strict  definition  of: 

Ty=Ta(  l  +  .61t/)  (55) 

where  q  =  specific  humidity.  In  this  research,  Ta  in  degrees  Kelvin  was  substituted  for 

I'.,. 


Integrating  equation  (53)  from  the  surface  (z  =  z0,  where  u  =  0)  to  an  arbitrary 
height  (z  =  z.  u  =  «,),  an  expression  for  a  stability-dependent  logarithmic  wind  velocity 
profile  is  obtained: 


u. 


l\  =  —  (ln-f-yj-r)) 


(56) 


Determination  of  z,  L  requires  an  iterative  process,  since  Tv.  is  also  a  function  of 
stability.  With  all  z  L's  initially  set  to  zero,  and  the  humidity  roughness  scale  (roa)  and 
temperature  roughness  scale  (zct)  set  equal  to  10— 6,  an  iterative  approach  (Appendix  C) 
based  on  the  stability  flux  profile  relation  of  Businger  et  al.  (1971)  is  utilized  to  determine 
z  L.  and  Yq  are  assumed  equal  and  formulated  as  (Large  and  Pond,  1981): 


Yh  =  Y.: 


x  =  (1  -  16)  -j-  )i-i 

Y  =  2  ln[  1  ] 

T  =  -7(f) 


T<0 

t^u 


p7» 


The  air  at  the  air-sea  boundary  is  assumed  saturated. 

Tire  stability  function.  $m.  is  described  by  the  Businger-Dycr  formula  for  unstable 
stratifications  (Businscr  et  al..  1971)  as: 


4>,n  =  ( 1 


is- 


(5Sl 


The  expression  lor  4',,,  is  (Paulson.  1970): 

1  4-  Q,  .  1  rj)  ,  .  i  r 

Ym  =  2  ln[ - ]  +  ln[ - ]  -2  tan~  V  +  - 


(59) 


which  is  an  empirical  representation  of  comprehensive  observations  over  smooth  terrain. 

B.  DETERMINATION  OF  AN  ATMOSPHERIC  DRAG  COEFFICIENT  FOR  THE 
SURF  ZONE 

Atmospheric  shear  velocity  and  L  are  independent  of  height  in  the  surface  boundary 
layer  (Panofsky  and  Dutton.  19S4)  and  are  assumed  to  apply  for  the  sled  transect  lo¬ 
cations  in  the  surf  zone  located  within  approximately  5<><i  m  of  the  surface  wind  stress 


measurements.  Atmospheric  shear  velocity  and  L  for  the  sled  runs  are  linearly  interpo¬ 
lated  from  the  hourly  pier  data  for  the  temporal  midpoints  of  each  sled  run. 

Boundary  layer  investigators  universally  convert  drag  coefficients  to  be  applicable 
at  10  m  height  to  allow  intercomparison  between  different  investigations.  For  equation 
(56)  to  be  valid,  a  stability-dependent  wind  speed  at  10  m  (F'10)  is  required.  i'l0  is  ob¬ 
tained  from  equation  (56)  as: 


L’10  is  first  calculated  from  the  hourly  pier  data.  Then  two  drag  coefficients  for  a 
height  of  10  m  at  the  pier  are  determined  —  the  stability-dependent  CJ10  (from  equation 
11  with  i]0  —  Uc  substituted  for  W)  and  the  stability-independent  or  neutral  drag  coef¬ 
ficient,  Q„10  (from  equations  11.  52,  and  56): 


Qio  ~ 


Pa  (^10  -  Li 


6' /in  i  n 


r  i.  f  '^f1  n_; 
L  I  *  JJ 


I  he  surface  boundary  condition  (  U0  ).  represented  by  a  mean  current,  was  taken  into 
account  to  determine  the  wind  speed  relative  to  the  moving  water  surface.  L'0  is  assumed 
zero  for  the  end  of  the  pier.  C0  is  represented  by  the  mean  surf  zone  current  for  the  sled 
locations. 

Wind,  wave  height,  atmospheric  stability,  and  initial  pier  Qdata  versus  time  foi  the 
seven  days  of  sled  deployment  are  shown  in  Figure  20  and  Appendix  A  (Figure  44 
through  Figure  46).  Drag  coefficients  and  pier  and  sled  wind  speeds  are  adjusted  to  10 
m  elevation.  Sled  anemometer  winds  are  represented  by  the  first,  mid.  and  last  two- 
minute  mean  wind  of  the  run.  The  clear  arrowhead  depicts  the  hourly  mean  wind  vector 
for  the  pier  wind  and  the  solid  arrowhead  depicts  a  34-minute  mean  wave  direction  sec¬ 
tor  for  a  location  near  the  pier  end.  The  point  of  the  arrowhead  indicates  the  time  for 
which  the  vector  is  valid.  Wind  and  wave  directions  are  relative  to  the  shoreline  which 
is  represented  by  the  y-a\is. 
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16  OCT  86 


Figure  20.  Wind  speed,  wave  height,  atmospheric  stability,  and  initial  Cd  data  for 
October  15-16,  1986,  Beach  is  assumed  parallel  to  the  y-axis.  Open 
arrowheads  denote  wind  direction  and  closed  arrowheads  denote  wave 
incident  angie.  Winds  are  adjusted  to  z  =  10  m. 
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The  absence  of  hourly  10  L  (z/L  stability  indicator  for  z  =  10  m)  data  and  wind  data 
indicates  no  boundary  layer  information  was  available.  Absence  of  hourly  10  L  data 
when  a  coincident  wind  speed  is  depicted  indicates  a  10  L  value  more  negative  than  the 
limit  plotted. 

Boundary  layer  stability  is  a  primary  concern  in  interpreting  the  wind  data.  There 
were  no  stable  atmospheric  conditions  during  SUPERDUCK.  Q„10  was  plotted  so  as  to 
remove  the  effect  of  stability  in  an  effort  to  relate  the  drag  coefficient  variation  to  an¬ 
other  physical  parameter,  such  as  wave  height,  wind  direction,  or  wind  speed.  Q„,0  is 
the  standard  drag  parameter  utilized  for  wind  stress  comparison  between  different  in¬ 
vestigators  and  serves  as  a  basis  for  examining  wave  influence  on  wind  drag.  Neither 
wave  height,  wave  direction,  nor  wind  direction  exert  obvious  influence  on  the  drag  co¬ 
efficient  as  shown  in  Figure  20  and  Appendix  A  (Figure  44  through  Figure  46). 

Increased  instability  places  the  atmosphere  in  a  convective  instability  or  local  free 
convection  regime  where  buoyancy  forces  become  equal  or  greater  in  magnitude  than 
mechanical  mixing  forces.  Wyngaard  and  Cote  (1971)  conducted  measurements  of  ver¬ 
tical  velocity  and  temperature  under  unstable  conditions  and  observed  local  free  con¬ 
vection  to  occur  at  z  L  <  -0.5.  Using  the  Businger  et  al.  (1971)  stability  flux  profile 
relation.  Garrett  (1977)  showed  that  as  instability  increases,  not  only  docs  Q  increase, 
but  roughness  becomes  significantly  more  important.  The  Businger  et  al.  (1971  >  iterative 
approach  is  inapplicable  in  a  local  free  convection  regime.  Therefore  an  empirical  ap¬ 
proach  is  deemed  more  appropriate  for  surf  zone  Q  determination  under  local  free 
convection. 

Sled  runs  of  z  L  <  -0.5  generally  occurred  during  periods  of  low  onshore  wind  speed. 
Therefore  the  drag  coefficients  in  Figure  20  and  Appendix  A  (Figure  44  through 
Figure  46)  depicted  for  low  wind  speeds  and  very  unstable  conditions  are  inaccurate. 
Onshore  winds  of  low  speed  and  diabatic  (unstable)  conditions  arc  representative  of  a 
sea  breeze.  A  sea  breeze  is  a  local  coastal  wind  caused  by  the  temperature  difference 
when  the  sea  is  colder  than  the  adjacent  land.  It  usually  occurs  on  relatively  calm  and 
sunny  days  (Iluschke,  1959).  Review  of  surface  weather  analysis  maps  (Appendix  D) 
and  wind  speed,  wind  direction,  and  stability  (Table  7  and  Table  9)  indicate  sea  breeze 
conditions  for  most  of  the  sled  data  run  periods. 

Measurements  of  surface  wind  stress  above  the  swash  zone  under  the  effect  of  a  sea 
breeze  were  made  by  simultaneous  wind  and  temperature  profiles  over  a  shore  transect 
near  Fort  Walton  Beach.  I  L  (Hsu,  1970).  The  sea  breeze  in  the  surface  boundary  laser 
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was  determined  to  be  in  tbe  atmospheric  free-convection  regime.  A  least-squares  fit  to 
the  data  ( figure  21)  yields: 

Qio  X  10'=  10-06-  1.30(17,0  -  r„)  (6?) 

where  L\0  is  the  mean  wind  speed  at  10  m  above  the  sea  surface.  It  is  noted  that  the 
drag  coefficient  under  the  sea  breeze  condition  is  higher  than  that  under  more  stable 
conditions  which  is  in  agreement  with  e.g..  Roll  (1965). 


Figure  21.  An  empirically-derived  stability-dependent  drag  coefficient  obtained  un¬ 
der  a  sea-breeze  condition  (from  Hsu.  1970). 

To  determine  appropriate  drag  coefficients  for  the  sled  transects  based  on  the 
SUPERDL'CK  wind  stress  data,  the  following  procedure  is  employed: 

•  When  the  sled's  anemometer  failed  (Runs  18-1.  1S-2.  and  all  runs  on  21  October), 
a  linearly-interpolated  hourly  mean  f'l0  from  the  pier  is  used  as  the  sled  i'l0. 

•  A  linearly-interpolated  Q,0  from  the  pier  is  assigned  for  the  temporal  mid-point 
of  the  sled  run. 


•  When  the  actual  z  L  was  more  negative  than  -1  or  the  wind  speed  was  less  than 
6.5  ms-1.  Cd !0  is  determined  from  the  empirical  relation  of  equation  (63)  w  here  the 


sled's  L'10  is  determined  from  equation  (60)  and  the  z  represents  the  run-varying 
height  of  the  sled  anemometer  above  mean  water  le\el  for  that  particular  run. 

•  When  z  L  was  more  negative  than  -1  and  the  wind  speed  greater  than  6.5  ms'1, 
QI0  is  linearly  interpolated  between  the  value  preceding  and  following  the  one  in 
question.  This  is  felt  to  be  a  better  approximation  than  not  using  stability. 

C.  ADJUSTMENT  OF  Q  FOR  THE  SURF  ZONE  ENVIRONMENT 

Since  the  surface  wind  stress  was  measured  outside  the  surf  zone,  a  Q  adjustment 
is  required  to  more  realistically  reflect  the  surf  zone  environment.  Early  investigators 
of  oceanic  Q  s  sought  a  dependence  on  wind  speed.  Subsequent  investigators  (e.g..  Hsu. 
1974,  1986a;  Byrne,  1982;  Gecrnaert  et  ah,  1986)  began  to  search  for  a  sea  state  de¬ 
pendence  for  Q,  since  surface  gravity  waves  are  assumed  to  act  as  roughness  elements 
( z0 )  for  the  momentum  transfer  process.  Most  of  these  investigations  w  ere  done  for  open 
ocean  conditions  and  therefore  are  not  directly  applicable  to  the  surf  zone. 

1.  Uniqueness  of  surf  zone  air-sea  interaction 

The  uniqueness  of  surf  zone  air-sea  interaction  as  opposed  to  open  ocean  air-sea 
interaction  can  be  summarized  by  three  factors  -  shoaling  of  the  bottom,  wave  refrac¬ 
tion.  and  wave  breaking. 

a.  Shoaling  effect 

Momentum  transfer  is  based  on  the  relative  wind  speed  difference  between 
the  absolute  wind  and  wave  speeds.  During  shoaling,  there  is  a  wave  steepening  and 
celerity  decrease  (relative  wind  increase)  of  surface  gravity  waves  as  they  shoal. 
Gecrnaert  et  al.  (1987b)  state  that  Q  increases  with  decreasing  values  of  the  ratio  c,,  u.. 
where  c0  is  the  phase  speed  of  the  dominant  locally-generated  gravity  wave.  Thus,  larger 
than  expected  drag  coelficients  are  observed  over  the  sea  during  shoaling  conditions 
(Smith.  19S0:  Gecrnaert  et  ah,  1986). 

b.  Refaction 

In  most  cases,  the  wave  propagation  direction  is  at  some  angle  (up  to 
180°)  to  the  predominant  wind  direction  due  to  refraction  of  the  waves  in  shoaling  water. 
Over  the  open  ocean.  Q  is  found  to  be  larger  if  the  wind  is  blowing  at  an  angle  different 
than  the  wave  direction  (Large  and  Pond,  19SI;  Smith.  1980).  Gecrnaert  et  al.  (1986) 
proposed  wind  stress  enhancement  in  advance  of  a  cold  front  due  to  surface  gravity  wave 
convergence,  Swell  generated  behind  the  cold  from  would  travel  through  the  front  and 
interact  with  waves  generated  ahead  of  the  front  and  propagating  parallel  to  the  front. 
This  interaction  would  result  in  a  chaotic  sea  of  larger  wave  amplitudes,  steeper  eras  ity 
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waves,  and  more  wave  breaking  --  very  similar  to  a  robust  surf  zone.  This  enhancement 
was  observed  by  Davidson  et  al.  (1988a,  1988b)  during  the  FASIXEX  1986  experiment. 
c.  Wave  Breaking 

Moderate  to  high  surf  conditions  undeniably  result  in  increased  surface 
roughness,  which  would  result  in  increased  drag  coefficients.  Although  some  work 
(Amorocho  and  DeVries,  1980;  Melville,  1977)  has  quantified  relations  between  breaking 
waves  and  increased  drag  coefficients,  these  investigations  have  dealt  with  the  open 
ocean  where  the  onset  of  breaking  waves  is  produced  by  a  dramatic  increase  in  wind 
speed. 

2.  Surf  Zone  Q  Adjustment  Formulation  and  Results 

The  wind  stress  data  from  the  end  of  the  FRF  pier  is  assumed  to  include  both 
the  shoaling  and  refraction  effects  on  Cd  as  one  transitions  from  the  open  ocean  to  the 
surf  zone.  The  increased  Ca  effect  due  to  shoaling  from  the  end  of  pier  to  the  surf  zone, 
as  well  as  the  additional  refraction  0(15  —  25°  maximum),  is  considered  to  be  small 
compared  to  the  breaking  wave  effect. 

The  approximate  mean  distance  from  the  swash  line  to  the  end  of  the  pier  is  500 
m.  Each  day  at  0700  EST,  the  distance  from  the  mean  swash  line  to  the  outer  extent 
of  the  surf  zone  was  "stepped  off'  along  the  pier.  These  rough  values  for  surf  zone  width 
on  sled  deployment  days  are  listed  in  Table  S.  The  surf  zone  nc .  .r  extended  to  the  end 
of  the  pier,  although  on  1 1  October,  the  farthest  surf  zone  extent  nearly  did.  With 
onshore  winds,  the  increased  roughness  would  not  take  effect  until  one  was  a 
measureable  distance  within  the  surf  zone.  Therefore,  the  increased  roughness  due  to 
the  breaking  waves  of  the  surf  zone  is  not  reflected  in  the  end  of  the  pier  wind  stress 
data. 

The  purpose  of  this  research  is  not  to  explicitly  characterize  the  Q  for  the  surf 
zone,  but  to  intelligently  select  an  appropriate  Q  with  which  to  address  longshore  sur¬ 
face  wind  stress  as  measured  by  a  single  anemometer  in  the  surf  zone.  The  increased 
surface  roughness  due  to  surf  zone  breaking  waves  may  be  conservatively  likened  to  the 
chaotic  seas  observed  by  Davidson  et  al.  (198Sa.  1988b)  in  advance  of  storm  fronts. 
From  a  composite  of  wind  stress  measurements  taken  from  multiple  chaotic  sea  cases 
found  in  advance  of  storm  fronts,  the  drag  coefficient  was  found  to  be  increased  by  ap¬ 
proximately  33%  over  that  expected  for  non-chaotic  seas  and  equal  wind  speed 
(Davidson  et  al.,  1988a.  19SSb).  Therefore  an  increase  of  33%  in  the  drag  coefficient 
(determined  from  equation  61)  was  deemed  appropriate  for  all  runs  except  those  runs 
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Table  8.  SUPERDUCK  SURF  ZONE  WIDTH  FOR  SLED  DEPLOYMENT 
DAYS  (from  Crowson  et  al.,  1 9SS). 


October 

Surf  zone  width  (m) 

11 

500 

12 

320 

15 

120 

16 

100 

17 

80 

IS 

80 

21 

100 

occurring  under  strong  diabatic  conditions  'where  equation  (6?)  was  utilized.  Hsu's 
(1970)  relationship  is  empirical  and  therefore  the  C.d  values  determined  by  equation  (63) 
would  already  have  the  increased  roughness  considered.  Drag  coefficient  values  for  each 
run  and  procedural  step  are  listed  in  Table  9.  All  the  values  are  of  O(10-3). 

Surface  wind  stress  is  determined  from  equation  (11)  using  the  C/s  listed  in 
Table  9.  Atmospheric  density  is  assumed  constant  at  1.25  kgm~3  and  6  is  obtained  from 
the  end-of-pier  wind  direction.  L'l0  —  l0  is  substituted  for  W  where  C10  is  the  stability- 
dependent  mean  windspeed  measured  by  the  sled  anemometer  for  the  entire  data  run 
(appro.Nimately  34  minutes  ),  corrected  to  an  elevation  of  10  nr.  U0  is  the  mean  surf  zone 
current  measured  by  the  sled  current  meters.  The  term  L'10  -  l  '0  is  obtained  from 
vectorial  subtraction. 


Table  9.  ATMOSPHERIC  DRAG  COEFFICIENTS  (C0l„  x  10')  FOR  EACH 
SLED  DATA  Rl'N 


Sled 

Run 

Sled 

10  ~ 

Pier 

*-  •10  “  *-  o 

C. 

lo 

L 

Tier 

Qio 

‘Sied 

Qio 

-;Sied 

Q/:o 

111 

(ms-1) 

14.44 

(ms-1) 

11.13 

(ms-1) 

0.56 

-0.025 

2.49 

3.31 

112 

13.69 

11.33 

0.57 

-0.023 

2.53 

3.36 

113 

13.2S 

11.25 

0.55 

-O.023 

2.41 

3.21 

121 

10.70 

7.2S 

0.30 

-0.068 

1.72 

122 

9.62 

S.15 

0.32 

-0.O4S 

1.54 

2.05 

123 

9.14 

7.86 

0.29 

-0.055 

1.41 

1.88 

125 

6.34 

4.79 

0.15 

-0.256 

1.07 

1.82 

126 

4.34 

0. 1 6 

-0.230 

1.48 

2.88 

127 

5. S3 

4.36 

0.09 

-0.640 

0.49 

2.48 

151 

9.49 

S.78 

0.30 

•Uj;  / 

1.1S 

1.57 

152 

s.oij 

7.76 

0.24 

-0.6<  10 

0.95 

1.2o 

153 

".S3 

7.14 

0.27 

-0.44S 

1.44 

1.92 

154 

7.57 

6.33 

0.14 

-1  410 

0.55 

1.92 

162 

12.51 

10.99 

0.37 

-o.l  2S 

1.15 

1.53 

163 

11.14 

11.30 

0.39 

-0. 1 1 6 

1.22 

1.62 

164 

10.42 

11.25 

0.39 

-0.123 

1.22 

1.62 

165 

9J6 

10. '2 

0.39 

-o.l  30 

1.33 

1.77 

166 

9.S9 

9.67 

0.36 

-0.158 

1.40 

1.86 

171 

6.25 

6.37 

0.20 

-o.4'3 

1.08 

1.94 

172 

".21 

4.77 

0.2' > 

-0.471 

1.31 

l.”4 

173 

.^0 

3.11 

0.21 

-( ).4~6 

1.35 

1.8o 

1~4 

6  >o 

5.27 

o.21 

-i  i.476 

1.39 

1.85 

1  "5 

5. 87 

5.55 

0.18 

-O.U64 

1.10 

2.43 

176 

5.63 

5.  IS 

0. 1 1 

-1.414 

0.49 

2.74 

1" 

5.S6 

4.S5 

o.OS 

-2.5o5 

o.3o 

2.44 

J~S 

5.S6 

4.54 

o.oS 

-2.604 

0.33 

2.44 

1S1 

lo.34 

10.61 

o.35 

-0.163 

1.11 

1.4S 

1S2 

10.39 

10.87 

0.40 

-0. 1 3o 

1.34 

1 .78 

1  S3 

11.74 

11.2o 

0.39 

-o.  140 

117 

1.56 

184 

12.17 

11.42 

0.4O 

-0.12" 

1.2o 

1.60 

185 

13.60 

1 1 .50 

0.42 

-o.l  1 1 

1.32 

1  "6 

1S6 

13.94 

1 1.80 

0.43 

-0.09  1 

1.32 

1.76 

211 

3.94 

3.96 

O.OS 

-2.147 

o.46 

4.94 

212 

4.49 

4.24 

0. 1 5 

-0.72S 

1.20 

4.22 

213 

3.66 

3.51 

0.14 

-0.690 

2.10 

5.30 

214 

2.26 

2.24 

0.13 

-0.775 

3.38 

7.12 

215 

2.36 

2.25 

o.oS 

-2.068 

1.25 

6.99 

216 

2.90 

2.89 

o.o3 

-14.286 

0.12 

6.29 

217 

1.94 

1.93 

0.04 

-6.645 

0.49 

".54 

2 1 S 

1.61 

1.57 

o.05 

-2.939 

1.38 

7.97 

219 

1.56 

1.39 

0.07 

-1.67S 

3.69 

8.05 

'eqn. (61) 

:eqn. 

(63)  3  33n 

o  increase  due  to  surf  zone  roughness 
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VII.  LONGSHORE  MOMENTUM  BALANCE 

The  longshore  momentum  balance  is  calculated  based  on  field  measurements  of 
wind,  waves,  and  currents  acquired  across  a  barred  beach.  Spatially-dependent  bed 
shear  stress  coefficients  and  the  relative  importance  of  the  various  terms  in  the  momen¬ 
tum  balance  are  calculated.  Wind  force  to  wave  force  ratios  are  calculated  from  the  field 
data  and  compared  to  the  results  of  the  wind  forcing  model  in  Chapter  III.  This  chapter 
concludes  with  an  error  analysis  of  these  calculations. 

The  slope  array  did  not  perform  as  required  for  determination  of  the  wave-induced 
longshore  momentum  flux  (S^)  in  equation  (1).  Pre-  and  post-calibration  values  for  the 
slope  array  differed  as  much  as  10%.  Assuming  that  breaking  wave  turbulence  outside 
the  surf  zone  is  negligible,  the  gain  required  to  match  the  slope  array's  radiation  stress 
(S„(/))  to  the  current  meter's  radiation  stress  (Syx(f))  was  0(30  -  60%).  Therefore,  slope 
array  data  were  judged  in  error  and  were  not  used  in  this  analysis.  Studies  addressing 
turbulent  momentum  flux  were  suspended. 

The  slope  array  loss  necessitates  a  rewriting  of  the  longshore  momentum  balance 
equation  as: 


cMv 


c: 


+ 


CSyi 


cx 


b  ,  rj 

-Ty  +  Ty 


(64) 


The  turbulent  momentum  flux  gradient  is  combined  with  the  wave-induced  momentum 
.flux  gradient  to  form  a  combined  momentum  flux  gradient  as  measured  by  the  current 
meters.  The  Sv[  term  implies  integration  over  frequency  (equation  80).  Determination 
of  surface  wind  stress  is  addressed  in  the  Chapter  VI.  Calculation  of  the  remaining  terms 
in  equation  (64)  will  be  addressed  separately  in  this  chapter. 

A.  TEMPORAL,  OR  ACCELERATION  TERM  ( cMyl?t ) 

The  time-averaged  depth-integrated  mean  momentum  per  unit  area  in  the  longshore 
direction  (My)  is  composed  of  a  steady  ( My )  and  unsteady  (My)  flow  component: 


\fy  =  My  +  Mv'  =  pV{tj  +  h)  +  J  p  ■/  dz 


(65J 
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The  steady  component  is  calculated  every  64  seconds  with  the  term,  >/  +  h  ,  being  taken 
as  the  mean  water  level  (MWL)  measured  by  the  pressure  sensor  for  that  period  of  time. 
This  resulted  in  thirty-two  mean  values  of  64-second  duration  for  My  . 

The  unsteady  flow  contribution  to  the  mean  momentum,  My,  is  attributed  to  the 
mass  flux  of  the  wave  motion: 


My'  =  p  v'  dz  (66) 

where  v’  is  only  known  at  z  =  0  from  wave  theory.  By  expanding  in  a  Taylor  series 
about  z  =  0  and  retaining  only  the  first  order  term,  My  may  be  expressed  independent 
of  height  in  the  crest-trough  region: 


m;  =  f  p  v'(0)  dz  =  f)  v!'(0) »/ 

Using  linear  theory  and  substituting  for  v'(0)  and  >/: 


My  =  -~pga2  1 


C 


A 

k 


E_  A 

C  k 


(67) 


(6S) 


where  a  is  wave  amplitude.  Rewriting  equation  (68)  in  terms  of  a  finite  frequency  band, 
a  significant  wave  incident  angle  (a(/))  ,  total  wave  energy  (E(f)).  and  wave  celerity  (C(f)) 
yields: 


My  =  ^  sin  *{f)  df  m) 

where 

C(J)  =  tanh  kh  (70) 

w(/) 

Throughout  this  analysis,  wave  number  (k)  is  calculated  non-iteratively  by  the  formu¬ 
lation  of  Wu  and  Thornton  (1986). 

Spectra  of  "significant  wave  angle"  (Appendix  F)  are  calculated  as  (Higgins  ct  al.. 


1981): 


•  'VA  SyX{f) 

sjn  «(/)  coso(/)  = 


(71) 


The  energy  spectrum  is  calculated  from  velocity  measurements.  Time-averaged 
depth-integrated  kinetic  energy  (KE)  may  be  written  as: 

KE  =  0.5  pj_^(w2  +  v2  +  w2)dz  (72) 

where  w  is  the  component  of  the  flow  velocity  in  the  vertical  direction.  To  the  first  order 
of  approximation,  there  is  equipartition  of  energy  between  potential  energy  and  kinetic 
energy  for  linear  progressive  waves  (Kinsman,  1965).  The  total  energy,  E(f)>  is  calculated 
as  twice  the  kinetic  energy: 

E(f)  -  2 1  //,(/)  | 2  [Gu(/)  +  Gv(/)]  (73) 

where 


I  *,(/)!’ - ,’CSmh,:t/,)3-  174) 

Ak{_  cosl i'k(h  4-  r,.,)] 

and  Gu  and  Gv  are  the  one-sided  spectral  density  functions  for  the  horizontal  velocity 
components  in  the  x  and  y  directions  and  zm  is  the  measurement  depth.  Calculations 
of  equations  (69)  through  (73)  indicated  that  the  time-varying  momentum  of  the  waves 
for  this  data  set  was  one-m  two  orders  of  magnitude  less  than  the  time-varying  mo¬ 
mentum  of  the  steady  flow.  1  nerefore  the  temporal  term  may  be  simplified  as: 

cMv  AA1V 

- —  ~  -  (75) 

ct  Ai 

where 

My  ^  My  =  p  V(>]  +  h)  (76) 

and  At  =  time  step  over  which  the  finite  difference  is  calci  lated.  and 
)/  +  /!  =  MWL  as  sensed  by  the  sled  s  pressure  sensor. 


A  value  of  My  is  first  determined  for  each  ensemble  of  64  seconds  in  the  data  run. 
Since  the  temporal  scale  of  the  longshore  current  is  unknown,  a  range  of  ensembles 


I 


-!6 
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(comprised  of  1  to  1 1  segments  of  64-second  duration)  are  defined  using  a  running  mean 
interval.  A/'.  Then  a  finite  difference  is  calculated  from  these  larger  ensembles  for  a 
specified  A t.  Finally,  a  mean  temporal  term  is  calculated: 


k 


7=1 


for  the  k  increments  of  AMyIAt.  This  mean  temporal  term  is  assigned  as  the  temporal 
term  for  the  sled  run.  Calculations  of  the  mean  temporal  term  (equation  77)  with  and 
without  a  running  mean,  as  well  as  variation  of  the  time  step  (with  subsequent  change 
in  the  k  increment),  made  negligible  difference  for  this  data.  Therefore  A i  =  64  seconds 
is  used.  The  selection  of  At'  is  discussed  in  paragraph  D  below. 

B.  RADIATION  STRESS  GRADIENT  (dSyx  T  /  dx) 

Longshore  momentum  flux,  as  defined  in  equation  (2),  is  calculated  from  current 
meter  measurements  and  applying  linear  wave  theory: 

SyXif)  =  \H2(J)\ 2  cuv{f)  (7S) 

where  |  H2(J)  1 2  is  a  linear  transfer  function: 

(79, 

4A'[  cosh lk{h  +  zm)3 

and  Cuv  is  the  real  part  of  the  cross-spectrum  of  u’  and  v\  It  is  assumed  in  the  applica¬ 
tion  of  the  transfer  function  (equation  79)  to  the  measured  covariance  that  the  vertical 
distribution  of  turbulence  in  the  water  column  is  the  same  as  in  the  waves,  and  is  de¬ 
scribed  by  the  transfer  function  which  is  a  function  of  depth  and  frequency 
(Figure  22). 

The  total  radiation  stress  (SyJ),  where  total  indicates  frequency  summation,  is  de¬ 
termined  from: 

ft 

Syx  =  X  VJO  A f  (80) 

h 
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Figure  22.  Linear  wave  theory  transfer  function  used  to  calculate  from  cur¬ 
rent  velocity,  showing  decreasing  frequency  dependence  with  depth 
across  a  barred  beach. 


where  f  and  f2  represent  the  lower  and  upper  frequency  limits  of  the  sea-swell  band  un¬ 
der  investigation  and  A f  represents  the  resolution  of  the  spectral  estimates.  S> J  is  the 
quantity  required  for  the  momentum  balance  of  equation  (64). 

Analysis  of  the  longshore  component  of  momentum  flux  ( Syx(f) )  was  accomplished 
using  Fast  Fourier  Transforms  and  ensemble  averaging.  Ensemble  averaging  requires 
segmentation  of  a  time  series  into  several  shorter  time  series  of  equal  length,  and  calcu¬ 
lation  of  a  sample  spectrum  for  each  of  these  sub-series.  The  mean  of  these  spectra  is 
a  smoothed  spectral  estimate.  The  advantages  of  ensemble  averaging  are  minimal  com¬ 
puter  storage  and  computational  ease. 

The  data  were  initially  demeaned  and  detrended.  The  34.1  minute  data  runs  were 
segmented  into  32  sub-series  of  length  64  seconds.  The  Syx(f)  used  in  equation  (80)  is 
the  Bartlett  smoothed  spectral  estimate  (Jenkins  and  Watts,  1968)  from  32  sample 
spectra.  All  spectral  estimates  (Syx,  rj,  etc.)  have  63  (v  =  2m  -  1)  degrees  of  freedom,  a 
spectral  resolution  (A/)  of  0.015  Hz,  and  Chi-squared  95%  confidence  limits  equal  to 
0.73  and  1.45  times  the  spectral  estimate  itself  (Jenkins  and  Watts,  196S).  Since  sea  and 
swell  are  of  interest,  lower  and  upper  frequency  limits  are  established  at  0.06  and  0.44 
Hz,  respectively.  These  limits  are  used  whenever  frequency  summation  was  required  for 
any  parameter. 

Radiation  stress  spectra  may  be  negative  at  one  frequency  while  positive  at  another 
frequency,  as  in  the  case  of  sea  and  swell  approaching  from  different  quadrants.  SVX{J) 
spectra  transitioning  from  positive,  through  zero,  to  negative,  and  vice  versa,  may  not 
be  accomodated  with  a  log-linear  plot,  therefore  Syx(f)  spectra  are  plotted  on  linear  plots. 
Related  spectra  are  also  plotted  in  this  manner  for  consistency. 

The  calculation  of  the  gradient  of  radiation  stress  requires  numerical  differencing 
over  an  appropriate  length  scale: 


dsj  ASj 


c.x 


Ax 


(SI) 


The  spatial  step  (Ajc)  was  initially  designated  as  the  distance  between  the  front  and  rear 
current  meters  after  adjustment  for  sled  yaw  angle.  A  radiation  stress  gradient  (equation 
81)  determined  with  this  definition  of  A.r  is  labeled  to  have  been  determined  by  the  "lo¬ 
cal"  method.  The  Ax  distance  on  the  sled  was  approximately  2  m  and  the  "local"  finite 
differencing  scheme 
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A  second  approach  to  the  determination  of  the  radiation  stress  gradient,  called  the  "in¬ 
tegrated"  method,  utilized  the  distance  between  sled  run  positions  for  Ax.  This  distance 
was  0(20-40  m)  and  this  finite  differencing  scheme  is: 
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(83) 


where  the  numbers  1  through  4  denote  the  two  pair  (front  and  rear)  of  current  meter 
data  available  from  two  sled  run  positions  (see  Figure  23).  The  estimate  of  the  S}J 
gradient  is  based  on  four  independent  current  meter  locations  covering  slightly  different 
transect  segments  and  slightly  different  Ax's.  These  segment  and  Ax  differences  are 
generally  less  than  10%  of  the  distance  between  sled  mast  locations.  It  is  therefore  as¬ 
sumed  that  the  four  estimates  represent  four  sample  means  of  the  true  Sy[  gradient  be¬ 
tween  the  two  sled  locations.  From  the  Central  Limit  Theorem,  the  standard  deviation 
of  sample  means  is  written  as: 


or  o asJ 
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a  g-\  . T 

dx 
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(84) 


where  n  is  the  number  of  sample  means.  Thus  the  standard  deviation  of  the  sample 
mean  is  reduced  by  50%  with  the  use  of  four  sample  means  rather  than  one.  This  fact 
will  be  used  in  determining  the  measurement  error  for  cSyx  T  /  cx  in  paragraph  E.  Se¬ 
lection  of  the  "local"  method  versus  the  "integrated"  method  proved  to  be  of  significant 
importance. 

C.  SEA  SURFACE  ELEVATION  (rj) 

Although  not  specifically  required  for  the  momentum  balance,  sea  surface  elevation 
spectra  provided  a  useful  reference  against  which  to  compare  radiation  stress  spectra. 

Sea  surface  elevation  spectra  ( Gv(f)  )  are  determined  from  linear  wave  theory  using 
data  from  both  the  sled's  pressure  sensor  and  current  meters: 
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Figure  23.  Finite  differencing  scheme  for  cS}X  T  /  cx.  from  integrated  method. 


g„(/)  =  \m\2  Gpjj) 


where 


\H3(f)\2  = 


cosn  (A/t) 
cosh 2[/c(/t  4-  zj] 


and  GpiJ)  is  the  one-sided  energy  density  spectrum  for  the  pressure  head  (ph)  measured 
by  the  pressure  sensor.  The  pressure  head  is  determined  from  the  hydrostatic  equation: 

P  ~  Patm  zo-n 


where  p  is  the  absolute  pressure  sensed  by  the  pressure  gage  and  palm  is  the  atmospheric 
pressure. 

The  formulation  for  determining  Gn{j)  from  current  meter  velocities  is: 


GJf)  =  I //400 1 2  lGu(J)  +  GM 


where 


sinh(2/c/fi 

Ha(j)  I  =  - rz - — 

2 gk  cosVlk(h  +  2m)l 


Grjif)  produced  by  these  two  formulations  are  in  good  agreement  for  the  three 
vertically-stacked  current  meters  and  the  pressure  sensor,  thereby  providing  confidence 
for  the  t]  spectra.  Sea  surface  spectra  produced  in  Appendix  F  are  formulated  from 
equation  (85). 

Wave  height  parameters  are  calculated  from  pressure-derived  t]  spectra: 
tfrms  »  (8  oh-T  (90) 


^mo  ^ 


where  the  total  variance  is  expressed  as: 


»  £  G„(f)  Af 
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D.  RESULTS  OF  THE  MOMENTUM  BALANCE 

Field  measurements  of  wind,  waves,  and  currents  are  used  to  calculate  the  local 
longshore  momentum  balance  at  various  locations  across  a  barred  beach.  From  the 
momentum  balance,  bed  shear  stress  coefficients  and  the  relative  importance  of  the 
various  terms  comprising  the  balance  are  calculated.  Mean  values  are  determined  rela¬ 
tive  to  the  bar.  Wind  force  to  wave  force  ratios  are  calculated  from  the  field  measure¬ 
ments  and  compared  to  the  results  of  the  wind  forcing  model  in  Chapter  Ill. 

1.  Bed  Shear  Stress  Coefficient 

A  spatially-dependent  bed  shear  stress  coefficient  (cy)  may  be  determined  from 
equations  (18)  and  (64)  as: 

8Svxt  „  d\fv 

c,  -  - °x-  ■  =L - £2-  (93) 

p\ut\v 

Resultant  spatially-dependent  bed  shear  stress  coefficients  are  listed  in  Table  10  and  il¬ 
lustrated  in  Figure  24.  The  column  titled  "c/l)"  are  the  c/s  determined  utilizing  the  "lo¬ 
cal"  method  for  the  radiation  stress  gradient.  These  coefficients  are  representative  for 
the  sled  run  position  itself.  The  surface  wind  stress  and  temporal  variability  terms  were 
34.1 -minute  means  taken  only  from  the  sled  run  being  investigated  in  anticipation  of 
gaining  insight  into  the  temporal  variability  of  all  terms  in  the  longshore  momentum 
balance.  A  running  mean  of  12.8  minutes  was  selected  for  determination  of  both  the 
temporal  term  and  the  radiation  stress  based  on  radiation  stress  statistics  of  Freilich  and 
Pawka  (1987).  They  showed  that,  given  n  independent  realizations  of  a  time  scries  of 
u(t)  and  v(t),  that  the  Gaussian  approximation  for  true  Syx  statistics  is  valid  for  n  >  11 
independent  realizations.  An  ensemble  length  of  64  seconds  and  a  minimum  number  of 
ensembles  (12)  was  selected  for  the  analysis.  Unfortunately,  this  procedure  resulted  in 
significant  c/l)  variability.  Utilization  of  a  34. 1 -minute  mean  6Syx  T  /  dx  made  negligible 
difference  in  resultant  c/l).  Two  possible  reasons  for  the  large  c/l)  variability  may  be 
that  the  radiation  stress  is  much  more  variable  across  the  surf  zone  on  a  small  spatial 
scale  than  previously  imagined  and  the  Ax  is  too  small  to  resolve  the  difference  between 
two  large  numbers  (see  equation  ($2)  and  paragraph  2  of  the  Error  Analysis  section  in 
this  chapter). 

The  column  titled  "c/i)"  are  those  c/s  determined  using  the  "integrated"  method 
for  the  radiation  stress  gradient.  These  coefficients  are  representative  for  the  portion  of 
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the  bathymetric  profile  between  sled  run  positions.  The  wind  stress  and  temporal  vari¬ 
ability  terms  are  the  34.1-minute  means  from  the  two  appropriate  sled  run  means.  The 
"beach"  position  along  the  transect  is  defined  as  the  intersection  of  the  bathymetric 
profile  with  the  MWL  as  determined  by  the  sled  pressure  sensor  on  the  run  closest  to 
the  beach.  Extremely  weak  longshore  currents  inshore  of  run  17-6's  location  explain  the 
variable  cAi)' s.  They  are  disregarded.  The  bed  shear  stress  coefficients  determined  by  the 
"integrated"  method  exhibit  less  variability  than  those  from  the  "local"  method  and  are 
all  of  G(10-3).  The  "integrated"  method  is  considered  to  be  a  reasonable  method  for  the 
determination  of  the  radiation  stress  gradient. 

2.  Relative  Contributions  of  the  Longshore  Momentum  Balance  Terms 

Relative  contributions  of  wave  forcing  (cSyJ  /  dx),  wind  forcing  (rfi,  and 
temporal  variability  of  mean  momentum  (SMy  j  ct )  to  the  total  momentum  balance  are 
determined  by  first  defining  their  total  contribution  (TOT)  as  the  cumulative  sum  of 
their  absolute  values: 


TOT 
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dMy 

6t 


(94) 


Their  individual  relative  contributions  (rc)  are  then  determined  by  dividing  by  the  total 
contribution: 


waverc  = 


cSyx 

cx 


TOT 


windrc  = 


TOT 


(95) 

(96) 


(97) 


Relative  contributions  from  the  wave,  wind,  and  temporal  terms,  as  formulated 
in  equations  (95)  to  (97)  and  utilizing  the  "integrated"  method,  are  calculated.  Results 
are  presented  in  order  of  consecutive  runs  (Table  1 1 )  and  by  relative  location  along  the 
transect  with  respect  to  the  bar  (Table  12).  Relative  locations  are  delineated  as  "off¬ 
shore"  of  the  bar,  "on  the  bar"  and  immediately  before  the  bar,  and  in  the  "trough"  re- 
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Table  10.  SUPERDUCK  BED  SHEAR  STRESS  COEFFICIENTS  (OCTOBER 
15-18,  1986) 


October  15.  1986 

October  17,  1986 

Run  Cy(/)  cj(i) 

Run  Cj{[) 

c/J) 

15-1  * 

17-1  0.002 

0.006 

0.003 

15-2  0.008 

17-2  0.003 

0.002 

0.005 

15-3  * 

17-3  0.004 

* 

0.003 

15-4  0.012 

17-4  0.026 

0.001 

0.005 

Beach 

17-5  0.099 

17-6  0.119 

0.003 

October  16.  1986 

October  18,  19S6 

Run  Cj(F)  Cj{i) 

Run  Cj{[) 

Cji‘1 

16-2  0.001 

18-1  0.001 

0.003 

0.003 

16-3  0.004 

18-2  * 

0.004 

0.004 

16-4  0.003 

1S-3  0.007 

0.001 

0.001 

16-5  0.0C1 

18-4  0.043 

0.001 

0.001 

16-6  0.002 

18-5  -.001 

0.002 

0.001 

Beach 

Beach 

•  indicates  value  <  5.  x  10-4 

•  cj{[)  are  determined  by  the  "local"  method 

•  c/i)  arc  determined  by  the  "integrated" 

method 

•  relative  contributions  determined  by  the  "integrated"  method 

gion.  The  mean  wind  stress  remained  nearly  constant  throughout  each  transect  of  the 
surf  zone.  However,  the  temporal  term,  as  averaged  over  34.1  minutes,  was  variable  in 
sign  and  magnitude.  The  location  of  maximum  wave  breaking  along  the  transect  can 
be  identified  (Table  11)  where  wave  forcing's  relative  contribution  is  maximum.  Two 
peaks  in  wave  forcing  correspond  to  two  breaking  areas  -  waves  breaking  offshore  of  the 
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bar  and  waves  reforming  and  then  breaking  on  the  nearshore  beach  face.  The  relative 
contribution  of  wind  stress  inside  the  surf  zone  generally  increases  in  value  as  one  pro¬ 
ceeds  onshore  from  the  location  of  maximum  breaking,  due  to  the  decrease  in  wave 
forcing.  The  relative  contribution  of  wind  stress  can  then  decrease  again  inshore  at  the 
location  where  waves  re-break  on  the  nearshore.  The  relative  contribution  of  the 
temporal  term  appears  to  be  largest  in  the  trough  where  wave  forcing  is  reduced  and 
possibly  where  the  effect  of  longer  period  waves  (e.g.,  surf  beat  or  edge  waves)  is  more 
pronounced. 

The  mean  relative  contributions  for  the  wave,  wind,  and  temporal  terms  (0.86, 
0.08,  and  0.06,  respectively)  are  essentially  identical  for  both  "offshore"  and  "on  the  bar" 
(Table  12).  However  in  the  "trough",  the  mean  relative  contributions  (0.64,  0.17,  and 
0.19,  respectively)  are  significantly  different.  The  relative  importance  of  wave  forcing  is 
reduced  by  25%,  whereas  the  relative  importance  of  surface  wind  stress  is  doubled,  and 
the  temporal  term  relative  importance  is  tripled. 

A  comparison  is  made  between  the  wind  force  to  wave  force  ratios  produced  by 
the  wind  forcing  model  in  Chapter  III  to  those  determined  for  the  SLPERDLCK  data. 
The  wind  to  wave  force  ratios  in  Table  3  and  Table  4  are  based  on  a  planar  beach, 
therefore  only  the  offshore  side  to  the  SLPERDL'CK  bar  (tan/?  *  0.03)  is  considered. 
The  SLPERDL'CK  force  ratios  are  obtained  from  ratios  of  the  relative  contributions 
of  surface  wind  stress  to  radiation  stress  gradient  in  Table  12  for  the  "offshore"  and  "on 
the  bar"  categories.  The  wind  force  to  wave  force  ratio  for  15-18  October  field  data  is 
0.09.  The  mean  wind  speed  was  9.6  ms-1  and  the  mean  wave  incident  angle  for  the 
"offshore"  region  was  12°.  Lsing  linear  interpolation  in  Table  3  and  Table  4  results  in 
a  wind  to  wave  force  ratio  of  0.1 2.  Thus  the  field  measurements  are  in  good  agreement 
with  the  results  of  the  wind  forcing  model  in  Chapter  III,  with  differences  attributable 
to  non-fully  arisen  seas  and  the  presence  of  swell  from  a  different  direction  than  the  wind 
waves. 

Bed  shear  stress  coefficients  are  calculated  for  cases  without  the  surface  wind 
stress  and  without  both  surface  wind  stress  and  the  temporal  term.  The  percentage 
change  in  cy  ranges  from  near  zero  to  100%.  with  a  mean  change  of  10-30%.  Thus  under 
certain  conditions  for  this  data  set,  the  surface  wind  stress  and  temporal  terms  can  be 
of  first  order  in  the  momentum  balance. 

A  mean  Cy  based  on  relative  postion  along  the  transect  with  respect  to  the  bar 
is  determined.  Bed  shear  stress  coefficients  determined  from  the  "integrated"  method  arc 
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Table  11.  RELATIVE  CONTRIBUTIONS  OF  THE  LONGSHORE  MOMEN¬ 
TUM  BALANCE  TERMS  FOR  SUPERDUCK 


Run 

cSyx  T 1  cx 

n 

dMy/Ct 

15-1 

0.97 

0.02 

0.01 

15-2 

0.93 

0.04 

0.04 

15-3 

0.77 

0.17 

0.06 

15-4 

0.79 

0.14 

0.07 

16-2 

0.S7 

o.os 

16-3 

0.95 

0.03 

16-4 

0.77 

0.19 

16-5 

0.70 

0.18 

16-6 

0.S1 

0.02 

0.63 

0.17 

0.2o 

0.94 

0.03 

0.04 

0.8S 

0.03 

0,09 

0.89 

0.05 

0.07 

0.72 

0.16 

0.12 

0.S3 

0. 1 3 

0.04 

0.93 

0.05 

0.02 

0.S1 

0.12 

0.07 

0.46 

0.33 

0.2) 

0.21 

0.19 

0.60 

•  Relative  contributions  determined  by  "integrated"  method.  Run  number  indi¬ 
cates  ofFshore  run  of  the  required  run  pair. 

•  Runs  17-6  through  17-8  had  negligible  currents  of  variable  direction.  Their  cj 
results  are  disrecarded. 


Table  12.  RELATIVE  IMPORTANCE  OF  MO¬ 
MENTUM  BALANCE  TERMS  FOR  LO¬ 
CATIONS  RELATIVE  TO  THE  BAR 


Offshore 

Run 

<55,.  r  T 1  ex 

n 

cAIJct 

15-1 

0.97 

0.02 

001 

16-3 

0.95 

0.03 

0.02 

16-2 

0.87 

0.08 

0.04 

17-1 

0.63 

0.17 

0.20 

17-2 

0.94 

0.03 

0.04 

18-1 

0.83 

0.13 

0.04 

18-2 

0.93 

0.05 

0.02 

Mean 

0.87 

0.07 

0.05 

On  the  Bar 

Run 

dS,r  T lex 

rl 

c'MJct 

15-2 

0.93 

0.04 

0.04 

16-4 

0.77 

0.19 

0.03 

17-4 

0.S9 

0.05 

0.07 

17-3 

0.88 

0.03 

0.09 

18-3 

0.81 

0.12 

0.07 

Mean 

0.86 

0.09 

o.06 

Trough 

Run 

<55,  r  T 1  ('x 

rl 

eMjct 

15-4 

0.79 

0.14 

0.07 

15-3 

0.7  7 

0.17 

0.06 

16-5 

0.70 

0.18 

0.12 

16-6 

0.S! 

0.02 

0.17 

17-5 

0.72 

0.16 

0.12 

18-4 

0.46 

0.33 

0.21 

18-5 

0.21 

0.19 

0.60 

Mean 

0.64 

0.17 

0.19 

assigned  and  a  mean  Cj  determined  for  each  of  these  relative  locations  (Table  13).  The 
result-  indicate  a  decrease  in  mean  Cf  as  one  proceeds  from  offshore  to  onshore  with 
mean  Cy  values  equal  to  0.004  ("offshore").  0.002  ("on  the  bar"),  and  0.001  (trough  ). 


The  shoreward  decreasing  cy  was  initially  thought  to  be  the  result  of  decreasing  bed 
roughness.  As  previously  stated,  bed  roughness  is  a  function  of  both  sand  characteristics 
and  bed  form.  Beaches  generally  exhibit  a  shoreward  gradation  of  sand  grain  size  from 
large  to  small.  Waves  tend  to  sort  the  sand  with  the  largest  sand  residing  in  the  most 
energetic  area  of  the  surf  zone.  This  gradation  would  be  commensurate  with  the  maxi¬ 
mum  cy  found  just  inside  the  surf  zone  and  a  decreasing  cy  further  onshore.  Unfortu¬ 
nately,  field  measurements  of  bed  roughness  were  overlooked  during  the  SUPERDUCK 
experiment  and  cannot  be  used  to  substantiate  this  hypothesis.  Interviews  with 
SUPERDUCK  divers  were  conducted  several  months  after  the  experiment,  however 
their  information  about  bed  roughness  was  not  specific  enough  to  draw  conclusions. 
An  October  1982  sampling  of  sand  grain  size  across  the  FRF  beach  showed  shoreward 
increasing  grain  size  (Birkemeier  et  al.,  1985).  It  is  reasonable  to  assume  that  grain  size 
gradation  would  not  change  over  a  few  years,  therefore  bed  roughness  is  discarded  as  the 
reason  for  the  Cy  trend.  Field  and  laboratory  studies  have  shown  that  bed  shear  stress  is 
enhanced  when  the  near-bottom  flow  field  is  composed  of  mean  current  and  oscillatory 
(i.e.  wave)  components  (Grant  and  Madsen,  1979,  1982;  Grant  et  al.,  19S4).  Thus 
shoreward  decreasing  wave  action  would  result  in  shoreward  decreasing  cy. 

Table  13.  MEAN  BED  SHEAR  COEFFICIENTS  DETERMINED  FOR  LO¬ 


CATIONS  RELATIVE  TO  THE  BAR 


Trough 

On  the  Bar 

Offshore 

15-4 

.001 

15-2 

.002 

15-1 

.006 

15-3 

<.001 

16-4 

.001 

16-3 

.004 

16-5 

.001 

17-4 

.005 

16-2 

.003 

16-6 

.002 

17-3 

.003 

17-1 

.003 

17-5 

.003 

1S-3 

.001 

17-2 

.005 

18-4 

.001 

18-1 

.003 

18-5 

.001 

18-2 

.004 

Mean  Cy= 

.001 

.002 

.004 

The  variation  in  cy  due  to  varying  At'  in  the  temporal  term  was  examined  by 
analyzing  sample  groups  where  At'  was  specified  as  a  running  mean  of  (approximately) 


1,  5,  10,  12,  14,  18,  and  22  minutes.  The  mean  variation  in  cydue  to  varying  A r'  for  the 
SUPERDUCK  data  was  negligible. 

E.  ERROR  ANALYSIS 

An  error  analysis  is  conducted  to  assign  appropriate  error  bounds  on  the  resultant 
spatially-dependent  Cy  values.  Approximated  errors  are  the  use  of  a  linear  wave  theory 
transfer  function  to  integrate  the  velocity  covariance  measurement  over  the  water  col¬ 
umn  and  errors  in  stationarity  due  to  tidal  variation.  Three  sources  of  error  which  are 
specifically  quantified  are  incident  wave  angle  error,  propagation  of  measurement  error, 
and  finite  differencing  error  for  the  radiation  stress  gradient. 

1.  Incident  Wave  Angle  Error 

Radiation  stress  (equation  3)  is  sensitive  to  the  incident  wave  angle  (a).  The 
determination  of  an  accurate  a  is  dependent  on  several  other  angles.  The  Zeiss  elec¬ 
tronic  surveying  system  was  required  to  measure  the  sled  yaw  angle.  The  Zeiss  error 
(Table  14)  is  assumed  negligible  in  the  a  determination.  The  adjustments  (Chapter  V) 
to  translate  the  sled  instrument's  three-dimensional  coordinates  into  the  FRF  coordinate 
system,  and  the  measurements  by  the  portable  Zeiss  system  which  made  them  possible, 
are  considered  to  have  negligible  error.  Three  angular  sources  of  error  which  cannot  be 
neglected  are  beach  angle  orientation  error,  current  meter  orientation  error,  and  current 
meter  directional  gain  error. 

a.  Beach  Angle  Orientation  Error 

An  analysis  is  conducted  to  assess  sensitivity  to  beach  angle  orientation  er¬ 
ror.  The  beach  angle  orientation  is  calculated  from  a  mean  bathymetric  contour  deter¬ 
mined  in  a  least  squares  sense  as  described  in  Chapter  V.  For  this  analysis, 
monochromatic,  linear,  unidirectional  wind  waves  are  assumed  for  a  straight-and- 
parallel  planar  beach.  Using  Snell's  Law,  the  wave  incident  angle  (a2)  at  the  inshore 
current  meter  may  be  determined  given  the  wave  incident  angle  (a,)  at  the  offshore  cur¬ 
rent  meter: 


where  wave  celerity  (C)  is  determined  from  equation  (70).  The  radiation  stress  gradient 
(from  equation  3)  in  the  cross-shore  direction  may  be  stated  as  a  finite  difference  be¬ 
tween  radiation  stress  calculated  at  the  two  on-offshore  current  meters  where  the  sub- 
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Figure  24.  tyand  percent  breaking  wave  (Q)  values  as  a  function  of  offshore  distance 
(October  15-18,  1986).  {cj  is  above  and  Q  is  below  the  profile  line). 


91 


script  "rear"  refers  to  the  rear  or  offshore  current  meter  and  "front"  to  the  front  or 
inshore  current  meter  of  the  sled: 


8Syx  lEn  Sin(a  ±  Ao.e)  cos(a  ±  Aae)]rear  -  [£n  sin  (a  ±  Aac)  cos(a  ±  Aa,)3from 

— —  *  - 7 -  99) 

cx  Ax 

where  A»e  is  the  beach  angle  orientation  error.  Total  energy  (E)  is  calculated  from: 

E  —  jPgH?™  (100) 

Root-mean-square  wave  height  (Hrms)  for  monochromatic  shallow  water  waves  can  be 
related  to  depth  (h)  in  the  inner  surf  zone  by  equation  (39).  For  SL'PERDLCK,  ratios 
of  Htms  to  h  ranged  from  0(0.3)  outside  the  surf  zone  to  0(0.4)  inside  the  surf  zone. 

The  radiation  stress  gradient  (equation  99)  versus  orientation  error  for  rep¬ 
resentative  depth,  slope,  and  incident  wave  angles  is  illustrated  in  Figure  25.  The  three 
angles  plotted  for  each  depth  are  centered  about  a  typical  incident  wave  angle  encount¬ 
ered  during  SL'PERDL'CK  for  that  particular  depth.  For  example,  an  incident  wave 
angle  of  25°,  a  depth  of  3.5  m,  and  a  beach  orientation  angle  error  of  -5°,  would  result 
in  an  overestimation  of  the  radiation  stress  gradient  by  10%.  On  the  other  hand,  a  ra¬ 
diation  stress  gradient  measured  in  a  shallower  depth  with  a  reduced  incident  wave  angle 
(h  =  1.5  m,  a  =  15°)  and  the  same  error  angle  of -5°  would  result  in  an  overestimation 
by  30%.  Note  that  the  overall  governing  parameter  is  incident  wave  angle.  Therefore, 
for  a  fixed  beach  orientation  angle  error,  the  relative  radiation  stress  error  increases  with 
decreasing  wave  angle,  and  hence  with  decreasing  depth.  For  SL'PERDL'CK,  errors  in 
mean  beach  orientation  angle  are  estimated  to  be  less  than  ±2°  and  1.0  m  was  the 
minimum  depth  in  which  data  were  acquired.  Therefore  radiation  stress  gradient  errors 
for  3.5,  2.5,  and  1.5  m  depths  would  be  +  5%,  ±  7-8%,  and  ±  10%,  respectively.  A  10% 
radiation  stress  error  can  be  equated  to  approximately  a  10%  cf  error,  as  long  as  wind 
stress  and  temporal  variability  of  mean  momentum  are  small.  As  the  magnitude  of  wind 
stress  and  temporal  variability  increase,  the  effect  of  a  fixed  radiation  stress  error  on  cj 
becomes  less.  Varying  beach  slope  from  0.01  to  0.10  or  varying  Aw  front  2  m  ("local" 
method)  to  25  m  ("integrated"  method)  had  negligible  effect  on  the  relative  radiation 
stress  gradient  error  due  to  an  incorrect  beach  orientation  angle.  Therefore  cy  error  due 
to  inaccurate  beach  orientation  angle  is  estimated  to  be  less  than  10%. 
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Figure  25.  Sensitivity  of  dSyx  /  cx  to  beach  angle  orientation  error  for  vary  ing  inci¬ 
dent  wave  angles  and  depths  and  tan  /?  =  0.03. 
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b.  Current  Meter  Orientation  Error 

Orientation  of  the  current  meters  axes  to  the  sled  frame  axes  was  considered 
to  be  within  1  —  2°  (Chapter  V).  Assuming  conservation  of  the  radiation  stress  gradient 
outside  the  surf  zone,  all  current  meters  were  individually  oriented  by  a  previously  de¬ 
scribed  selection  process  (Table  6).  This  numerical  rotation  is  considered  to  have  elim¬ 
inated  any  current  meter  orientation  error. 

To  assess  the  effect  of  other  rotations  which  minimized  the  radiation  stress 
gradient  and  which  had  realistic  spectra  shapes,  rotation  subsets  within  Table  6  combi¬ 
nations  were  also  examined.  For  example, 

•  in  run  16-2,  rotating  both  current  meters  by  a  lesser  amount  instead  of  just  the  rear 
current  meter  by  -2.86°,  or 

•  in  run  15-1.  rotating  the  front  current  meter  by  0.27°  instead  of  the  rear  current 
meter  by  -0.29°, 

yielded  no  substantial  difference  in  the  magnitude  of  Cj -  along  that  transect.  In  other 
words,  by  using  the  concept  of  conservation  of  Syx  to  initially  orient  the  current  meters, 
the  resultant  cj s  were  insensitive  to  the  manner  in  which  the  orientation  was  done. 

c.  Current  Meter  Directional  Gain  Error 

Current  meter  directional  gain  error  is  the  error  introduced  by  using  a  single 
gain  and  bias  factor  to  calibrate  current  flows  for  all  directions.  The  current  meters  were 
initially  calibrated  in  a  laboratory  flow  tank  where  gains  and  biases  were  determined  for 
meter  probe  angles  of  0°,  90°,  180°,  and  270°.  These  angles  correspond  to  the  four 
sensor  probes  protruding  from  the  current  meter  ball.  Using  the  gain  and  bias  factors 
for  the  u  and  v  sensors  determined  for  these  four  probe  angles,  the  current  meter  was 
then  checked  against  a  controlled  flow  of  0(0.5  ms"1)  for  meter  probe  angle  increments 
of  15°  .  The  standard  deviations  from  the  controlled  flow  for  sled  current  meters  1,  2, 
and  4  were  0.012,  0.012,  and  0.009  ms-1,  respectively.  This  error  is  included  within  the 
measurement  error  for  current  velocities  addressed  in  the  next  section. 

2.  Measurement  Error 

To  ascertain  the  effect  of  direct  measurement  error  on  the  final  determination 
of  cp  an  error  propagation  analysis  is  conducted  for  both  the  "local"  and  "integrated" 
methods  for  Cy  determination.  Standard  deviations  (o)  of  directly  and  indirectly  meas¬ 
ured  variables  used  in  equation  (93)  were  estimated  or  determined  from  laboratory  cali¬ 
bration  results  (Table  14). 

The  coefficient  of  variation  (Barry',  1978),  or  fractional  standard  deviation 
(Young,  1962),  is  defined  as: 
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Table  14.  MEASUREMENT  ACCURACY 


Direct  Measurement 

parameter 

q_ 

V*£ol 

tvpical  value 

water  density 

±  2  kgm-3 

±0.2 

1022  kgm-3 

Zeiss  vert  control1 

±  0.03  m 

±  1.5 

2  m 

Zeiss  horiz  control1 

±  0.03  m 

+  0.075 

40  m 

c  m  velocity 

+  0.05  ms-1 

±  5.0 

u  =  0.25  ms-i 

time 

wind  stress 

insignificant  error 
±  0.04  Jm~3 

±  40.0 

v  =  1.  ms-1 

0.1  Jrrr3 

MWL  (Paro) 

±  0.0003  m 

+  0.015 

2  m 

wind  speed 

±  0.35  ms-1 

±  5.0 

7  ms-1 

atmos.  press. 

+  2  mb 

+  0.2 

1010.0  mb 

atmos.  temp. 

±2°K 

±  0.7 

288°K 

Indirect  Measurement 

atmos.  drag  coeff. 

+0.5  x  10-3 

±  30. 

1.5  x  10~3 

atmos.  density 

±  0.008  kgm-3 

±  0.6 

1.25  kgm-3 

•  where  ±  refers  to 

an  approximate  standard  deviation. 

1  from  Birkemeier  and  Mason  (1984) 

yx= 


(101) 


where  A'  is  the  arithmetic  mean  of  replicate  measurements  of  a  variable.  For  the  coef¬ 
ficient  of  variation  determination  (Table  14),  the  strict  X  has  been  replaced  by  typical 
(order  of  magnitude)  values,  which  are  listed. 

Calibration  of  the  sled  current  meters  was  performed  before  and  after  the  ex¬ 
periment  by  two  different  calibration  laboratories.  Pre-calibration  gains  and  biases  are 
listed  in  Table  15.  Post-calibration  gains  varied  by  ±  3.0%  and  0.01  -  0.04  ms-'  in  bias 
or  offset.  A  conservative  estimate  of  measurement  standard  deviation  for  the  current 
velocities  is  ±  5%  to  include  the  unmeasured  environmental  source  of  error  between  field 
and  laboratory  conditions  as  well  as  the  uncertainty  (+  3%)  associated  with  using  a 
single  gain  factor  for  all  frequencies  (Guza  and  Thornton.  1980)  and  directions. 


The  general  rules  of  measurement  error  propagation  can  be  stated  for  four 


arithematic  operators  as  (Barry,  197S): 

,  /  2,2,2,  .2 
^sumordiff  =  ±  \/aA  +  aB  +  aC  +  -  +  an 

(102) 

product  -  ±  ABC  ...  +  (^)2  +  i^f  +  ., 

+ 

T  n  > 

(103) 

_  _  J_A_  !(  °A  \2  1  /  aB  ,2 

a  quotient  ~  —  g  \  '  A  ^  ^  ft  ’ 

(104) 

CTpower  =  i  O A 

(105) 

where  the  standard  deviations  aA  aB  etc.,  refer  to  the  ±  error  of  the  parameters  (e.g., 
A,  B)  being  analyzed  (Young.  1962). 


Table  15.  SLED  CURRENT  METER  CALIBRATION 


CM 
Ser.  H 

Sled 

Lctn  Ret 

Sled 
CM  # 

PCM 

Channel 

Vel 

Gain  Bias 

Correlation 
Coeff.  (r2) 

SI  14 

rear  mid 

1 

6 

u 

0.951  0.000 

0.99944 

5 

V 

0.916  -0.008 

0.99970 

S762 

front  lower 

4 

13 

u 

1.156  0.015 

0.99S99 

12 

V 

1.147  -0.012 

0.99S7S 

S354 

sideupper 

2 

8 

u 

1.457  -0.030 

0.99949 

7 

V 

1.457  0.017 

0.99955 

•  current  meters  were  Marsh-McBirnev,  Model  55 J 

•  calibrations  by  U.S.G.S,  Gulf  Coast  Hydroscience  Center,  NSTL,  MS 

•  formulation:  velocity  in  ms-1  =  (gain  x  voltage)  +  bias 


The  wind  stress  measurement  error  has  the  largest  Vx  of  all  the  parameters  in 
Table  14.  Vx  can  reach  as  high  as  40%,  depending  on  the  manner  and  conditions  under 
which  the  measurement  was  taken.  Conditions  during  SUPERDUCK  were  conducive 
to  accurate  measurements  of  wind  stress.  These  conditions  were: 

•  over-water  measurements  taken  during  onshore  winds, 

•  at  the  end  of  a  500  m  pier, 
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•  at  a  height  of  18.7  m  MSL  which  is  11.1  m  above  the  pier  deck. 

•  using  the  eddy  correlation  method,  and 

•  using  a  Campbell  CA27T  Sonic  Anemometer  as  backup  for  the  Gill  anemometer. 


Standard  deviation  of  atmospheric  density  ( pa )  is  estimated  from  equations 
(103),  (104),  and  the  equation  of  state  for  a  perfect  gas  (Haltiner  and  Martin,  1957). 
Therefore,  using  a  subjectively-determined  and  conservative  value  of  ±  30%  for  the  wind 
stress  measurement  error,  a  Q  coefficient  of  variation  of  30%  is  determined  from 
equations  (11)  and  (103). 

Rewriting  equation  (93)  as  : 


c/  = 


r 

dx 

p\ut\v 


dMy 

ct 

P I U,  |  V 


(106) 


or 

cf  =  -A  +  B  -  C  ,  (107) 


the  measurement  error  in  Cj,  or  the  standard  deviation  of  cf,  is  the  result  of  the  additive 
o's  of  the  three  terms  on  the  right-hand-side  of  equation  (106),  abbreviated  as  A,  B,  C 
in  equation  (107).  Applying  equations  (102)  to  (104),  results  in  standard  deviations  for 
the  four  terms  in  equation  (107)  as  listed  in  Table  16. 

Standard  deviation  determination  of  the  radiation  stress  gradient  (oA),  requires 
o  of  (c'SyJ /  dx),  which  in  turn  requires  o  of  Sy[.  Sy[  or  Syx  arc  calculated  from  equation 
(7S)  using  the  directly  measured  parameters  u',  v',  p,  and  h  .  An  upper  bound  for  o  of 
SyJ  can  be  formulated  by  approximating  for  shallow  water  waves: 


r  r J2 

Syx  =  ph  I  Cuvdf  =  *phcovuv  <,  phouov  (108) 

where  covuv  is  the  covariance  of  the  u  and  v  velocities  and  is  bounded  by  (e.g.  Davis, 
1986): 

covuy  <.  Ou<7y  (109) 

Therefore  cs  t  can  conservatively  be  approximated  as: 
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(110) 


p 


*5  T 


~  + 


s  T 

°yx 


VfTT*'  +  +  <X)! 


(T^  and  are  negligible  compared  to  <xu  and  av  (Table  14)  and  may  be  neglected, 
One  may  then  write: 


A5„ 


=  ±  AS 


r 


\/5S 


2 

rear 


,  7\2 
l^sJ/ronr 


dll) 


a&sj_ 

Ax 


AS/ 

A.y 


)2  +  (^f)2 


(112) 


The  standard  deviation  of  A*  is  much  less  than  the  standard  deviation  of  A Sy[  and 
therefore  may  be  neglected.  The  standard  deviation  determined  by  equation  (112)  for 
the  "integrated"  method  is  reduced  by  a  factor  of  one-half  to  account  for  that  method's 
increased  sample  size  (see  equation  S4).  Standard  deviation  of  the  radiation  stress  gra¬ 
dient  for  each  of  the  sled's  data  runs,  using  both  the  "local"  and  "integrated"  methods 
for  Cf  determination,  were  calculated.  The  mean  a  for  radiation  stress  gradient  of  all  four 
transects  is  listed  in  Table  16. 

The  standard  deviation  for  the  wind  stress  (rj)  was  determined  from  equations 
(11)  and  (103),  and  Table  14.  Standard  deviation  of  6Myt  ct  was  determined  by  as¬ 
suming: 


d.\tv  A.Vv  A  (pDV) 

—Tf  *  -7—  *  - - -  (113) 

ct  At  At  K  ’ 

since  \ly!  4  My,  and  then  applying  equations  (102)  through  (104)  and  Table  14  values. 

A  typical  mean  Cj-  value  of  0.003  was  used  in  equation  (101)  to  determine  the  Cf 
coefficient  of  variation  in  Table  16.  Because  of  the  relatively  large  incident  wave  angles 
and  the  measurement  precision  utilized  in  the  SUPERDLCK  experiment,  propagation 
of  measurement  error  resulted  in  reasonable  Cj  values  for  both  methods  (Table  16).  The 
"integrated"  method's  oC/  was  superior  to  the  "local"  method  because  of  the  larger  rela¬ 
tive  error  for  the  radiation  stress  gradient  (equation  1 12)  experienced  when  A.v  is  small. 
The  integrated  method's  error  bounds  for  the  cy  range  in  Table  13  are  +  4-15%. 
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Table  16.  ERROR  PROPAGATION  ANALYSIS 


Term  (eqn  107) 

o_ 

j)  Cf  Method 

A 

0.0015 

24 

local 

0.00013 

10 

integrated 

B 

0.00003 

31 

C 

0.00007 

70 

cf 

0.0015 

50 

local 

0.00015 

5 

integrated 

where  term 

U. 

C>.> 

& 

* 

T 1  dx 

g 

e\Llci 

r 

pl 

U,\v 

p  1 W,  1  V 

P 1  u,  |  V 

3.  Finite  Differencing  Error 

Thornton  and  Guza  (19S1)  examined  the  error  involved  in  representing  an  ar¬ 
bitrary  continuous  function  by  discrete  measurement  points.  Using  a  Taylor  series  ex¬ 
pansion  for  Syx: 


Syx(x  +  Ax)  =  SYX(x)  +  Ax 


?Syx(x)  (Ax)2  d2Svx(x) 


cx 


+ 


■>! 


Ax 


+ 


(114) 


and  then  solving  for  the  true  radiation  stress  gradient: 

A.r2  v  “yx 


cS,r  A  Syx  a  v.  <?2S„r  d3Sv 


yx _ 

cx  Ax 


Ax 

21 


cx2 


3! 


(115) 


cx 


or 


cSyx  AS^X 

— —  —  — - -  =  Error 

cx  Ax 


(116) 

I 


where  the  sum  of  the  second  and  infinite  terms  on  the  right-hand  side  of  equation  (115) 
is  the  error  of  the  finite  difference  space  approximation. 

Relative  error  (R)  may  be  defined  as: 


I 


L 

I 

I 

I 


» 


R  = 


Error 


dx 


(H7) 


Solving  for  R  where  only  the  first  error  term  in  equation  (115)  is  considered,  Syx  for 
monochromatic  shallow  water  linear  waves  inside  the  surf  zone  may  be  expressed  as: 

Syx  =  £  sin  a  cos  a  * /ID5'12  (118) 

where 

E  -  P5tD\ 

V  =  K mA 

D  =  >/  +  h,  and 

A  =  complicated  constant 


Substituting  equation  (118)  into  (117): 


R* 


l 

4 


Ax 


tan  (1 
h 


(119) 


where  D  «  h.  (Thornton  and  Guza.  1981) 

This  approximation  (equation  119)  is  adequate  providing  Ax  and  tan  /?  are  rea¬ 
sonably  small.  However,  Ax  and  tan  0  are  not  always  small  in  this  application.  There¬ 
fore  R  is  redefined  to  include  a  sufficient  number  of  terms  on  the  right-hand  side  of 
equation  (115)  to  allow  convergence  to  within  5%  of  R.  Resultant  computations  re¬ 
vealed  that  equation  (119)  was  sufficiently  accurate  (within  5%)  for  most  of  the 
SL'PERDL'CK  data.  However,  using  the  "integrated"  method  between  the  most  inshore 
data  run  and  the  beach  required  additional  iterations  (up  to  five)  to  stabilize  R. 

The  relative  error  in  equation  (117)  is  biased  in  one  direction  because  the  term 
of  largest  magnitude  is  written  as: 


Ax 


c2S 


>x 


2! 


cx 


oc 


iH 

ch 

cx 


=  slope 


(120) 
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and  assuming  negligible  bottom  curvature.  Therefore  the  sign  of  R  is  fixed  as  the  same 
sign  as  beach  slope.  The  finite  approximation  for  the  radiation  stress  gradient  may  be 
written  as: 


ASyjf 

Ax 


(1  +  Ri- 


dS* 


cx 


(121) 


where  R  represents  the  overestimation  of  dSyx  /  dx  for  positive  slopes  and  the  underes¬ 
timation  for  negative  slopes.  Finite  difference  relative  error  for  the  four  transects  of  this 
research  are  listed  in  Table  17  and  Table  18.  They  are  listed  for  both  methods  of  co¬ 
determination.  The  larger  error  for  the  "integrated"  method  value  between  the  most  in¬ 
shore  data  run  and  the  beach  was  due  to  the  combined  effect  of  increased  Ax,  decreased 
depth,  and  steeper  slope,  all  which  combine  to  increase  R.  The  relative  error  is  signif¬ 
icantly  reduced  when  using  the  "local"  method  due  to  its  much  smaller  A .r.  "Integrated" 
method  s  R  ranges  from  -27%  on  the  nearshore  side  of  the  bar  to  nearly  200%  on  the 
nearshore  slope.  However,  except  for  the  gross  error  on  the  nearshore  slope,  'inte¬ 
grated"  method's  relative  error  is  less  than  50%.  This  would  cause  c,  error  due  to  finite 
differencing  to  be  less  than  50". o  (assuming  negligible  wind  stress  and  steady  state  con¬ 
ditions  in  equation  93). 

4.  Combined  Error 

The  various  sources  of  error  are  now  combined  into  a  total,  or  combined,  error 
relative  to  the  bar.  To  establish  a  conservative  combined  error,  finite  differencing  error 
is  assigned  a  range  of  error  rather  than  a  bias  and  unquantified  errors  are  estimated  at 
±  20%.  Mean  finite  difference  error  for  the  trough  does  not  include  the  gross  error, 
0(200"  o),  between  the  most  inshore  sled  run  position  and  the  beach.  Therefore  the  final 
Cf  value  for  the  trough  does  not  reflect  the  nearshore  slope. 

The  maximum  combined  error  is  determined  by  adding  the  four  maximum  pos¬ 
sible  errors  listed  in  Table  19.  However  it  is  not  reasonable  to  define  maximum  error 
in  this  manner  since  each  of  the  four  errors  have  a  range  of  uncertainty  associated  with 
them  and  they  are  uncorrelated.  It  is  very  unlikely  that  the  errors  would  all  be  a  maxi¬ 
mum  and  in  the  same  direction  at  the  same  time  (Barry,  1978).  Hence,  assuming  the 
errors  are  independent  of  one  another,  the  combined  error  is  obtained  from  equation  102 
and  results  in  final  momentum  balance  Cj  values  of  0.004  ±  0.0010.  for  offshore  the  bar, 
0.002  ±  0.0006  for  on  ton  and  immediately  before  the  bar,  and  0.001  ±  0.0003  for  the 
trough  region. 
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Table  17.  FINITE  DIFFERENCE  RELATIVE  ERRORS 


October  15  1986  Relative  Error  (R) 

Sled  Run 

zlr  (ml 

tan  /? 

MWL  (m) 

R(%) 

15-1 

1.99 

0.0199 

2.80 

1. 

27.99 

0.0238 

2.29 

23. 

15-2 

1.99 

0.0287 

1.78 

2. 

24.24 

0.0241 

1.72 

27. 

15-5 

1.99 

-0.0031 

1.67 

0. 

44.65 

-0.0127 

1.99 

-20. 

15-4 

1.99 

0.0554 

2.30 

4. 

42.03 

0.0549 

1.15 

197. 

Beach 

October  16  1986  Relative  Error  (R) 

16-2 

1.99 

0.0158 

3.33 

1. 

50.52 

0.0290 

2.61 

46. 

16-3 

1.99 

0.0349 

1.89 

3. 

17.69 

0.0341 

1.70 

28. 

16-4 

1.99 

0.0309 

1.51 

3. 

12.99 

0.0198 

1.52 

13. 

16-5 

1.99 

0.0145 

1.54 

1. 

33.49 

-0.0081 

1.72 

-12. 

16-6 

1.99 

0.0029 

1.91 

0. 

33.88 

0.0564 

0.96 

196. 

Beach 

Table  18.  FINITE  DIFFERENCE  RELATIVE  ERRORS 


October  17  1986  Relative  Error  (R) 

Sled  Run 

Ax{m} 

tan  /? 

MWL  (m) 

R  (%) 

17-1 

1.99 

0.0182 

3.85 

1. 

34.60 

0.0224 

3.22 

19. 

17-2 

1.99 

0.0397 

2.59 

2. 

26.80 

0.0364 

2.01 

*  39. 

17-3 

1.99 

0.0309 

1.43 

3. 

9.54 

0.0265 

1.31 

15. 

17-4 

1.99 

0.0221 

1.19 

3. 

4.09 

0.0232 

1.14 

6. 

17-5 

1.99 

0.0166 

1.09 

2. 

31.56 

-0.0153 

1.25 

-27. 

17-6 

1.99 

-0.0036 

1.40 

0. 

4.88 

0.0109 

1.48 

J. 

17-7 

1.99 

0.0137 

1.55 

1. 

6.27 

0.0276 

1.49 

9. 

17-8 

1.99 

0.0495 

1.43 

5. 

20.06 

0.0716 

0.72 

197. 

Beach 

October  18  1986  Relative  Error  (R) 

18-1 

1.99 

0.0113 

3.80 

0. 

4S.33 

0.0259 

3.09 

32. 

18-2 

1.99 

0.0405 

2.38 

3. 

IS. 85 

0.0369 

1.94 

29. 

18-3 

1.99 

0.0334 

1.49 

J. 

21.60 

0.025S 

1.31 

34. 

18-4 

1.99 

0.0083 

1.13 

i. 

27.67 

-0.0096 

1.30 

-15. 

18-5 

1.99 

-0.0223 

1.4S 

*  A-  • 

33.71 

0.0438 

0.74 

196. 

Beach 
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Table  19.  COMBINED  ERROR  RELATIVE  TO  THE  BAR  ON  THE  BED 
SHEAR  STRESS  COEFFICIENT  AS  DETERMINED  FROM  A  MO¬ 
MENTUM  BALANCE 


Trough 

On  the  Bar 

Offshore 

beach  angle  orientation  error 

±  8% 

±  10% 

±  10% 

measurement  error 

±  5% 

±  5% 

+  5% 

finite  differencing  error 

mean  R 

ASlt 

-  error  =  1  /( 1  +  R) 

-19% 

19% 

31°/0 

23% 

-16% 

-24% 

unquantified  errors  (est.) 

±  20% 

±  20°  o 

±  20°.  o 

combined  ey  error 

±  32% 

±  28% 

±  33°/o 

final  Cf 

0.001  ±  0.0003 

0.002  ±  0.0006 

0.004  ±  0.0010 
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VIII.  LONGSHORE  CURRENT  MODEL 


A.  INTRODUCTION 

A  first-order  longshore  current  model  for  a  barred  beach  is  proposed.  The  model  is 
based  on  balancing  the  wave-induced  radiation  stress  gradient  with  longshore  bed  shear 
stress,  turbulent  momentum  exchange,  and  longshore  surface  wind  stress.  Stationary 
wave  conditions  and  straight  and  parallel  bottom  contours  are  assumed.  The  waves  are 
described  as  random,  but  narrow-banded  in  frequency  and  direction.  Model  results  are 
compared  to  field  observations  and  bed  shear  stress  coefficients  are  determined  by 
model-fitting. 

B.  MODEL  THEORY 

There  are  no  field-verified  longshore  current  models  for  a  barred  beach.  Therefore 
a  model  which  has  been  field-verified  for  a  planar  beach  was  selected  for  this  analysis. 
The  model  (Thornton  and  Guza,  1986;  hereafter  referred  to  as  TG86)  is  based  on  a 
simple  balance  between  the  gradient  of  the  wave-induced  radiation  stress,  the  bed  shear 
stress,  and  the  gradient  of  turbulence-induced  radiation  stress.  The  model's  assumptions 
are  stated  in  the  previous  paragraph.  The  longshore  current  model  is  based  on  a  random 
wave  height  transformation  model  (Thornton  and  Guza,  1983;  hereafter  referred  to  as 
TG83).  Modifications  to  the  TG86  model  are  made  based  on  the  availability  of  addi¬ 
tional  data  from  SUPERDUCK  and  with  the  purpose  of  improving  and  broadening  the 
model  s  applicability  to  barred  beaches.  These  modifications  are  the  inclusion  of  surface 
wind  stress  in  the  momentum  balance,  retention  of  incident  wave  angle  dependency,  and 
the  use  of  a  breaking  wave  height  distribution  based  on  a  slope-dependent  y  and  an  im¬ 
proved  empirical  expression  for  the  breaking  wave  percentage  as  a  function  of  H,  h,  and 
y.  Whereas  TG86  found  turbulent  momentum  exchange  to  be  unimportant  for  their 
planar  beach  topography,  turbulent  momentum  exchange  is  shown  to  be  required  in 
order  to  model  current  velocities  in  the  trough  region  of  a  barred  beach.  The  TG86 
longshore  current  model  is  briefly  summarized  below  with  modifications  described  later. 

1.  Wave  Height  Transformation 

This  portion  of  the  model  is  an  extension  of  work  by  Battjes  and  Jannscn  (1978). 
The  wave  height  transformation  is  conceptualized  as  an  energy  flux  balance  given  by: 
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<£„> 


c(£Cg  cos  at) 
c.v 


n::j 


where  Cg  cos  a  is  the  shoreward  component  of  wave  group  velocity,  o.  is  the  mean  inci¬ 
dent  wave  angle  of  the  narrow-banded  waves,  and  <  eb  >  is  the  ensemble-averaged  dis¬ 
sipation  due  to  wave  breaking.  Analysis  of  SL'PERDUCK  data,  as  well  as  results  from 
TG83,  indicated  that  the  ensemble-averaged  dissipation  due  to  bottom  friction  <  cy->  is 
negligible  compared  with  <  cb  >,  and  therefore  <  tj>  is  omitted  in  equation  (122). 

Based  on  analysis  of  field  data  from  Torrey  Pines  Beach,  CA,  the  Rayleigh  wave 
height  probability  density  function  (pdf).  /»(//),  was  shown  to  qualitatively  describe  the 
measured  wave  heights  across  that  planar  beach.  The  Rayleigh  pdf  is  described  by: 


where  H  is  wave  height.  The  Rayleigh  distribution  describes  all  waves. 

Only  the  breaking  waves  are  dissipative,  resulting  in  changes  in  the  momentum 
flux,  i.e.  wave  forcing.  Therefore  a  distribution  describing  the  breaking  waves  alone  is 
required.  Wave  height  probability  density  functions  were  calculated  from  field  meas¬ 
urements  taken  at  Ft.  Ord's  Soldiers  Beach,  Monterey.  CA.  These  pdfs  were  used  to 
more  realistically  model  random  wave  height  transformation  across  the  surf  zone. 

Within  the  surf  zone,  breaking  wave  height  distribution  C/’tl/F)J  is  approxi¬ 
mated  by  weighting  the  Rayleigh  pdf: 

Po(H)  =  W(H)  p{H)  (124j 

where  the  weighting  function  W(H)  <  1.  it  is  noted  the  breaking  ware  distribution  is 
not  a  pdf.  but  that  the  area  of  pn{1T)  gives  the  percentage  of  the  waves  that  are  breaking. 
A  simple  form  of  W(H)  which  approximates  the  fraction  of  waves  that  break  from  on¬ 
shore  to  saturation  conditions  is: 

W( H )  =  |^  ~~j~  (125) 


where 


(120) 


10b 


for  saturated  wave  breaking  conditions.  Based  on  the  Ft.  Ord  data,  a  modification  to 
equation  (1251  which  more  heavily  weights  the  larger  waves  is  given  by: 


\\ 


<h>-[W  [-+[*?! 


^  i 


Wave  breaking  dissipation  is  modeled  as  a  simple  periodic  bore: 


? 

=  4  PS 


(BHf 


(128) 


where  B  is  a  breaker  coefficient  of  0(1)  and  f  is  the  mean  wave  frequency  which  is  the 
peak  frequency  for  narrow-band  wave  spectra.  The  B  coefficient  accounts  for  differences 
in  various  breaker  types  and  is  a  function  of  the  proportion  of  foam  region  on  the 
breaker  face  (Figure  26).  The  peak  frequency  in  the  SUPERDLCK  Syx  spectra  is  used 
for  /. 


Figure  26.  Periodic  bore  used  to  describe  spilling  breakers,  (from  Thornton  and 
Guza,  19S3). 


The  ensemble-averaged  dissipation  due  to  wave  breaking  is  calculated  by 
multiplying  the  dissipation  for  a  single  broken  wave  of  height  H  by  the  probability  of 
wave  breaking  at  each  height  [pfc(//)D.  and  integrating  over  all  H: 
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B 1  f 
<t»>  =  —psj 


/•<*> 


H3pt>(H)  dH 


1 1 29) 


By  substituting  from  equations  (124)  and  (127)  and  evaluating  the  definite  integral, 
equation  (129)  may  be  written  as: 


<£i)>  = 


•>V* 

16 


pgBTf 


H5 

2.  3 

y  h 


l 


l 


(1  +  (Hrms/yh)2)512 


(130) 


Energy  flux  (ECg  from  equation  122),  is  approximated  using  linear  theory  lor  C„  at  fre¬ 
quency/,  and  energy  is  calculated  based  on  the  Rayleigh  pdf: 


E  =  y  prC  pgH;ms  (131) 

Substituting  equations  (130)  and  (131)  into  (122)  yields  a  first-order  dilferential  equation 
for  the  linear  wave  height  transformation  model: 

77  -j  wrfms  cs cos «  =  77  1 1 32) 

y  ” 

A  spatially-dependent  wave  height  [//rw(.r)3  can  then  be  determined  for  an  arbitrary 
bottom  profile  across  the  surf  zone  by  numerical  integration,  given  an  initial  value  of 
Hrms.  The  Hrm j  profile  is  used  as  input  for  the  longshore  current  calculation. 

2.  Longshore  Current 

The  longshore  current  distribution  across  the  surf  zone  is  based  on  conservation 
of  momentum.  Rewriting  equation  (3).  wave-induced  radiation  stress  may  be  written: 

,  A 

p  r~ i — '  _  ■*'  sin  a  / 1 1  ^  -v 

S)X  =  ECg  cos  a — — —  (133) 

where 


i 


i 


4 


i 


sin  a 
C 


constant 


sin 

Co 


(134.) 


by  Snell's  Law  of  linear  wave  refraction.  The  subscript  "o"  refers  to  any  arbitrary  lo¬ 
cation.  Therefore  the  wave-induced  radiation  stress  gradient  may  be  written  as 


IDS 


C. X 


Syx  ~ 


sin  y.,. 


ox 


(EC,  cos  a)  = 


sin  yr 


<cb> 


(135) 


where  <  cb  >  is  substituted  from  equation  (122).  Assuming  steady  state  conditions  and 
negligible  wind  stress,  the  longshore  momentum  balance  is  written  as: 


sin  y_o 

Q 


<£t> 


(136) 


TGS6  used  equation  (136)  and  model-fitting  to  solve  for  the 'horizontal  profile 
of  longshore  current  across  a  planar  beach.  They  obtained  a  hierarchy  of  solutions  by 
choosing  combinations  of  linear  or  non-linear  bed  shear  stress,  with  and  without  turbu¬ 
lent  momentum  exchange.  Turbulent  momentum  exchange  was  found  not  to  be  im¬ 
portant  for  their  planar  beach  topography.  They  concluded  that  a  non-linear  bed  shear 
stress  without  turbulent  momentum  exchange  was  their  preferred  choice  for  general  ap¬ 
plications,  based  on  its  fit  to  their  data.  A  mean  c,-  for  their  four  days  of  data  was  0.006 
±  0.0007. 

C.  MODIFICATION  OF  MODEL  BASED  ON  FIELD  DATA 

Modifications  to  the  TGS6  longshore  current  model  are  the  inclusion  of  surface 
wind  stress,  incident  wave  angle  dependency,  and  use  of  a  breaking  wave  weighting 
function  based  on  a  slope-dependent  ;  and  an  improved  empirical  expression  for  break¬ 
ing  wave  percentage  as  a  function  of  H.  h,  and  y.  The  use  of  turbulent  momentum  ex¬ 
change  is  found  to  be  important  for  modeling  longshore  current  velocities  in  the  trough. 
The  inclusion  of  surface  wind  stress  and  incident  wave  angle  dependency  are  covered  in 
the  Model  Results  section.  The  other  modifications  are  addressed  next. 

1.  A  Slope  Dependent  y 

Since  the  TGS3  model  was  developed  for  planar  beaches,  y  was  assigned  a  con¬ 
stant  value  of  0.42  (equation  39)  and  is  the  ratio  Hrmslh  for  saturated  wave  breaking 
conditions.  However,  Sallengcr  and  Holman  (1985)  found  >■  to  be  a  function  of  beach 
slope  based  on  data  from  several  field  experiments.  They  determined  an  empirical  re¬ 
lation  for  a  slope-dependent  y  based  on  their  data: 

y  =  3.2  tan/?  +  0.30  (137) 


The  ratio  Hr„.Jh  is  compared  to  beach  slope  (Figure  27)  for  both  the  SUPER- 
DUCK  data  and  the  data  used  by  Sallenger  and  Holman  (19S5).  The  solid  line  in 
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Figure  27  is  Sallenger  and  Holman's  (1985)  empirical  relation  for  y.  Most  of  the 
SUPERDUCK  data  has  the  same  magnitude  of  scatter  as  previous  data  sets  about  the 
solid  line,  with  even  negative  slopes  (found  on  the  inshore  side  of  the  bar)  in  general 
agreement.  The  exception  is  for  SUPERDUCK  data  sets  located  in  the  trough  and  on 
the  foreshore  slope  (Runs  11-2,  15-4,  17-6,  17-7,  17-8,  21-4).  They  have  a  significantly 
lower  value  of  Hmslh  and  are  indicative  of  non-saturated  breaking  wave  conditions. 
Obviously  a  slope-dependent  y  relationship  is  required  for  a  barred  beach.  Based  on 
Figure  27  and  the  general  fit  of  the  SUPERDUCK  data,  equation  (137)  is  used  for  the 
SUPERDUCK  analysis. 

2.  Breaking  Wave  Weighting  Function 

The  weighting  function,  W(H),  used  in  TG83  was  based  on  only  four  data  runs 
for  a  single  transect  of  a  beach  with  a  tidal  plateau.  To  improve  on  the  formulation  of 
W(H),  breaking  wave  pdfs  are  calculated  for  the  SUPERDUCK  data.  Time  series  of 
pressure  head,  [ph{t)~\,  are  linearly  transformed  into  >j(i)  time  series  by  first  calculating 
the  complex  Fourier  spectra  of  ph(t),  then  multiplying  by  a  linear  transform  function 
[H(f)l  band  pass  filtering  between  0.03  -  1.00  Hz,  and  finally  applying  an  inverse 
Fourier  transform  to  the  >f(f)  function. 

The  transfer  function  H(f)  (from  equation  3S)  is  held  constant  at  the  0.5  Hz 
value  for  the  interval  0.5  Hz  to  1.0  Hz  so  as  to  include  the  additional  higher  frequency 
variance  for  the  Hrns  calculations,  but  not  to  overly  bias  them  with  the  exponentially 
increasing  H(f)  values.  Filtering  is  accomplished  by  zeroing  out  the  Fourier  amplitude 
coefficients  of  the  undesired  frequencies  before  applying  the  inverse  Fourier  transform. 
The  entire  34.1  min  record  is  transformed  at  one  time.  As  an  example,  the  initial  three 
minutes  of  Run  15-2’s  ph{t)  time  series  is  compared  with  its  transformed  >/(/)  series  in 
Figure  2S.  The  transfer  function  H<f)  and  the  G^(/)  spectra  determined  from  the  entire 
tj(t)  series  are  also  included.  Wave  periods  and  heights  are  determined  from  the  >7 (/)  time 
series  by  the  down-crossing  method.  Root-mean-square  wave  height  is  determined  sta¬ 
tistically  from  the  calculated  ti(t )  time  series  as: 

,v 

Hms  “  [jZ^],/3  (,38) 

/«  1 

The  statistical  Hrms  from  the  rj{t)  series  are  compared  with  the  energy-based  Hrms 
(equation  90)  determined  from  both  the  Gn(J)  spectra  discussed  in  Chapter  VH  and  the 
Gv(f)  spectra  calculated  from  the  rj(t)  record  (Figure  28).  Root-mean-square  wave 


Figure  27.  Ratio  HrmJh  versus  beach  slope  for  various  data  sets.  Circled  data 
represent  trough  and  foreshore  slope  locations. 
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height  from  all  three  sources  compared  favorably  [within  (10%)],  with  the  difference 
indicating  a  measure  of  non-linearity  in  the  surface  elevation  record. 

As  discussed  in  Chapter  IV,  visual  observations  of  wave  dissipation  (i.e. ,  wave 
breaking)  at  the  sled  were  electronically  marked  on  the  data  set  records.  These  marks 
are  then  superimposed  on  the  calculated  t](t)  time  series  to  identify  which  waves  were 
breaking  during  each  34.1  min  data  run.  From  this  procedure,  wave  heights  and  periods 
for  all  waves,  and  a  subset  for  just  breaking  waves,  are  obtained.  Pdfs  of  all  waves  and 
breaking  waves  for  15-17  October  are  illustrated  in  Figure  29  and  Appendix  A 
(Figure  47  through  Figure  49).  Since  the  data  were  measured  sequentially  and  not  si¬ 
multaneously,  wave  heights  are  normalized  by  the  mean  offshore  wave  height  to  account 
for  any  variability  in  the  incident  wave  conditions.  The  pdf  s  are  normalized  to  the  mean 
offshore  wave  height  measured  by  a  WAVERIDER  buoy  located  6  km  offshore. 
Equipment  malfunction  precluded  breaking  wave  identification  on  other  days  of  the  ex¬ 
periment. 

Breaking  wave  data  from  SUPERDU'CK,  as  well  as  the  four  data  runs  from  Ft. 
Ord  (TG83),  are  compared  in  Figure  30.  The  ordinate  is  the  fraction  of  waves  which 
are  breaking  (Q)  and  the  abcissa  is  the  ratio,  ]IrrrJyh,  where  y  is  slope-dependent  and 
calculated  from  equation  (137).  Thus  for  saturated  conditions,  both  the  ordinate  and 
abcissa  are  equal  to  unity.  The  breaking  wave  fraction  is  obtained  by  dividing  the 
number  of  waves  identified  as  breaking  by  the  total  number  of  waves.  The  total  number 
of  waves  in  a  data  run  was  approximately  400.  The  Ft.  Ord  data  are  indicated  by  an  's’, 
SUPERDUCK  data  acquired  in  the  trough  are  solid  circles,  and  all  other  SU  PERDU CK 
data  are  indicated  by  open  cirlcles.  Three  breaking  wave  weighting  functions  are  illus¬ 
trated.  The  dotted  line  is  the  TG83  weighting  function  (equation  127).  The  chain-dot 
line  is  from  a  theoretical  expression  (Battjes  and  Janssen,  1978;  Battjes  and  Stive.  19S5) 
where  the  cumulative  probability  distribution  of  all  wave  heights  (breaking  or  non¬ 
breaking)  is  the  Rayleigh  distribution,  but  abruptly  terminated  at  the  breaking  wave 
height: 

Q  “  ["T/t]1"2  +  1  (139) 

Neither  of  these  two  functions  provide  a  reasonable  fit  to  the  data.  However,  a  function 
which  substantially  improves  the  fit  to  the  data  is: 


112 


~  2 


LJ 

5  o 

u~> 

U~> 

U 

2  -C  -1 


TIME  (min) 


2n 


RUN  1 5 —  2  filtered  0.03  ~  l.OOhz 

^rms  C.73 

W  0  ^vJ-J\J -Jy kJ'-J Sj'^p \jArJ*\yh'v 

P 


-2  J 


13 


TIME  (min) 


Figure  28.  E.xample  ph[t)  and  ij(t)  time  series  »itli  corresponding  GJJ)  spectra. 


Figure  29.  Wave  height  and  breaking  nave  height  distributions  (hatched  area)  nor¬ 


malized  to  offshore  Hrmi  for  16  October.  The  Rayleigh  distribution  is 
the  solid  line  and  empirical  breaking  wave  distributions  given  by 
equations  (127)  and  (140)  are  the  dotted  and  dashed  lines,  respectively. 
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and  is  delineated  by  a  dashed  line.  This  function  (equation  140)  has  no  physical  justi¬ 
fication  but  is  merely  an  empirical  expression  derived  to  better  fit  the  non-trough  data. 
It  is  noted  that  none  of  the  three  functions  fit  both  the  trough  and  non-trough  data. 
The  right  half  of  the  right-hand  side  of  equation  (140)  is  taken  from  equation  (127)  and 
is  used  to  more  heavily  weight  the  larger  waves  towards  breaking.  Observations  show 
that  the  largest  waves  at  a  particular  location  are  more  likely  to  break  (TG83).  This 
skewing  of  the  pdf  is  substantiated  by  the  breaking  wave  observations  (hatched  areas) 
seen  in  the  pdf  figures.  Very  little  wave  breaking  occurs  until  hrmslyh  reaches  0.8,  then 
the  percentage  of  breaking  waves  increases  rapidly.  The  addition  of  the  SLPERDl  CK 
data  suggests  a  more  rapid  increase  in  the  breaking  wave  percentage  within  a  smaller 
Hrmslyh  band  -  almost  to  the  point  of  describing  a  "step"  function  indicative  of  all  waves 
breaking  at  a  fixed  breaker  point. 

The  ratio  HrmJyh  in  Figure  30  has  a  certain  amount  of  error  associated  with  it. 
The  error  in  y  is  the  result  of  an  imperfect  fit  to  the  data  in  Figure  27.  Tidal  fluctuations 
also  affect  the  h  term.  Thus  there  are  runs  which  have  Hrmslyh  ratios  greater  than  one, 
vet  were  unsaturated. 

It  is  also  noted  that  during  field  identification  of  the  breaking  waves,  breaking 
wave  percentages  were  visually  estimated  to  be  ol'tcn  in  the  range  between  60-90%.  Yet 
subsequent  data  analysis  reveals  the  percentages  to  be  much  lower.  This  may  perhaps 
be  similar  to  the  phenomena  associated  with  visually  estimating  wave  heights  from  ships, 
where  it  has  long  been  recognized  that  the  human  eye  perceives  an  average  of  the  highest 
1  3  wave  heights  and  not  the  true  mean  wave  height.  This  lead  to  the  designation  of  ship 
wave  height  observations  as  "significant  wave  heights"  and  not  mean  wave  heights.  In 
a  similar  manner,  SLPERDUCK  visual  observations  of  breaking  wave  percentages  were 
overestimates.  Thus  visual  observations  of  saturated  wave  conditions  from  other  ex¬ 
periments  may  need  to  be  reviewed. 

D.  MODEL  RESULTS 

1.  Wave  Height  Transformation 

As  input  to  the  wave  transformation  model,  IIrms  and  f  measured  during  the 
orientation  runs  outside  the  surf  zone  (Runs  15-1,  16-2,  17-1,  18-1)  are  used  as  initial 
conditions.  Mean  wave  frequency  is  assigned  as  the  frequency  of  the  spectral  peak  in 
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Figure  30.  Breaking  nave  fraction  (Q)  versus  the  ratio  Hrmsjyh. 
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sea  surface  elevation.  Mean  incident  wave  angle  is  then  assigned  as  the  angle  for/ from 
the  significant  wave  spectra  for  that  day.  TG83  and  TGS6  assumed  small  wave  incident 
angles  and  set  the  cos  a  term  in  equation  (122)  to  unity.  Because  of  the  large  incident 
wave  angles  during  SL'PERDL'CK.  incident  wave  angle  dependency  is  retained 
throughout  the  analysis  by  inclusion  of  the  cos  y.  term.  A  forward  stepping  scheme  is 
used  to  evaluate  equation  (122): 

{ECg  cos  a)2  =  (ECg  cos  a),  +  (  <  >  )1A.v  (141) 

All  variables  in  equation  (141)  are  known  except  B,  which  is  determined  by  model-fitting 
the  calculated  Hrms  to  the  field-measured  Hrms.  Model-fitting  is  accomplished  by  mini¬ 
mizing  the  least  square  error  between  observed  and  modeled  data.  Results  of  the 
model-fitting  are  shown  in  Figure  31  and  Appendix  A  (Figure  50  through  Figure  52). 
Field  Hrn:s  is  represented  by  an  "x"  for  each  of  the  sled  runs  and  is  determined  from 
pressure  sensor  measurements.  The  spatial  step  size.  A.v.  was  varied  from  0.5  m  to  10.0 
m  with  negligible  effect  on  the  Hrms  profile.  Therefore  A.v  =  5.0  m  was  arbitrarily  se¬ 
lected  to  allow  suitable  cross-shore  definition  along  with  reduced  computational  time. 

A  constant  B  coefficient  of  0(1)  is  determined  by  a  least  square  fit  for  all  four 
days  (Table  20).  Field  and  modeled  Urrns  show  good  agreement  except  for  the  trough 
area  on  15  and  16  October.  Root-incan-square  wave  height  is  undennodeled  by  about 
35°  o  for  these  exceptions.  Spatially-variable  B  coefficients  were  also  determined  which 
gave  exact  fits  to  the  field  data.  However  their  drastic  spatial  changes  in  B  appeared 
unrealistic  since  the  breaker  type  was  observed  to  be  similar  for  most  of  the  transect, 
except  on  the  very  nearshore  slope  inshore  of  the  last  sled  runs.  A  constant  B  coefficient 
for  a  transect  was  deemed  appropriate. 

2.  Longshore  Current  Model  Results 

Longshore  current  is  modeled  by  including  surface  wind  ..'•ess  in  the  momentum 
balance  equation  of  equation  (136): 


<cb>  =  - 


J>  +  - 

T>-  +  lv 


—  V 

c.v 


and  using  various  combinations  of  the  following  approaches: 

•  linear  or  non-linear  formulation  of  r\., 

•  with  and  without  turbulent  mixing. 

•  spatially-variable  or  constant  c,  and  N. 


Figure  31.  Model-generated  Hrms  and  V  compared  v\ith  field  observations  versus 
distance  for  invariant  cj  and  B  and  including  turbulent  momentum  ex¬ 
change,  16  October.  Model  results  are  solid  lines  and  field  observa¬ 
tions  are  "x"'s.  Elevation  is  referenced  to  MSL. 


In  all  cases,  the  spatially-in variant  B  coefficients  listed  in  Table  20  are  used.  As  with  the 
B  coefficient  determination,  Cy  and  N  are  model-fitted  by  a  least  square  error  criteria. 

A  spatially-invariant  Cj  with  a  non-linear  (and  presumably  more  realistic)  t*  and 
without  turbulent  mixing  results  in  poor  agreement  between  modeled  and  field-measured 
longshore  current  velocities  (Figure  32  and  Appendix  A,  Figure  53).  Although  the 
general  shape  of  the  observed  longshore  current  profile  is  modeled,  there  are  large  dis¬ 
crepancies  between  observed  and  modeled  current  velocities,  especially  in  the  trough  re¬ 
gions.  Current  velocities  in  the  trough  region  are  in  error  by  as  much  as  a  factor  of  7. 
Order  of  magnitude  for  cy  is  similar  to  the  momentum  balance  determination  of  cy 
(Chapter  VII)  for  all  4  transects. 

Allowing  ty  to  vary  spatially  provides  some  improvement  in  the  model-fit. 
However  major  discrepancies  still  occur  on  the  nearshore  slope  where  a  steep  increase 
in  velocity  is  modeled  along  with  significant  underestimation  of  velocity  in  the  trough. 
The  spatially-variant  cy  solution  is  dismissed  because  of  these  discrepancies. 

Although  TG86  determined  turbulent  momentum  exchange  to  be  unimportant 
for  their  planar  beach  data,  turbulent  momentum  exchange  is  investigated  for  this  anal¬ 
ysis  as  a  possible  mechanism  to  transfer  momentum  across  the  bar  and  into  the  more 
placid  trough  region  where  wave  breaking  is  reduced.  Turbulent  momentum  exchange 
is  expressed  as  the  last  term  in  equation  (142)  and  is  approximated  using  an  eddy- 
viscosity  formulation: 


V 


J-*  puw'dz 


(M3) 


where  D  is  total  depth  and  the  eddy  viscosity  coefficient  (c)  is  parameterized  by  a 
dimensionless  and  adjustable  coefficient  (N)  from  equation  (21).  A  first  approach  for 
the  inclusion  of  turbulent  momentum  exchange  was  to  determine  the  minimum  least 
square  error  for  spatially-invariant  Cy  and  N  across  a  reasonable  range  of  Cy  and  N  values. 
Bed  shear  stress  coefficients  were  varied  from  0.0005  to  0.0200  and  N  was  varied  from 
0.000  to  0.016.  A  three-dimensional  depiction  of  their  solution  surface  is  illustrated  in 
Figure  33  for  a  typical  data  run.  Minimum  least  square  error  is  represented  by  a  peak 
in  the  solution  surface  closest  to  the  top  of  the  three-dimensional  box.  Note  that  there 
are  many  combinations  of  Cy  and  N  that  minimize  the  least  square  error  to  nearly  the 
same  degree.  Thus  too  many  degrees  of  freedom  result  in  multiple  mathematical  sol¬ 
utions  with  a  corresponding  loss  of  insight  for  the  physics  involved. 
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Figure  32.  Model-generated  V  compared  to  field  observations  versus  offshore  dis¬ 
tance  for  constant  cy  no  turbulent  momentum  exchange,  15-16 
October.  Model  results  are  solid  line  and  field  observations  are  "x'"s. 
Error  is  mean  least  square  error. 


Example  of  a  least  square  error  solution  surface  for  model-fitted 
longshore  current  for  various  cj-  and  N  combinations.  Results  show  an 
almost  constant  ridge  for  optimal  cj  and  N  indicating  a  non-unique 
solution. 


MODELED  L/S  CURRENT  ERROR  FI5  F(cr,  N) 


The  spatially-invariant  cf  approach  is  therefore  considered  a  baseline  upon 
which  to  use  turbulent  momentum  exchange  to  improve  the  model  results.  The  proce¬ 
dures  are  to  first  determine  an  invariant  cy,  and  then  to  vary  X  until  the  least  square  error 
is  minimized.  This  results  in  a  substantial  improvement  in  model  results  over  previous 
approaches.  The  18  October  data  run  required  an  additional  Cy  iteration  after  a  value 
of  X  was  determined.  The  model  result  prior  to  this  additional  iteration  is  represented 
by  a  dashed  line  in  Figure  52.  With  longshore  momentum  being  transferred  over  the 
bar  and  into  the  trough,  the  previous  velocity  underestimation  in  the  trough  has  been 
significantly  reduced.  Velocities  are  modeled  within  20%  of  measured  values,  with  the 
exception  of  the  trough  region  on  17  October.  This  case  is  unique  and  is  addressed  be¬ 
low.  Interval  bars  in  Figure  31  and  Appendix  A  (Figure  50  through  Figure  52)  repre¬ 
sent  10%  and  20%  of  the  observed  value  and  are  provided  solely  as  a  visual  aid  for 
interpreting  the  data.  Mean  least  square  error  ranges  between  7-15%  for  all  four 
transects. 

Current  velocities  in  the  trough  (Runs  17-6,  17-7,  and  17-8)  on  17  October  are 
much  lower  (and  actually  of  opposite  direction  to  the  longshore  current  offshore  the  bar) 
than  on  the  other  three  days.  Guza  et  al.  (1986)  showed  that  it  is  SyJ,  and  not  the 
structure  of  Syx(f),  that  is  important  for  longshore  current  response.  It  is  therefore  pos¬ 
sible  that  the  low  current  velocities  in  the  trough  on  17  October  are  the  result  of  bi¬ 
directional  wave  forcing  where  breaking  swell  and  sea  are  driving  longshore  currents  in 
opposing  directions  and  cancelling  each  other  out.  Radiation  stress  spectra  (Appendix 
F)  show  weak  wave  forcing  for  the  trough  region  on  17  October.  Radiation  stress 
spectra  for  runs  17-1  and  17-2  show  sea  and  swell  approaching  from  different  directions 
with  relatively  large  incident  angles  (18°  and  —9°,  respectively).  The  cross-over  wave 
period  separating  swell  from  sea  is  6.7  seconds.  A  scatter  plot  of  breaking  wave  heights 
versus  breaking  wave  periods  for  a  trough  run  (17-6)  shows  that  approximately  half  of 
the  breaking  waves  consisted  of  swell  and  half  consisted  of  sea  (Figure  34).  A  scatter 
plot  for  Run  17-4  (located  offshore  the  bar)  does  not  show  as  much  breaking  wave  ac¬ 
tivity  for  the  swell.  Thus  the  minimum  combined  velocity  in  the  trough  on  17  October 
and  the  poor  model  results  for  that  particular  day  and  location  are  possibly  a  result  of 
opposing  longshore  currents. 

3.  Bed  Shear  Stress  Coefficients 

Spatially-invariant  Cy's  are  determined  from  model-fitting  using  a  non-linear 
formulation  of  t*  (TG86)  and  turbulent  momentum  exchange  (Table  20).  The  ratio  of 
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Figure  34.  Scatter  diagram  of  breaking  wave  periods  and  heights  for  17 

October.  Run  17-4  is  offshore  the  bar  and  17-6  is  in  the  trough. 
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Table  20.  MODEL-FITTING  RESULTS  FOR  B,  cf  AND  N 


B 

X 

Error 

cf\ 
c ? 


15  Oct 

16  Oct 

0.74 

0.78 

0.006 

0.015 

0.146 

0.071 

0.003 

0.003 

0.003 

0.002 

17  Oct 

18  Oct 

U0 

1.06 

0.004 

0.016 

0.078 

0.114 

0.009 

0.003 

0.004 

0.002 

Error  is  mean  least  square  error  for  the  transect 
'transect  mean  Cy  determined  by  model-fitting 


2 

transect  mean  cy  determined  by  momentum  balance 


wind  force  to  wave  force  for  the  SUPERDL'CK  data  was  0.09  (see  Chapter  VI 1),  there¬ 
fore  the  inclusion  of  surface  wind  stress  in  this  model  did  not  have  significant  impact  on 
the  magnitude  of  Cy  from  model-fitting  for  this  particular  data  set.  However,  as  discussed 
earlier,  the  inclusion  of  surface  wind  stress  under  certain  wind,  wave,  and  beach  slope 
conditions  can  have  significant  impact.  Bed  shear  stress  coefficient  values  for  15-18 
October  are  0.003,  0.003,  0.009,  and  0.003,  respectively.  X  coefficient  values  are  0.006, 
0.015,  0.004,  and  0.016,  respectively.  These  Cy  values  are  in  good  agreement  (see 
Table  20)  with  the  mean  c/ s  determined  for  each  day  by  the  momentum  balance 
method.  This  agreement  is  encouraging  and  lends  credibility  to  the  ty  results  considering 
the  complex  physics  being  addressed  and  the  different  and  independent  methods  em¬ 
ployed. 
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IX.  SUMMARY  AND  CONCLUSIONS 

Previous  investigators  of  longshore  currents  have  used  many  simplifying  assump¬ 
tions  which  often  preclude  meaningful  and  realistic  results.  The  assumptions  of  a  planar 
beach,  a  steady  and  depth  uniform  flow,  invariant  bed  shear  stress,  and  negligible  wind 
stress  are  quantitatively  examined  by  measuring  the  relative  importance  of  each  term  in 
a  local  longshore  momentum  balance  based  on  measurements  acquired  during  a  field 
experiment  at  a  barred  beach.  Wind  and  wave  forcing  of  longshore  currents  is  investi¬ 
gated  with  a  numerical  model  and  from  field  measurements.  Measurements  of  pressure, 
current  velocities  at  three  levels,  wind  speed,  and  sea  surface  slope  were  acquired  from 
an  instrumented  sled  positioned  at  various  locations  across  a  transect  of  a  barred  beach. 
Unfortunately  the  slope  measurements  were  found  during  data  analysis  to  be  unreliable 
and  were  disregarded. 

Three-dimensional  depictions  of  the  current  structure  over  a  nearshore  bar  are  de¬ 
veloped.  The  cross-shore  velocities  [0(O.2ms-1)3  indicate  a  relative  onshore  flow  in  the 
upper  third  of  the  water  column  and  a  relative  offshore  flow  in  the  lower  third  of  the 
water  column,  irrespective  of  bar  location  and  in  agreement  with  cross-shore  theory. 
The  longshore  velocities  [O(1.0ms_I)3  indicate  nearly  depth  uniform  flow  with  a  slight 
velocity  increase  with  elevation.  This  depth  uniform  flow  confirms  the  appropriateness 
of  the  depth-integration  aspect  of  radiation  stress  for  a  barred  beach.  The  longshore 
vertical  current  structure  did  not  exhibit  a  logarithmic  profile. 

A  numerical  model  is  developed  to  compare  wind  and  wave  forcing  for  vary  ing  wind 
conditions  and  beach  slope.  The  finite-depth  TMA  spectral  wave  model  is  used  as  input 
for  the  energy  spectrum.  Waves  are  assumed  to  be  locally  generated  and  fully  arisen. 
W'ind  force  is  proportional  to  surf  zone  width  and  therefore  inversely  proportional  to 
beach  slope.  Wave  force  is  independent  of  surf  zone  width  and  beach  slope.  Wind  force 
is  found  to  increase  more  rapidly  than  wave  force  for  an  increase  in  wind  speed  or  an 
increase  in  the  wind  and  wave  incident  angle  relative  to  beach  normal.  Ratios  of  wind 
to  wave  force  are  calculated.  For  a  given  wind  direction,  a  wind  speed  increase  of  10  ms-1 
results  in  a  50%  increase  in  this  ratio.  In  other  words,  wind  force  becomes  increasing 
significant  relative  to  wave  force,  as  wind  speed  and  incident  wind  direction  increase. 
For  obliquely  onshore  winds  and  a  steep  beach  slope  (tan  /?  =  0.10),  wind  force  is  gen¬ 
erally  less  than  10%  of  the  wave  force  and  thus  may  be  considered  negligible.  For  ob- 
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liquely  onshore  winds  and  a  gentle  beach  slope  (tan  /?  =  0.01).  wind  force  can  be 
significant,  with  wind  to  wave  force  ratios  ranging  from  15%  to  100%  for  winds  between 
10-30  ms-1  and  wind  and  wave  directions  from  10—  70°.  Thus  under  certain  wind  and 
wave  conditions,  wind  force  can  be  a  first  order  term  in  the  longshore  momentum  bal¬ 
ance.  For  offshore  winds,  longshore  currents  are  totally  wind-driven.  In  summary,  the 
model  shows  that  wind  force  should  be  considered  as  a  forcing  term  for  longshore  cur¬ 
rents  and  that  it  is  inappropriate  to  assume  that  wind  force  can  be  omitted. 

There  are  no  field  studies  addressing  surface  wind  stress  for  robust  surf  zones  (Hm0 
greater  than  C.5  m).  Surf  zone  atmospheric  drag  coefficients  are  determined  from  surface 
udnd  stress  measurements  acquired  immediately  outside  the  surf  zone  and  wind  speed 
measurements  acquired  by  the  sled  anemometer  inside  the  surf  zone.  Stability- 
dependent  drag  coefficients  are  calculated  using  the  Busingcr  et  al.  (1971)  iterative  ap¬ 
proach.  A  conservative  increase  of  33%  for  Q  was  made  to  account  for  the  significantly 
increased  surface  roughness  of  the  surf  zone  due  to  breaking  waves  based  on  the  results 
of  Davidson  et  al.  (19S8)  which  indicate  an  increased  Q  for  chaotic  seas.  During  periods 
of  atmospheric  convective  instability  and  reduced  surface  roughness  (low  wave  heights), 
Hsu’s  (1970)  empirical  relationship  for  Cd  is  used.  Resultant  Q's  range  from  1.5  x  10-3 
to  3.4xl0~3  for  non-convective  instability  and  from  4.6  x  10~3  to  8.25  x  10-3  for 
convective  instability.  These  Q  magnitudes  are  consistent  with  previous  field  measure¬ 
ments  of  Cd  over  the  open  ocean  as  well  as  the  theoretical  Cd  values  of  Geernaert  et  al. 
(1987)  for  shallow  water  depths. 

All  terms  in  the  longshore  momentum  balance  are  calculated  from  field  measure¬ 
ments  of  pressure,  current  velocity,  and  wind  speed.  The  sled  was  initially  positioned 
beyond  the  surf  zone  for  a  complete  data  run.  This  additional  data  run  allowed  a  daily 
numerical  calibration  of  the  current  meters  during  data  analysis  based  on  conservation 
of  dSyx/dx  outside  the  surf  zone.  This  procedure  is  essential  for  obtaining  realistic 
physical  results  and  is  a  distinguishing  feature  of  the  research.  Longshore  momentum 
due  to  unsteady  flow  was  determined  to  be  less  than  10%  of  the  longshore  momentum 
due  to  steady  flow  for  this  data  set.  Relative  contributions  of  the  radiation  stress  gra¬ 
dient  ( 3Syx  I  8x),  temporal  term  (< ?My/ct),  and  surface  wind  stress  (rp  to  the  longshore 
momentum  balance  were  calculated  and  analyzed  by  location  relative  to  the  bar.  Mean 
relative  contributions  for  the  wave,  wind,  and  temporal  terms  were  0.86,  0.08,  and  0.06, 
respectively,  for  both  offshore  and  on  top  the  bar.  However  in  the  trough,  the  mean 
relative  contributions  were  significantly  different  (0.64,  0.17,  and  0.19,  respectively),  l  or 
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the  trough  region,  the  relative  importance  of  wave  forcing  is  reduced  by  25%,  whereas 
the  relative  importance  of  surface  wind  stress  is  doubled  and  the  temporal  term  relative 
importance  is  tripled.  All  three  terms  are  of  first  order  in  the  trough.  Although  surface 
wind  stress  was  nearly  constant  across  the  surf  zone,  the  relative  contribution  of  surface 
wind  stress  increases  shoreward  from  the  location  of  maximum  breaking  due  to  the  de¬ 
crease  in  wave  forcing.  The  relative  contribution  of  the  temporal  term  appears  to  be 
largest  in  the  trough,  possibly  due  to  the  presence  of  surf  beat  or  edge  waves.  Wind  to 
wave  force  ratios  calculated  from  the  field  data  (0.09)  are  in  agreement  with  those  de¬ 
termined  from  the  model  (0.12)  for  the  same  wind  conditions. 

Spatially-dependent  bed  shear  stress  coefficients  are  determined  between  sled  run 
locations.  The  spatial  step  size  (At  of  the  radiation  stress  gradient)  is  determined  by 
error  analysis  to  be  too  small  when  assigned  the  2  m  distance  between  the  current  meters 
on  the  sled.  It  had  originally  been  planned  that  the  slope  array,  w’hich  did  not  work 
properly,  would  have  provided  the  necessary  precision  to  measure  the  radiation  stress 
gradient  over  the  2  m  distance.  Reasonable  <y  results  are  calculated  when  the  spatial  step 
size  is  increased  to  the  distance  between  run  positions  0(20  -  40  m).  Resultant  c/s  range 
from  less  than  0.001  to  0.006.  Spatially- variable  mean  cj s  are  found  with  values  of  0.004 
±  0.0010  for  offshore  the  bar,  0.002  ±  0.0006  for  on  top  the  bar,  and  0.001  ±  0.0003  for 
the  trough  (excluding  the  nearshore  slope).  A  trend  for  a  shoreward  decreasing  Cy  is 
noted,  possibly  due  to  the  shoreward  decrease  in  wave  action.  A  mean  cf  for  the 
nearshore  slope  is  not  specified  since  error  analysis  shows  that  a  significant  ovcrcsti- 
mation  of  Cy  can  occur  for  the  nearshore  slope  due  to  the  combined  eflect  of  decreased 
depth  and  steeper  slope  upon  the  finite  difference  approximation.  The  percentage 
change  in  cy  due  to  neglecting  surface  wind  stress  of  the  temporal  term  in  the  momentum 
balance  ranges  from  near  zero  to  100%.  with  a  mean  change  between  10-30%.  Thus  for 
this  data  set  and  under  certain  conditions,  the  surface  wind  stress  term  and  the  temporal 
term  can  be  of  first  order  in  the  momentum  balance. 

Longshore  current  velocities  at  a  barred  beach  are  modeled  to  within  ±  20%  of  the 
observed  values  at  SL'PERDUCK  by  modifying  the  Thornton  and  Guza  (1986)  model 
to  include  w'ind  stress,  wave  incident  angle  dependency,  slope-dependent  y,  and  an  im¬ 
proved  breaking  wave  weighting  function.  Mean  least  square  error  between  modeled 
and  observed  current  velocity  ranges  from  7-15%.  The  inclusion  of  turbulent  momen¬ 
tum  exchange  to  transport  momentum  into  the  more  placid  trough  region  is  required 
and  is  physically  reasonable.  Spatially-invariant  cy  values  determined  by  model-fitting 


for  15-18  October  are  0.003,  0.003,  0.009,  and  0.003,  respectively.  These  cj  values  agree 
well  with  mean  c}-  values  determined  for  each  day  by  the  momentum  balance  method. 
This  agreement  is  encouraging  and  lends  credibility  to  the  cj  results  considering  the 
complex  physics  being  addressed  and  the  different  and  independent  methods  employed. 
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APPENDIX  A. 


DATA  ANALYSIS  FIGURES  FOR  DAYS  OTHER  THAN 
OCTOBER  16,  1986 
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36.  Mean  current  structure  on  12  October  1986.  Numbers  represent  data 
runs,  short  dotted  lines  are  MWL's,  and  distance  alongshore  is  arbi¬ 
trary. 
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;  37.  Mean  current  structure  on  15  October  1986.  Numbers  represent  data 
runs,  short  dotted  lines  are  MWL's,  and  distance  alongshore  is  arbi¬ 
trary. 
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38.  Mean  current  structure  on  17  October  1986.  Numbers  represent  data 
runs,  short  dotted  lines  are  MWL's,  and  distance  alongshore  is  arbi- 


Figure  39.  Mean  current  structure  on  18  October  1986.  Numbers  represent  data 
runs,  short  dotted  lines  are  MWL's,  and  distance  alongshore  is  arbi¬ 


trary. 
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Figure  40.  Mean  current  structure  on  21  October  1986.  Numbers  represent  data 
runs,  short  dotted  lines  are  MWL's,  and  distance  alongshore  is  arbi¬ 
trary. 


Figure  42.  Offshore  Syx  spectra  and  S}J  from  'puv'  gage  621  for  17  October 
1986.  Numbers  indicate  data  runs.  Intervals  are  95%  confidence. 


Figure  43.  Offshore  Syx  spectra  and  S>xr  from  'puv'  gage  621  for  18  October 
1986.  Numbers  indicate  data  runs.  Intervals  are  95%  confidence. 
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Figure  44,  Wind  speed,  wave  height,  atmospheric  stability,  and  initial  C,  data  for 
October  11-12,  1986.  Beach  is  assumed  parallel  to  the  y-axis.  Open 
arrowheads  denote  wind  direction  and  closed  arrowheads  denote  wave 
incident  angle.  Winds  arc  adjusted  to  z  =  10  m. 
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Figure  45.  Wind  speed,  wave  height,  atmospheric  stability,  and  initial  Cd  data  for 
October  17-18.  1986.  Beach  is  assumed  parallel  to  the  y-axis.  Open 
arrowheads  denote  wind  direction  and  closed  arrowheads  denote  wave 
incident  angle.  Winds  are  adjusted  to  z  =  10  m. 


140 


\\\  *  u 


Figure  46.  Wind  speed,  wave  height,  atmospheric  stability,  and  initial  Cj,  data  for 
October  21,  1986.  Beach  is  assumed  parallel  to  the  y-a.\is.  Open 
arrowheads  denote  wind  direction  and  closed  arrowheads  denote  wave 
incident  angle.  Winds  are  adjusted  to  z  =  10  m. 
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Figure  47.  Wave  height  and  breaking  wave  height  distributions  (hatched  area)  nor¬ 
malized  to  offshore  Hrms  for  15  October..  The  Rayleigh  distribution  is 
the  solid  line  and  empirical  breaking  wave  distributions  given  by 
equations  (127)  and  (140)  are  the  dashed  and  dotted  lines,  respectively. 
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Figure  48.  Wave  height  and  breaking  wave  height  distributions  (hatched  area)  nor¬ 
malized  to  offshore  Hrms  for  17  October.  The  Rayleigh  distribution  is 
the  solid  line  and  empirical  breaking  wave  distiibutions  given  by- 
equations  (127)  and  (140)  are  the  dashed  and  dotted  lines,  respectively. 
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Figure  49.  Wave  height  and  breaking  wave  height  distributions  (hatched  area)  nor¬ 
malized  to  offshore  Hrm  for  IS  October.  The  Rayleigh  distribution  is 
the  solid  line  and  empirical  breaking  wave  distributions  given  by 
equations  (127)  and  (140)  are  the  dashed  and  dotted  lines,  respectively. 
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Model-generated  Hrms  and  V  compared  with  field  observations  versus 
distance  for  invariant  c,  and  B  and  including  turbulent  momentum  ex¬ 
change,  15  October.  Model  results  are  solid  lines  and  field  observa¬ 
tions  are  "x"'s.  Elevation  is  referenced  to  MSL. 


Figure  51.  Model-generated  Hrms  and  V  compared  with  field  observations  versus 
distance  for  invariant  Cf  and  B  and  including  turbulent  momentum  ex¬ 
change,  17  October.  Model  results  are  solid  lines  and  field  observa¬ 
tions  are  "x"'s.  Elevation  is  referenced  to  MSL. 
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B=  1.06 


Figure  52.  Model-generated  Hrms  and  V  compared  with  field  observations  versus 
distance  for  invariant  cj  and  B  and  including  turbulent  momentum  ex¬ 
change,  18  October.  Model  results  are  solid  lines  and  field  observa¬ 
tions  are  ”x‘"s.  Elevation  is  referenced  to  MSL.  Model  result  prior  to 
the  last  Cy  iteration  is  represented  by  the  dashed  line. 
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tance  for  constant  <y  and  no  turbulent  momentum  exchange,  17-18  Oc¬ 
tober.  Model  results  are  solid  line  and  field  observations  are  "x  "s. 


Error  is  mean  least  square  error. 
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APPENDIX  B. 


CURRENT  VECTORS  OVERLAYED  ON  BATHYMETRY 
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Figure  54.  Sled  and  alongshore  array  mean  current  vectors  for  II  October  1986 
overlayed  on  bathymetry  .  Velocity  scale  is  at  upper  right. 
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Figure  55.  Sled  and  alongshore  array  mean  current  vectors  for  12  October  1986 


overlayed  on  bathymetry.  Velocity  scale  is  at  upper  right. 
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Figure  57.  Sled  and  alongshore  array  mean  current  vectors  for  16  October  1986 
overlayed  on  bathymetry  .  Velocity  scale  is  at  upper  right. 


overlayed  on  bathymetry.  Velocity  scale  is  at  upper  right. 


Figure  59.  Sled  and  alongshore  array  mean  current  vectors  for  18  October  1986 
overlayed  on  bathymetry.  Velocity  scale  is  at  upper  right. 


Figure  60.  Sled  and  alongshore  array  mean  current  vectors  for  21  October  1986 
overlayed  on  bathymetry.  Velocity  scale  is  at  upper  right. 
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APPENDIX  C.  ATMOSPHERIC  STABILITY  DETERMINATION 

The  following  set  of  equations  represent  an  iterative  procedure  to  determine  z  L 
based  on  the  stability  flux  profile  relation  of  Businger  et  al.  (1971). 


Initial  equations: 

Tjew  ~  Tair  4-  35(  ln(.0l  RH)) 
AT-  Tair  -  T surface  ~  -0098Zr 
Tfiet =  Tair  +  273.15 


Assign  all  -~'s  =  0 
Determine  V  from  iteration: 


T  «■ 


ATk 


T  =  2  ln[ 


1  +x2) 


] 


Y  =  -7(-jr) 


T<0 


Q  —  Qsi^dew  +  273. 1 5)(.01  RH) 
Qo  =  Qs(TSea  +  273.15) 


ft  =  .62510”  n  =  23.84  -  -  5.03  log (TKW) 

g«  (g-gp)* 

ln(^-)-V*(f- )q 

T*  -  7*  +  6.1  x  10-4(rar2‘) 


(Cl) 

(C2) 

(C3) 


(C4) 


(C5) 

(C6) 

(C7) 

(C8) 

(C9) 

(CIO) 
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Continue  iterations  from  equation  (C4)  until  ( ~ )  converges 


Definitions: 

zT  =  temperature  measurement  height  (m), 
zq  -  humidity  measurement  height  (m), 

RH  =  relative  humidity  (%), 

Td  =  dewpoint  temperature  (degrees  C) 
z0T  =  temperature  roughness  length  (2  x  10~6) 
Zoq  -  humidity  roughness  length  (2  x  10-6) 

'I'*  =  integral  diabatic  term  for  temperature 
'Vq  =  integral  diabatic  term  for  humidity 
Q  =  specific  humidity  (g  ka) 


Figure  61.  Surface  weather  analysis  for  0700  E.S.T.  J5  October  1986. 
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Figure  62.  Surface  weather  analysis  for  0700  E.S.T.  16  October  1986. 
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APPENDIX  E.  ANALYSIS  PARAMETERS  FOR  ALL  SLED  DATA  RUNS 


Run 

start 

time 

(hhmmss) 

h 

(m) 

^rms 

(m) 

fm 

(Hz) 

c  r 

•v 

{Jm-') 

EMW 

(deg) 

Zanem 

(m) 

tan  /? 

11-1 

123504 

2.99 

0.92 

0.0938 

24.33 

1.4 

5.92 

.000 

11-2 

134359 

2.80 

0.85 

0.0938 

18.86 

1.2 

6.11 

.000 

11-3 

143700 

1.95 

0.76 

0.0938 

18.51 

1.5 

6.95 

.106 

12-1 

074600 

2.97 

1.12 

0.0938 

-134.39 

-5.1 

5.95 

.059 

12-2 

092453 

2.03 

0.73 

0.0938 

13.23 

1.2 

6.88 

.000 

12-3 

100300 

2.29 

0.64 

0.0781 

-14.03 

-1.6 

6.64 

.000 

12-5 

145000 

2.93 

1.21 

0.07S1 

-56.10 

-1.8 

6.00 

.000 

12-6 

155700 

3.03 

0.93 

0.0781 

-28.95 

-1.6 

5.89 

.000 

12-7 

164200 

1.94 

0.82 

0.0781 

13.61 

0.9 

6.97 

.106 

15-1 

103257 

2.79 

0.92 

0.1563 

150.03 

9.1 

6.13 

.029 

15-2 

120000 

1.78 

0.76 

0.1563 

60.60 

5.2 

7.15 

.029 

15-3 

124700 

1.67 

0.57 

0.1563 

16.76 

2.5 

7.25 

.029 

15-4 

133100 

2.30 

0.59 

0.2031 

33.45 

5.0 

6.62 

.053 

16-2 

122200 

3.33 

0.98 

0.1875 

250.54 

14.7 

5.60 

.012 

16-3 

131821 

1.89 

O.SO 

0.1875 

151.02 

12.1 

7.03 

.034 

16-4 

140600 

1.51 

0,57 

0.1875 

57.20 

8.9 

7.41 

.034 

16-5 

144900 

1.54 

0.52 

0.1875 

47.99 

9.0 

7.39 

.015 

16-6 

153600 

1.91 

0.49 

0.1875 

25.94 

KKkXI 

7.01 

.005 

17-1 

092040 

3.85 

0.72 

0.1719 

61.53 

6.3 

5.08 

.017 

17-2 

100036 

2.59 

0.69 

0.1719 

52.84 

hes 

6.32 

.034 

17-3 

104900 

1.43 

0.57 

0.1094 

24.14 

7.50 

.034 

17-4 

113500 

1.18 

0.41 

0.1094 

27.96 

7.8 

7.74 

.034 

17-5 

122857 

1.09 

0.36 

0.1094 

70.37 

31.4 

7.83 

.034 

17-6 

134336 

1.41 

0.24 

0.0938 

1.45 

1.1 

7.52 

.017 

17-7 

141749 

1.55 

0.26 

0.1094 

3.99 

2.9 

7.3S 

.017 

17-8 

150024 

1.43 

0.32 

0.1094 

2.60 

1.2 

7.45 

.017 

18-1 

110730 

3.80 

0.93 

0.2031 

175.62 

11.8 

5.12 

.007 

18-2 

115000 

2.38 

0.84 

0.1250 

124.59 

8.6 

6.54 

.032 

18-3 

122555 

1.49 

0.67 

0.1094 

49.62 

5.3 

7.44 

.032 

18-4 

130900 

1.13 

0.38 

0.1094 

99.69 

45.0 

7.80 

.007 

18-5 

134937 

1.48 

0.28 

0.0469 

3.71 

2.2 

7.45 

-.025 

21-1 

110319 

3.55 

0.76 

0.0938 

-6.82 

-0.6 

5.31 

.021 

21-2 

121953 

1.65 

0.68 

0.0938 

2.38 

7.22 

.000 

21-3 

130448 

2.51 

0.49 

0.0938 

-17.72 

ra 

6.35 

-.029 

21-4 

134800 

2.33 

0.45 

0.0938 

2.08 

0.5 

6.54 

.014 

21-5 

144000 

2.17 

0.38 

0.0938 

11.54 

■ESI 

6.71 

.014 

21-6 

154200 

3.05 

0.76 

0.0938 

4.20 

0.3 

5.82 

.021 

21-7 

163059 

1.20 

0.59 

0.0938 

0.28 

0.0 

7.64 

.010 

21-8 

171000 

2.22 

0.45 

0.093S 

6.13 

1.5 

6.65 

-.029 

21-9 

174900 

1.51 

0.4S 

0.0781 

-7.72 

-1.6 

7.35 

.0S9 

•  start  time  is  referenced  to  Eastern  Standard  Time 

•  Hrms  is  measured  from  variance  of  the  pressure  sensor 

•  zanem  is  the  height  of  the  anemometer  above  the  mean  water  level  as  measured 
by  the  pressure  sensor  for  that  run 

•  x,y  coordinates  are  referenced  to  the  FRF  coordinate  system 

•  fm  is  the  peak  freq  in  the  sea  sfc  elevation  spectrum 

•  'EMW  is  the  "equivalent  monochromatic  wave"  incident  angle  (Thornton  and 
Guza,  19S6) 

•  Sy/  is  the  mean  Sy/  for  the  submerged  current  meters  on  that  run 

•  fm  is  the  peak  frequency  of  the  sea  surface  elevation  spectra 
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Run 

sled  x 
coord 
(m) 

sled  y 
coord 
(m) 

Palm 

(10 2  Pa) 

0 

(°  77V) 

H',« 

(ms-1) 

boa 

(deg) 

#  sub¬ 
merged 
c/m's 

11-1 

168.44 

899.90 

1023.6 

067. 

13.37 

0.5 

3 

11-2 

139.66 

905.65 

1023.8 

066. 

12.84 

0.5 

3 

11-3 

116.11 

907.30 

1023.9 

065. 

12.72 

0.5 

3 

12-1 

247.15 

900.44 

1024.1 

049. 

10.00 

0.5 

12-2 

19S.94 

904.57 

1023.8 

047. 

8.82 

0.5 

12-3 

154.69 

904.92 

1023.7 

047. 

8.71 

0.5 

12-5 

19S.55 

904.42 

1023.6 

052. 

5.66 

0.5 

12-6 

159.41 

902.77 

1023.5 

047. 

5.04 

0.5 

12-7 

125.21 

904.17 

1023.6 

041. 

■HD 

0.5 

15-1 

244.01 

887.14 

1017.3 

037. 

9.37 

4.0 

EBI 

mxmmw 

S91.67 

1017.2 

037. 

8.11 

4.0 

15-3 

191.78 

S95.64 

1017.2 

039. 

7.81 

4.0 

msn 

■BOI 

901.89 

1017.3 

040. 

7.80 

4.0 

16-2 

269.28 

1157.39 

1016.S 

030. 

12.28 

-0.5 

16-3 

218.76 

1162.25 

1016.9 

027. 

nmai 

-0.5 

16-4 

2<M.06 

1165.24 

nm 

10.90 

-0.5 

16-5 

1SS.07 

1 166.64 

1017.1 

026. 

10.22 

-0.5 

16-6 

154.58 

1170.37 

1017.3 

026. 

10.04 

-0.5 

17-1 

264.71 

II 

1019.1 

ool. 

-0.5 

17-2 

230.10 

1182.31 

DEfll 

276. 

7.24 

-0.5 

17-3 

203.30 

1184.07 

1019.3 

15o. 

7.06 

-0.5 

17-4 

193.76 

1 1 84.49 

1019.4 

o37. 

6.55 

-0.5 

17-5 

189.67 

1184.39 

1019.5 

033. 

5.89 

-0.5 

17-6 

158.10 

11SS.54 

1019.9 

025. 

5.53 

-0.5 

17-7 

153.23 

1 189.60 

1020.1 

026. 

5.  SO 

-0.5 

17-8 

146.96 

1190.44 

1020.2 

029. 

5.8o 

-0.5 

18-1 

275.56 

1179.98 

1026.8 

023. 

inop 

1.5 

■TO  1 

■mu 

1178.08 

1027.1 

021. 

inop 

1.5 

18-3 

20S.39 

1  177.67 

1027.2 

01 S. 

MEDII 

1.5 

18-4 

186.79 

1 1 79.03 

lo27.4 

017. 

3  2.59 

1.5 

18-5 

159.11 

1 185.79 

1027.4 

018. 

13.75 

1.5 

21-1 

252.89 

1282.94 

1023.1 

353. 

inop 

0.8 

21-2 

204.39 

1280.80 

1022.5 

044. 

inop 

0.8 

21-3 

182.12 

1283.45 

1022.6 

035. 

inop 

0.8 

21-4 

173.04 

1283.81 

1022.7 

041. 

inop 

0.8 

21-5 

157.84 

1285.67 

1022.7 

056. 

inop 

0.8 

21-6 

259.15 

1295.27 

1022.S 

083. 

inop 

0.8 

21-7 

207.8S 

1292.45 

1022.8 

089. 

inop 

0.8 

21-8 

ISO.  63 

1291.45 

1022.8 

079. 

inop 

0.8 

21-9 

147.96 

1290.75 

1022.9 

083. 

inop 

0.8 

•  sled  x,y  coordinates  are  referenced  to  the  FRF  coordinate  system 

•  'boa'  is  beach  orientation  angle 

•  MKS  conversion  for  pressure:  1  mb  =  100  Pa 


162 

19.0 

1.59 

.1S75 

0.00 

0.78 

163 

14.9 

i 

.53 

.1875 

0. 1 3 

0.79 

167 


WMBMMBM 

1.46 

.1875 

0.06 

0.58 

165 

wmam 

1.41 

.  1 S75 

(I'M 

O.5o 

166 

9.5 

1.35 

— 

.1S75 

■BB 

o.45 

171 

7.6 

1.16 

.  1S~5 

0.00 

o.5o 

172 

S.8 

1.14 

.1875 

0.0 1 

o.59 

173 

mm 

1.11 

.1719 

0.19 

0.57 

174 

1.8 

1.15 

.1719 

0.19 

mi 

175 

-2.1 

1 .09 

.  1 250 

0.21 

0.38 

176 

mm 

1.09 

.0938 

0.17 

0.24 

177 

-1.9 

1.11 

.07S1 

0.06 

178 

0.2 

1.11 

.0781 

0.07 

0.28 

yifrr)  is  the  mean  a  of  the  rear  and  front  current 
meter  (if  they  were  both  submerged!  and  is  ob¬ 
tained  from  the  significant  wave  ancle  spectra  for 
peak  frequency  in  sea  surface  elevation  spectra 
(see  equation  71 ). 

1I0  is  determined  from  SUPERDICK  gage  630. 

/„,  of  Svx  is  the  peak  frequency  from  'he  rear  cur¬ 
rent  meter  SVJC  spectra 

//rm!  is  determined  from  the  tj  time  series  using 
equation  138 

Q  is  the  fraction  of  all  waves  which  are  breaking. 
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APPENDIX  F.  RADIATION  STRESS,  SEA  SURFACE  ELEVATION, 
SIGNIFICANT  WAVE  ANGLE,  AND  CURRENT  VELOCITY  SPECTRA 
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Figure  65.  Radiation  stress  (SJX),  sea  surface  elevation  (rj).  significant  wave  angle 
(a),  and  current  velocity  spectra  for  sled  run  15-1. 
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Figure  66.  Radiation  stress  (5^),  sea  surface  elevation  (7).  significant  wave  angle 
(£),  and  current  velocity  spectra  for  sled  run  15-2. 
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Figure  67.  Radiation  stress  sea  surface  elevation  (rj),  significant  wave  angle 
(a),  and  current  velocity  spectra  for  sled  run  15-3. 
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Figure  68.  Radiation  stress  (5^,),  sea  surface  elevation  (7),  significant  wave  angle 
(7),  and  current  velocity  spectra  for  sled  run  15-4. 
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Figure  70.  Radiation  stress  ( Syx ),  sea  surface  elevation  (q),  significant  wave  angle 
(a),  and  current  velocity  spectra  for  sled  run  16-3. 
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Figure  71.  Radiation  stress  (5>Jt),  sea  surface  elevation  ( tj ),  significant  wave  angle 
(a),  and  current  velocity  spectra  for  sled  run  16-4. 
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Figure  72.  Radiation  stress  (5>x),  sea  surface  elevation  ( 1 1),  significant  wave  angle 
($),  and  current  velocity  spectra  for  sled  run  16-5. 
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RUN  17-1  —  SUPERDUCK 

95%  x!  (0.73,1.45) 

Sy*  a 
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Figure  74.  Radiation  stress  (5)*),  sea  surface  elevation  (>/),  significant  wave  angle 
(a),  and  current  velocity  spectra  for  sled  run  17-1. 
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Figure  75.  Radiation  stress  ($>.,),  sea  surface  elevation  (»;),  significant  wave  angle 
(£),  and  current  velocity  spectra  for  sled  run  17-2. 
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Figure  76.  Radiation  stress  (S>t),  sea  surface  elevation  (rj),  significant  wave  angle 
(a),  and  current  velocity  spectra  for  sled  run  17-3. 
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Figure  77.  Radiation  stress  (5^,,),  sea  surface  elevation  (q),  significant  wave  angle 
(a),  and  current  velocity  spectra  for  sled  run  17-4. 
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Figure  78.  Radiation  stress  (S>Jt),  sea  surface  elevation  (r/),  significant  wave  angle 
(a),  and  current  velocity  spectra  for  sled  run  17-5. 
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Figure  80.  Radiation  stress  (S>x),  sea  surface  elevation  (tj),  significant  wave  angle 
(5).  and  current  velocity  spectra  for  sled  run  17-7. 
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Radiation  stress  ($„),  sea  surface  elevation  (17),  significant  wave  angle 
(a),  and  current  velocity  spectra  for  sled  run  17-8. 
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Figure  82.  Radiation  stress  (S>x),  sea  surface  elevation  (tj),  significant  wave  angle 
(£),  and  current  velocity  spectra  for  sled  run  18*1. 
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Figure  83.  Radiation  stress  (5>r),  sea  surface  elevation  (ij),  significant  wave  angle 
(a),  and  current  velocity  spectra  for  sled  run  18-2. 
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Figure  84.  Radiation  stress  (Syx),  sea  surface  elevation  (>;),  significant  wave  angle 
(a),  and  current  velocity  spectra  for  sled  run  18-3. 
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Figure  85.  Radiation  stress  (5/x),  sea  surface  elevation  (»/),  significant  wave  angle 
(a),  and  current  velocity  spectra  for  sled  run  18-4. 
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Figure  86.  Radiation  stress  (S>x),  sea  surface  elevation  (tj),  significant  wave  angle 
(£),  and  current  velocity  spectra  for  sled  run  18-5. 
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APPENDIX  G.  MEAN  CURRENT  VELOCITIES  AND  TOTAL 
RADIATION  STRESS  FOR  EACH  SLED  CURRENT  METER  AND  RUN 


Run 

C/M  1  (mid) 

C/M  2  (upper) 

C/M  4  (lower) 

V  (ym-*) 

V 

(ms-1) 

U 

(ms-1) 

V 

(ms-1) 

U 

(ms-1) 

V 

(ms-1) 

U 

(ms-1) 

C/M 

1 

C/M 

4 

C/M 

111 

-0.41  S 

0.367" 

-0.512 

0.266 

-0.354 

0.407 

112 

-0.514 

0.195 

-0.588 

0.155 

-0.405 

0.154 

113 

-0.256 

0.074 

-0.339 

0.130 

-0.164 

0.005 

121 

0.307 

0.252 

0.259 

0.325 

0.304 

0.173 

miiroi 

122 

0.165 

0.517 

0.166 

0.365 

O.Ot.9 

0.592 

s.ss 

13.01 

UMI 

123 

-0.246 

0.295 

-0.320 

0.31S 

-0.1 89 

0.193 

-21.60 

-S.81 

-11.67 

125 

0.295 

0.476 

0.2  SO 

0.471 

0.278 

0.449 

-66.46 

-39.30 

126 

0.219 

0.323 

0.195 

0.2S7 

0.204 

0.270 

-25.27 

-27.02 

127 

0.544 

0.1  IS 

0.49S 

0.476 

0.069 

5. 78 

25.20 

151 

-0.393 

0.06S 

-0.470 

0.10S 

-0.396 

BlUB 

150.00 

1 50.00 

150. 

152 

-0.786 

0.247 

-0.830 

0.141 

-0.798 

0.385 

59.47 

43.67 

7S.87 

153 

-0.773 

0.207 

-O.S82 

-0.045 

-0.723 

0.425 

16.S4 

16.67 

45.94 

154 

-0.719 

0.053 

-0.S24 

0.073 

-0.669 

0. 1 30 

36.62 

22.32 

41.41 

162 

-0.443 

fcixtiifc 

-0.532 

0.052 

-0.437 

0.074 

250.50 

250.50 

250.50 

163 

-0.985 

-1.090 

0.050 

-0.993 

0.076 

153.50 

143.10 

156.40 

164 

-1.209 

fcBBHSi 

-0.931 

-0.320 

-1.156 

0.111 

62.10 

52.29 

£ 

165 

-1.131 

-0.892 

-0.309 

-1.072 

0.036 

4S.40 

47.58 

fr- 

166 

-0.638 

-0.008 

-0.747 

-0.073 

-0.605 

0.070 

27.52 

25.64 

24.66 

171 

-0.219 

-0.024 

-0.301 

0.025 

-0.191 

0.054 

61.53 

61.53 

61.53 

172 

-0.270 

0.015 

-0.355 

0.078 

■39 

53.12 

51.66 

192 


173 

-0.414 

-0.037 

-0.280 

-0.291 

-0.382 

0.074 

25.65 

22.62 

* 

174 

-0.327 

-0.335 

-0.144 

-0.144 

-0.299 

-0.022 

15.34 

* 

175 

-0.299 

-0.485 

-0.114 

-0.041 

-0.376 

-0.150 

* 

10.67 

* 

176 

0.003 

-0.163 

-0.106 

-0.144 

0.023 

0.028 

0.95 

1.95 

177 

0.031 

-0.074 

-0.078 

-0.117 

0.047 

0.055 

4.03 

3.95 

* 

178 

0.038 

-0.045 

-0.091 

-0.133 

0.028 

0.065 

3.56 

1.63 

* 

181 

-0.402 

0.010 

-0.479 

0.060 

-0.376 

0.077 

175.60 

175.60 

175.60 

182 

-0.703 

0.050 

-0.795 

0.125 

-0.693 

0.106 

122.00 

122.60 

129.10 

183 

-0.981 

0.126 

-0.708 

-0.252 

-0.978 

0.227 

58.04 

41.20 

* 

184 

-0.778 

-0.290 

-0.198 

-0.092 

-0.896 

0.073 

* 

14.04 

* 

185 

-0.622 

-0.022 

-0.419 

-0.181 

-0.554 

0.123 

3.02 

4.39 

* 

186 

-0.651 

0.054 

-0.748 

0.075 

-0.656 

0.118 

211 

-0.023 

0.033 

-0.086 

0.107 

-0.012 

0.119 

-6.82 

-6.82 

-6.82 

212 

-0.002 

0.209 

-0.071 

-0.063 

0.010 

0.331 

20.  IS 

27.38 

213 

-0.120 

0.240 

-0.177 

0.218 

-0.098 

0.343 

-22.19 

-16.32 

-14.64 

214 

-0.256 

0.112 

-0.313 

0.096 

-0.276 

0.252 

0.83 

2.71 

2.67 

215 

-0.183 

0.106 

-0.285 

0.176 

-0.121 

0.190 

11.20 

13.86 

9.56 

216 

0.003 

-0.053 

-0.06S 

-0.007 

0.005 

0.064 

4.20 

4.20 

4.20 

217 

0.060 

-0.263 

-0.082 

-0.228 

0.083 

-0.036 

* 

3.44 

* 

218 

-0.150 

-0.037 

-0.204 

-0.140 

-0.156 

0.171 

7.63 

- 

3.83 

6.93 

219 

-0.125 

0.115 

-0.128 

-0.042 

-0.087 

0.173 

-11.35 

-4.10 

1.84 

*  indicates  current  meter  was  not  underwater  for  entire  data  run 
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